AD- A 168  004 


Dr.  H.  W.  Bergmann 

DFVLR  -  Institute  for  Structural  Mechanics 
Braunschweig,  West  Germany 


December  1985 


Final  Report  for  Period  January  1982  -  January  1985 


Approved  for  public  release;  distribution  unlimited. 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

FOR  THE  COMMANDER 


THEODORE  J.  R^HART,  Chief 
Materials  Engineering  Branch 
Systems  Support  Division 


UNCLASW  FI 


.*  1 1 1  i(,A'  ION  ..  U  r  HOfll  T  V 


;3  ►XAIiON.UCWNCHAOINCSCHtDULE 

N/A 


4  jriijANIi.^iON  RE  POR  T  MU  MBE  R  (S) 

None 


3  OlSTRIHUIIONAVAILAHlLITY  Of-  HI  HO H. T 

Approved  for  Public  Release 
Distribution  unlimited 


b.  MONt  roniNG  ORGANISATION  REPORT  NUMBEfllSi 

AFWAL-TR-85-41 50,  Vol  3 


j4  NAME  OP  PERFORMING  ORGANIZATION  |6b.  OFPICE  SYMBOL  I  7a.  NAME  OP  MONITORING  ORGANIZATION 

DFVLR  -  Institute  for  |  ,,t”p“cabl') 


AFWAL/Materials  Laboratory  (AFWAL/MLSE) 


iPaMOir.uiiiM 


5c.  ADDRESS  (City.  Slate  and  Z IP  Code I 

Braunschweig,  West  Germany 


U.  NAME  OF  FUNOING/SPONSORING  8b.  OFI 

ORGANIZATION  (Ifi 

European  Space  Technology  Comm  ttee 


3c.  ADDRESS  iCity.  State  and  ZIP  Code I 


7b.  ADDRESS  ICity.  State  and  ZIP  Code) 

Wright-Patterson  AFB,  OH  45433 


8b.  OFFICE  SYMBOL  9-  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable} 

ttee  N/A  ESTEC/Contract  No  4825/81 /ML/AK(SC) 


10.  SOURCE  OF  FUNDING  NOS. 


,  l»  »  v  •  «  ' 


PROGRAM 

PROJECT 

TASK 

ELEMENT  NO. 

NO. 

NO. 

N/A 

N/A 

N/A 

Brussels,  Belgium  element  no. 

11.  TITLE  (include  Security  Classification)  [)0velOpment  Of  N/A 

Fracture  Mechanics  Maps  for  Composite  Materials 


1J.  PERSONAL  AUTHOR(S) 

Editor  (Dr.  H.  W.  Berqmann 


13a  TYPE  OP  REPORT  13b.  TIME  COVERED  13.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

Final  from  Jan  82  to  Jan.  85.  1985  Dec 


13.  supplementary  notation 

Vol  1  (216) 

_ _ Vol  3  (265) 

17.  COSATI  COOES  V  I  18.  SUBJECT  TERMS  t Continue  on  reverie  if  necessary  and  identii 


15.  PAGE  COUNT 

see  below  (16) 


Vol  2  (204) 

Vol  4  (192 


FIELD 

GROUP 

11 

04 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Fiber  Reinforced  Composites,  Dynamic  Test,  Thermal  Cycling, 
Quality  Control,  Static  Test,  Stress  Concentrations 


9.  ABSTRACT  /Continue  on  reverse  if  necessary  and  identify  by  block  number )  v 

he  Institute  for  Structural  Mechanics  of  the  German  Aerospace  Research  Establishment  (DFVLR 
has  conducted  a  research  program  aimed  at  the  characterization  of  the  fracture  behavior  of 
carbon-fiber-reinforced  epoxy  resins.  The  research  effort  encompassed  a  comprehensive 
experimental  and  analytical  investigation  of  the  response  of  test  specimens  under  a  broad 
range  of  material,  loading  and  environmental  parameters.  By  combining  an  evaluation  of 
global  laminate  properties  with  an  Investigation  of  micro-  and  macroscopic  modes  of  failure, 
the  understanding  of  the  fracture  and  fatigue  behavior  of  carbonfiber-relnforced  composites 
was  enhanced.  It  is  expected  that  the  results  of  the  research  program  and  their  representa¬ 
tion  in  the  form  of  tables,  figures,  graphs  and  maps  will  assist  the  design  engineer  and 
lead  to  improved  engineering  concepts.  >  ;  i  • 


DISTRIBUTION  AVAILABILITY  OF  ABSTRACT 
-NCLASSlFlEO  UNLIMITED  !TX  SAME  AS  APT  C  OTIC  USERS  □ 

::j  name  of  responsible  individual 

THEODORE  J.  REINHART _ 

DD  FORM  1473,  83  APR  EDITION  OF  1  JAN 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

22b  TELEPHONE  NUMBER  22c  OFFICE  SYMBOL 

/Include  In-a  Coda 


1  513/255-3691 

EOITION  OF  1  JAN  73  IS  OBSOLETE. 


SECURITY  CLASSIFICATION  OF  THIS  I’ftlil 


DEVELOPMENT  OF  FRACTURE  MECHANICS  MAPS 


FOR  COMPOSITE  MATERIALS 


Final  Report 
-  Volume  3  - 


APPENDICES  A  -  E 


VOLUME  3 


CONTENTS 


APPENDIX  A 

Charakterisierung  von  Schadensfortschritten  in  CFK-Laminaten 
mittels  Schallemissionsanalyse 


APPENDIX  B 

On-Line  Measurement  of  Onset  and  Growth  of  Edge  Delaminations  in 
CFRP-Laminates  by  an  Optical  Grating  Reflection  Method 


APPENDIX  C 

Eine  effektive  Teilstrukturtechnik  angewendet  auf  Bruchprobleme 


APPENDIX  D 

Crack  Conditions  for  Statically  Loaded  Notched  UD-Laminates 
of  Carbon  Fiber  Reinforced  Plastic.  Part  2: 

Fitting  Curves  for  the  Test  Results 


APPENDIX  E 

Einfluss  von  Temperatur  und  Feuchte  auf  die  Schwingfestigkeit 
ungekerbter  CFK  -  Laminate 


APPENDIX  A 


Charakterisierung  von  Schadensf ortschri tten  in 
CFK  -  Laminaten  mittels  Schal lemi ssionsanalyse 


von 


J.  Block 
DFVLR 

Institut  fur  Strukturmechanik 
Braunschweig 


DEUTSCHE  FORSCHUNGS-  UND  VERSUCHSANSTALT 
FUR  LUFT-  UND  RAUMFAHRT  E.V. 

INST  I  TUT  F(jR  STRUKTURMECHANIK 
Gruppe  Zerstorungsf reie  Priifung 


r.  W*  c  Jr  * 

Sb^t 


Interner  Bericht 
IB  131-84/51 

Charakteri sierung  von  Schadensfortschritten  in 
CFK  -  Laminaten  mittels  Schallemissionsanalyse 


%*£j 


'/■  V. 


.v/.'.v.v, 


Vo- V-yJ 


a.  ■*  -  •*  -*  a 


.■  V  \*>V 


V  .  :**•-.*- a 


Braunschweig,  den  15.11.1984 


Institutslei ter  und  Gruppenleiter : 


Dr.-Ing.  H.W.  Bergmann 


Der  Bericht  umfaJJt  : 

32  Seiten  mit 
14  Bildern  und 
1  Anhang 

Der  Bearbeiter: 

Dipl.-Phys.  J.  Block 


•sw: 

«L* 

.»  -j*  •  *  V 

••vw. 


INHALTSVEkZEICHNMS 


Abstract  . . 

1.  Einfiihrung:  Schallemissionsanalyse  an  CFK-Laminaten 

2.  Schallemissions  -  Messtechnik  . 

3.  Probenlibersicht  . 

4.  Auswertung  der  Versuche  . 


4.1  Detektion  erster  matrixkontrollierter  Schadigungen  13 

4.2  Verfolgung  des  Schadensf ortschri tts  bis  zum  Bruch  17 

4.3  Faserbiindelbruche  und  Probenversagen  .  25 


5.  Zusammenf assung  der  Ergebnisse 

Literatur  . 

Anhang:  Datensatz-Beispiel 


Acceslon  For 

NTIS  CRA&I 
DUG  TAB 
Unannounced 
Justification 

By . . . 

Distribution/ 


Availability  Codes 

Dist 

Avail  and/or 
Special 

fi-l 

1 

f  Q  '  4  f  rjy 


-  M  -f.  •*. 


».  m 


>>>>: 

w-V-VaW 


IK 


jSSi&x 


v.>vv.; 

V 

"  v' 


.V, 


vtt'S'SS 

•  v’VV.Vd 


ABSTRACT 


Dieser  Bericht  enthalt  eine  Zusammenf assung  und  Deutung  der  Er- 
gebnisse  von  Schallemissionsversuchen  an  CFK-Laminaten,  die  sta- 
tisch  oder  quasi stati sch  im  Zugbereich  durchgefiihrt  wurden.  Er 
beschrankt  sich  auf  die  Auswertung  der  MeBergebnisse  an  Probesta- 
ben  aus  einem  Material  (Fasern:  T  300,  Harz:  914  C) ,  aber  mi  t 
variierender  Lagenfolge.  Der  Schwerpunkt  wurde  zum  einen  auf  das 
Einsetzen  der  Schal lemi ssionen,  die  den  feststellbaren  Beginn 
matrixkontrollierter  Schadigungen  bezeichnen,  zum  anderen  auf 
den  Bereich  vor  und  wahrend  des  Versagensbeginns  gelegt.  Das 
Fortrchrei ten  des  Schadi gungsprozesses  wurde  anhand  der  gemesse- 
nen  Schal lemi ssionen  ebenfalls  interpretiert .  Hier  sind  die 
i rreversiblen  Schadigungsverlauf e  der  einzelnen  Proben  individu- 
ell  exakt  verfolgbar,  laminatspezi f i sche  Charakteri stika  dagegen 
schwacher  ausgepragt. 


In  this  report,  acoustic  emission  test  results  on  CFRP  laminates 
under  static  or  guasi-static  tensile  load  are  summarized  and  in¬ 
terpreted.  The  scope  of  this  evaluation  was  restricted  to  slender 
test  specimens  of  T  300  fibers  and  914  C  epoxy  resin  with  varying 
stacking  sequences .  Emphasis  was  placed  upon  the  AE  initiation 
which  indicated  the  onset  of  resin-controlled  damage  mechanisms 
and  also  upon  the  emissions  occurring  shortly  before  and  during 
failure.  The  irreversible  process  of  damage  progression  with  in¬ 
creasing  load  could  be  investigated  in  detail  for  each  individual 
specimen,  whereas  the  repeatability  of  the  results  for  a  group  of 
specimens  of  the  same  type  was  less  exact. 


3 


einsetzen,  wenn  dieser  belastet,  also  letztlich  doch  geschadi .  tr 
wird.  Allerdings  braucht  die  Testbel astung  nicht  kritisch  ou 


Der  Nutzen  der  hohen  Auflosung  hangt  primar  von  der  Frage  ab,  >; 
wieweit  die  einzelnen  Signale  bestimmten  Schadensmechani  snv’n 
zugeordnet  werden  konnen.Die  Probleme  liegen  dabei  in  der  Benic':- 
sichtigung  der  komplizier  ten  Schallubertragungseigenschaf r  . 
des  Materials  CFK  ( hochgradige  Ani  sotropie ,  Dispersion,  Dampfu..  : 
usw .  )  .  Bis  jetr.t  ist  die  Unterscheidung  zwiscU 
Facer  (biindel  )bruchen  und  matrixkontrollierten  Mechanismen  •  • 
ei  nr.«lne  Scha  1 1  ere i  gni  sse  veri  f  izierbar  .  Eine  wei  tergehc-- 
Schadenskl  assi  f  i  kation  kann  aufgrund  makroskopi  seller  SE-Verb 
lung on  vorgenommen  werden. 

I  m  vorliegenden  Bericht  werden  derartige  SE-Vertei lungen 
Prob^rtaben  aus  verschiedenen  Laminaten,  aber  dem  gleichen  Ma  • 
rial,  ausgewertet.  Die  Auswertung  der  Haufigkeit  von  Scha  l '  •• 
ereigni  ???n  und  der  Werte  ihrer  gemessenen  Parameter,  vor  a  1  ’ 
Erm;  u  i evertei lungen,  zeigen  Korrelationen  zu  Materi a  lei  a-  •- 
schaften  und  Ee  1  astungsar t  auf.  Dabei  lassen  sich  Aussagen,  d  •  ■ 
fir  das  ieweilige  Lamina t  bzw.  die  Belastungsf orm  typisch  sin- 
von  yolchen  trennen,  die  die  indiv iduel len  Eigenschaf ten  d-  ?■ 
rriif lings  beschreiben. 

I n.-hesondere  eignet  sich  die  SEA  zum  Nachweis  erster  Schadigun  ■  •• 
w.i  e  den  plbt.zlichen  Beginn  des  Einspringens  von  Rissen  in  • 
Harz.  Mechanische  Auf  nehmer  (z.B.  Weqaufnehmer )  nehmen  d- 
F  i  sse  noch  nicht  wahr;  in  der  Energi evertei  lung  der  Scha  1 1  em  i  -•  - 
sicnen  zeichnen  sie  sich  dagegen  deutlich  ab.  Dasselbe  gilt  fi-> 
Schadicmncren ,  die  bei  konstanter  Belastung  durch  viskoelnr  ■- 
scdies  Kriechen  in  der  Probe  entstehen.  Auch  die  dem  Bruch  i:  ■ 
Frebe  mrangehende  Phase  kann  mittels  SEA  qenauer  analysiort  v> 
den.  Hier  spielen  bei  bestimmten  Laminaten  Faserbruche  e 
beherrschende  Rolle.  Diese  setzen  akustische  Enerqien  frei,  d 
sehr  vi  el  groBer  sir.d  als  bei  Kohasionsbruch  in  der  Mat- 

(Rinse,  Delaminationen)  Oder  Adhasi onsbruch  zwischen  Ease: . 

Matrix  (Debonding). 


*.  \  ■-  *.  \ , 


►  •  *  '  k.  '  »  k  .  .  , 
v"  O  V  *  *  ' 


der  Korrelation  mit  mechanischen  Phanomenen  und  damit  die  gewinn 
bare  Inf ormationsmenge  entscheidend  ab. 
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Bild  2.  Frequenzbereich  der  Schallemissionsanalyse 

Bild  2  zeigt  schematisch  den  Frequenzbereich  der  akustischen 
Emissionen  von  Festkorpern.  Er  ist  auBerordentlich  breit  und  be- 
ginnt  im  horbaren  Bereich,  wobei  die  Abgrenzung  gegeniiber  den 
noch  tieferen  Frequenzen  mechanischer  Schwingungen  des  Bauteils, 
die  man  nicht  mehr  als  Schall  bezeichnen  will,  mehr  oder  weniger 
willkurlich  bleibt.  Auch  nach  oben  gibt  es  keine  feste  Grenze, 
aber  extrem  hohe  Frequenzen  werden  normalerwei se  nur  von  einzel- 
nen  Schal lquanten  (Phononen),  nicht  jedoch  von  makroskopischen 
Schallwellen  erreicht.  Fur  die  Schallemissionsanalyse  an  CFK  ist 
ein  mittlerer  Bereich  von  etwa  100  kHz  bis  1-2  MHz  interessant. 
Unterhalb  davon  ist  der  Anteil  der  Storgerausche  aus  der  Umwelt 
recht  hoch;  bei  wesentlich  hoheren  Frequenzen  nimmt  die  Signalnm- 
plitude  rasch  ab.  Die  Eingrenzung  dieses  Bandes  erfolgt  mittels 
geeigneter  Filter.  Die  in  diesem  Bericht  diskutierten  Versuche 
wurden  mit  Sensoren  gemacht,  deren  Resonanzkurve  bei  150  kHz  ein 
Maximum  besitzt.  Fur  die  Frequenzanalyse  muB  die  Aufnahme  aller- 
dings  breitbandig  erfolgen. 

Bei  Angabe  der  Frequenz  ist  stets  das  Intensi tatsmaximum  der  Fre- 
quenzvertei lung  eines  Schallsignals  gemeint,  die  immer  eine  er- 
hebliche  Breite  aufweist.  Auch  besteht  das  Schallsignal  im  allge- 
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meinen  aus  mehreren  Wellenmoden  mit  jeweils  anderer  Dispersions- 
relation  usw. ,  so  daB  eine  analytische  Behandlung  der  Wellenaus- 
breitung  hier  nicht  in  Frage  kommt.  Jedoch  sind  die  durch  die 
starke  Anisotropie  des  Materials  CFK  hervorgeruf enen  Effekte  bei 
der  Interpretation  der  Ergebnisse  zu  berucksichtigen. 

Die  Schallemissionstechnik  ist  in  den  letzten  Jahren  entschei- 
dend  verbessert  worden.  Wahrend  Gerate  alterer  Bauart  meist  nur 
die  Ereignissummen  und  -raten  registrieren  oder  lediglich  die 
Uberschreitungen  einer  festgelegten  Schwelle  durch  die  Signal- 
spannung  zahlen  (sog.  counts),  tasten  modernere  Anlagen  die  aku- 
stischen  Signale  beziiglich  mehrerer  Parameter  ab,  wie  oben 
beschrieben.  Eine  groBe  Zahl  von  Versuchen  ist  von  Awerbuch  unter 
Assistenz  des  Verfassers  an  einer  Anlage  vom  Typ 

Dunegan  -  Endevco  DE  3000 

durchgefuhrt  worden  [1-11].  Fur  wei tergehende  Untersuchungen 
steht  seit  1982  die  in  Bild  3  dargestellte  mikrocomputergestiitzte 
Anlage 

Physical  Acoustics  PAC  3000/3004 

zur  Verfiigung.  Alle  Tests  mit  hoher  Auflosung  und  alle  Energie- 
messungen  wurden  mit  dieser  Anlage  durchgefuhrt. 

Die  mit  den  piezoelektrischen  Sensoren  (vgl.  dazu  z.B.[12])  auf- 
genomraenen  Schal lsi gnale  werden  rauscharm  verstarkt  und  dann 
beziiglich  der  in  Bild  1  dargestel  lten  Parameter  abgetastet: 

( Spitzen- ) Amplitude,  Ereigni sdauer ,  Counts  (Zahl  der  Schwellen- 
iiberschreitungen) ,  Anstiegszeit  und  Energieinhalt .  Die  Energie 
wird  durch  Abtastung  der  Signalform  bestimmt.  Die  Ereignisse  wer¬ 
den  numeriert,  auf  der  Zeitskala  datiert  und  mit  gleichzeitig 
gemessenen  Werten  mehrerer  externer  Parameter  korreliert.  AuBer- 
dem  werden  die  Lauf zeitdif ferenzen  eines  Signals  zu  den 
verschiedenen  Aufnehmern  (4  parallele  Kanale)  registriert. 

Die  Daten  werden  digitalisiert  und  auf  einen  Datenbus  gegeben, 
von  wo  sie  in  der  Regel  unmittelbar  auf  Floppydisks  geschrieben 
werden.  Von  dort  konnen  sie  zwecks  Aufbereitung  stets  wieder  ab- 
gerufen  und  zu  Listen,  grafischen  Verteilungen  usw.  weiterverar- 
beitet  werden.  Eine  begrenzte  Anzahl  von  Diagrammen  kann  man  auch 
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Bild  3.  Schallemissionsanlage  PAC  3000/3004 

on-line  wahrend  des  Versuchs  auf  dem  Bildschirm  verfolgen.  Zu- 
satzlich  konnen  noch  Frequenzinhalt  und  Signalform  der 
Eingangssignale  analysiert  werden,  wenn  diese  mit  dem  Transien- 
tenrekordereinschub  auf gezeichnet  und  mittels  FFT-Algori thmen 
(schnelle  Fouriertransf ormation)  ausgewertet  werden. 

Vor  der  Diskussion  von  Versuchsergebni ssen  muB  das  Problem  der 
eindeutigen  Korrelation  zwischen  Schal lausbrtichen  (bursts)  im 
Werkstoff  und  den  regi strierten  Ereignissen  diskutiert  werden. 
Es  ist  zu  fordern,  daB  jedes  registrierte  Schallereignis  genau 
einem  verursachenden  "MikroprozeB" ,  z.B.  einem  Faserbruch,  zuge- 
ordnet  werden  kann,  wobei  die  richtige  Interpretation  des 
Ereignisses  dann  Sache  der  weiteren  Analyse  ist. 

Die  Eindeutigkeit  dieser  Zuordnung  kann  bei  hohen  Datenraten  emp- 
findlich  gestort  sein,  wenn  namlich  die  Wellenpakete  von  ver- 
schiedenen  Quellen  sich  zeitlich  derart  iiberlappen,  daB  das  Sy¬ 
stem  sie  als  1  Ereignis  auffaBt.  Die  SE-Anlage  kennt  beztiglich 
der  Schallereignisse  nur  die  zwei  Zustande  "Ereignis"  und  "kein 
Ereignis",  die  formal  als  Zustandsfunktion  Z(t)  mit  den  moglichen 
Werten  1  und  0  bezeichnet  werden  soli  (siehe  Bild  4).  Uberschrei- 
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Bild  4.  Zur  Eindeutigkei t  der  Schallereignismessung 

tet  die  Signalspannung  U  die  eingestellte  Schwelle  (wird  also  ein 
Count  erzeugt),  so  wird  Z(t)=l  und  fallt  erst  dann  wieder  auf  0 
zuriick,  wenn  nach  dem  letzten  Count  noch  eine  gewisse  vorher  ein¬ 
gestellte  Totzeit  0  verstrichen  ist.  Liegen  nun  die  Emissionen 
sehr  dicht,  so  bleibt  fast  standig  Z(t)=l,  d.h.  die  Zahlrate  der 
Ereignisse  nimmt  stark  ab,  die  gemessene  Ereignisdauer  wachst 
sprunghaft  an. 

Betrachtet  man  ein  willkiirlich  aus  einer  Folge  von  N  Ereignissen 
herausgegrif f enes  (Nummer  n),  so  wird  dies  richtig  registriert 
genau  dann,  wenn  die  durch  die  anderen  Ereignisse  definierte  Zu- 
standsfunktion  an  seinem  Anfang  und  an  seinem  Ende  gleich  0  ist 
( Wahrscheinlichkeit  dafiir  sei  jeweils  p)  und  wenn  kein  anderes 
Ereignis  vollstandig  durch  das  n-te  uberdeckt  wird  (Wahrschein¬ 
lichkeit  q) .  Solange  die  N  Ereignisse  iiber  das  Zeitintervall,  in 
dem  sie  auftreten,  statistisch  gleichverteilt  sind,  ergibt  das 
Produkt  der  o.a.  Wahrscheinlichkeiten  den  Erwartungswert  $  fur 
eine  richtige  Korrelation,  also  fur  eine  korrekte  Ereigniszah- 
lung. 


Die  Wahrscheinlichkeit  p  kann  man  aus  der  im  betrachteten  Zeitin- 
tervall  T  gemessenen  Summe  der  Ereignisdauern  bestimmen.  Es 
ist  namlich 

p--  i-y(2d'*2d0  s1-f  2Di  IN:&1) 

i=1  i  =  n+1  i  =  1 

0  =  p2q  =  q  (1  (ED')-y  ZDit1) 

i:1  1=1 

(zu  den  Bezeichnungen  vgl.  Bild  4).  Die  Wahrscheinlichkeit  q  laBt 
sich  nur  mit  relativ  hohem  Rechenaufwand  abschatzen.  Da  sie  aber 
per  def.  hochstens  gleich  1  ist  und  $  dann  ein  Maximum  annimmt, 
ist  es  sinnvoll,  als  notwendiges  (aber  i.allg.  nicht  hinreichen- 
des)  Kriterium  fur  die  Korrektheit  der  Messung  zu  fordern,  daJ3 


^max  =  ^q-1  ~  1 


Die  Kontrolle  der  Korrelation  der  registrierten  Schallereigni sse 
mit  den  v^rursachenden  diskreten  Schadensereigni ssen  ist  unver- 
zichtbar.  Bei  der  Durchfuhrung  der  in  den  Berichten  [1-9]  doku- 
mentierten  Versuche  standen  entsprechende  Programme  noch  nicht 
zur  Verfugung.  Der  groBte  Teil  der  dort  aufgezeigten  Schallereig¬ 
ni  s-Summenkurven  ist  im  Bereich  kurz  vor  dem  Bruch  der  jeweiligen 
Probe  verfalscht,  wie  Nachmessungen  ergeben  haben.  Dort  ist  deut- 
lich  <  1  . 

Bei  den  mit  der  PAC  3000/3004  -  Anlage  durchgefiihrten  Messungen 
der  SE-Energie  konnte  dieser  Fehler  nicht  mehr  auftreten,  da  die 
im  Intervall  T  freigesetzte  Energiesumme 

N  t, ♦  D, 

e  =  £  (  /ultldt) 

1=1  I  =  t, 


von  der  Zahlung  der  Ereignisdauern  D^^  unabhangig  ist. 


-  10  - 

Die  lineare  Ortung  der  Schallquellen  langs  der  Probestablange  L 
kann  durch  Messung  der  Lauf zei tdi f f erenz  t^-t^  der  Schallsignale 
beziiglich  der  beiden  Sensoren  A  und  B  erfolgen  (vgl.  [13]).  Die 
Ortslinien  gleicher  Laufzeitdif ferenz  sind  bei  isotroper  Schall- 
geschwindigkeit  c  Hyperbeln  (siehe  Bild  5),  die  die  Probenlangs- 
achse  in  je  einem  Punkt  schneiden.  Dieser  Schnittpunkt  bleibt 
erhalten,  wenn  die  Ortlinien  wegen  der  tatsachlichen  starken  An- 
isotropie  der  Schallgeschwindigkeit  bei  CFK-Laminaten  seitlich 
auch  von  der  Hyperbelform  abweichen,  da  c  langs  der  Achse  AB  als 
konstant  angenommen  werden  kann. 


Sensor  A  Ortung  t  Sensor  8 


Bild  5.  Zur  linearen  Ortung  von  Schallereignissen 

Der  Ortungswert  l  in  Prozent  der  Lange  L  zwischen  den  Sensoren  A 
und  B  ergibt  sich  zu 


Bild  6  verdeutlicht  den  EinfluB  der  Dampfung:  Bei  angenommener 
Gleichvertei lung  der  Ereignisse  uber  der  Probenlange  nimmt  die 
Zahl  der  registrierten  Ereignisse  mit  zunehmender  Entfernung  x 
von  der  Probenmitte  ab,  so  daB  sich  statt  einer  horizontalen  eine 


durchgebogene  Vertei lungskurve  ergibt.  Bei  [ ±45° ] -Proben  wurden 
derartige  Vertei lungen,  wie  zu  erwarten,  auch  tatsachlich  gemes- 
sen.  Dieser  Effekt  hangt  von  der  Vertei lungsfunktion  f(U)  der 
"ungedampften"  Amplituden  U  sowie  vom  Betrag  der  Dampfung  langs 
L/2  +  x,  L/2  -  x  ab,  welche  die  Lage  der  Regi strierschwel len 
SA  B(x)  beziiglich  f(U)  bestimmt.  Nur  der  Bereich  U>SA(x)  in 
Bild  6  entspricht  den  tatsachlich  georteten  Schallereignissen. 


Mefllange  L 


I  ! 


!  f  ffuldu 

1  N  (x)  =  N  (0)  — - 


/  fluldu 

5(0) 


SABtx)  =  Schwellwert 
fur  Sensor  A,B 
in  Abhongigkeil  von 
54|x)  der  Dampfung 

k  fangs  L*x.  L-x 
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3.  PROBENUBERSICHT 


Es  wurden  Versuche  an  Probestaben  aus  folgenden  Laminaten1)  aus- 
gewertet : 

[0°]g  ungekerbt  und  gekerbt 

statische  Zugversuche 
langsame  zyklische  Belastung 

[0]g  off-axis-Winkel  0=2 . 5“ . . . 20° ; 

ungekerbt 

statische  Zugversuche 


[±45° ]. 


ro° ! so0 ]2s 


[02|M5 1021-4510190], 


ungekerbt 

statische  Zugversuche 
langsame  zyklische  Belastung 

und  Variationen  der  Lagenfolge;  unge¬ 
kerbt 

statische  Zugversuche 
langsame  zyklische  Belastung 

gekerbt  und  ungekerbt 
statische  Zugversuche 
langsame  zyklische  Belastung 


Die  Laminate  waren  aufgebaut  aus: 


Fasern  : 


Harz 


Torayca  T  300,  hochfest 


Ciba  914  C 


Zusatzlich  wurden  zum  Vergleich  einige  Reinharzproben  getestet. 


l)  Die  Laminate  halten  sich  in  dem  im  Proposal  des  Instituts  zum 
ESTEC/Contract  No . 4825/81/NL/AK  (SC)  angegebenen  Rahmen 
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4.  AUSWERTUNG  DER  VERSUCHE 


4.1  Detektion  erster  matrixkontrollierter  Schadigungen 

Die  ersten  Schadigungen  in  CFK-Laminaten  sind  generell  matrix- 
kontrolliert .  ( Faserversagen  tritt  unter  reiner  Zugbeanspruchung 
erst  bei  hoheren  Lasten  bzw.  unter  schwingender  Belastung  erst 
nach  vorherigem  pull-out  als  Knickbruch  auf).  Also  sind  Koha- 
sionsbriiche  im  Harz  (Mikrorisse)  oder  Adhasionsbruche  zwischen 
Harz  und  Faser  die  ersten  anhand  akustischer  Emissionen  fest- 
stellbaren  Schadigungen. 

Bei  Reinharzproben  wurden  diese  Emissionen  erst  viel  naher  der 
Versagensbeanspruchung  gemessen  als  bei  CFK-Laminaten:  die  im 
unverstarkten  Zustand  relativ  groBe  Bruchdehnung  der  Harze  kann 
im  Laminat  nicht  erhalten  bleiben,  da  durch  das  Vorhandensein  der 
Fasern  stets  ein  komplizierter  Spannungszustand  entsteht  und  die 
Duktilitat  der  Probe  grundsatzlich  vermindert  wird.  Dadurch 
komrnt  es  schon  bei  geringer  Belastung  des  Probestabes  zu  ortli- 
chetn  Versagen.  Ein  Bruchteil  der  dabei  f reigesetzten  Energie  wird 
in  die  Schallemi ssionssensoren  dissipiert  (zur  quanti tativen 
Abschatzung  siehe  Anhang). 

In  den  von  Awerbuch  und  d.Verf.  durchgefiihrten  Versuchen  war  je- 
weils  ein  AE  initiation- Punkt  gefunden  worden.  In  einigen  Kon- 
trollversuchen  wurde  nun  zunachst  iiberpruft,  ob  sich  diese  Punkte 
mit  der  verbesserten  PAC-Anlage  verifizieren  lieBen.  Dies  war  der 
Fall  auBer  bei  den  [0°|900]-  und  den  [ ±45° ) -Proben.  Die  in  Kap.3 
diskutierten  und  bei  der  alten  DE-3000-Anlage  auf getretenen  Kor- 
relationsfehler  in  der  Ereigni szahlung  spielen  bei  den  ersten 
Schadigungen  wegen  der  niedrigen  Datenrate  noch  keine  Rolle.  Da- 
her  konnten  die  damit  erzielten  MeBergebni sse  fur  die  Auswertung 
hinsichtlich  erster  Schadigungen  noch  herangezogen  werden: 


Nach  [1,  4,  5]  liegt  der  Einsatz  erster  Schallemissionen  bei  uni- 
di rektionalen  Probestaben  wie  folgt  (Werte  gerundet): 


Finit  Fult  ^init^^ult 


!°0]8 

Innenkerbe 

12 . 5mm 

2 

kN 

1°°  )8 

11 

1! 

7 . 5mm 

3 

kN 

[ 0  ° !  8 

11 

11 

2 . 5mm 

5.5 

kN 

[°’]8 

ungekerbt 

13.5 

kN 

40 

kN 

0.34 

t  50  }8 

off- 

axis 

12 

kN 

21 

kN 

0.57 

[7.5®] 

11 

B 

ti 

10 

kN 

16 

kN 

0.63 

[1°"]8 

n 

11 

9 

kN 

12.5 

kN 

0.72 

(15  'U 

rr 

n 

7.2 

kN 

8.3 

kN 

0.86 

)  00 

0 

O 

ii 

11 

5.4 

kN 

5.8 

kN 

0.93 

Da  bei  Versuchen  mit  [  0®  |  90°  ] -Proben  bei  unterschiedlichen  Tempe- 
raturen  und  Feuchtigkei ten  (Proben  576/1  ff .  in  den  Tests  AE012  - 
AE021)  Abweichungen  beim  Einsetzen  erster  Emissionen  gegeniiber 
den  in  [3]  dokumentierten  Messungen  gefunden  worden  waren,  wurde 
eine  entsprechende  MeBreihe  mit  hoher  Auflosung  gefahren: 


Test 

Probennr . 

Finit 

Fult 

Fini  t^ult 

AE217 

942/1 

2.58 

18.41 

0.14 

AE218 

942/2 

3.90 

19.00 

0.21 

AE219 

942/3 

2.62 

19.29 

0.13 

AE220 

942/4 

3.04 

18.67 

0.16 

AE221 

942/5 

2.97 

18.92 

0. 16 

AE222 

942/6 

3.40 

18.61 

0. 18 

AE223 

942/7 

2.86 

19.05 

* 

0.15 

AE224 

942/8 

5.2  (  ! 

•  ) 

19.54 

★ 

0.26  ( ! ) 

AE225 

942/9 

3.61 

19.20 

* 

0.19 

AE226 

942/10 

3.20 

19.35 

* 

0.17 

AE227 

942/11 

nicht 

auswertbar  wegen 

MeBf ehler 

AE228 

942/12 

2.77 

18.85 

* 

0.15 

AE229 

942/13 

3.45 

18.48 

* 

0.19 

Die  mit  *  bezeichneten  Werte  sind  die  Bruchlasten  bei 
Feiicit y-Versuchen  (langsame  zyklische  Belastung) 

Aus  dem  Standardlaminat  | O2 I +45 | O2 I -45 1 0 1  90 ] g  gefertigte  Proben 
mit  Innenkerbe  zeigen  Schallemissionen  ab  Belastungsbeginn;  ohne 


Amplitude 


[0/90], c  942/1 


.v. ... 


%  . 


■  .  lb  ft 


rel.  Last  F/F..I* 


Bild  7.  Erste  Schallereignisse  bei  quasi stati scher  Belastung 

Kerbe  treten  die  ersten  Emissionen  bei  F/Fult“0-3  auf  [6],  Bei 
|  ±45° j -Proben  ist  die  Festlegung  des  Schadigungsbeginns  anhand 
erster  Schallemissionen  fragwiirdig.  Die  in  [2]  dargestellten 
Werte  ( Emi ssionsbeginn  bei  3  -  3.5  kN,  Bruchlast  =  5  kN)  beruhen 
auf  der  hier  problemati schen  Wahl  der  Di skrimi natorschwel le .  Die 
Faser lagen  verdrehen  sich  aber  von  Anfang  an  gegeneinanden  und 
erzeugen  dabei  Schallemissionen  wachsender  Amplitude  (Eichver- 
such  AE 195-4 ,  270/10:  kein  scharfer  Ubergang  zwi schen 


Untergrundgerausch  und  diskreten  Emissionen  von  "echten"  Schadi¬ 
gungen)  . 


Der  Ort  der  ersten  Schadigungen  entspricht  bei  Probestaben  mit 
Innenkerbe  (bei  50%  der  Stablange)  im  wesentlichen  £=50.  Bei  un- 
gekerbten  Proben  ist  keine  Systematik  erkennbar,  wie  folgende 
Aufstellung  der  jeweils  ersten  10  Ereignisse  von  5  Tests  zeigt: 


Test-Probe  (Lam)  Ortung  ^  ( i=l , . . . , 10)  £  s 


;  i ingekerbt  | 

AE08 1 -n . gek . 

(1) 

04 

04 

29 

78 

35 

80 

83 

56 

66 

17 

45.2 

31.3 

AE2 13-596/1 

(2) 

90 

32 

58 

33 

26 

24 

17 

89 

26 

69 

46.4 

27.8 

AE223-942/7 

(3) 

16 

30 

41 

67 

99 

99 

43 

23 

76 

70 

56.4 

30. 1 

*  gekerbt  | 

AE23 1-263/1 1 

(1) 

52 

55 

54 

47 

55 

58 

20 

54 

55 

61 

51 . 1 

11.5 

AE230-274/14 

(2) 

51 

48 

48 

65 

49 

51 

47 

49 

49 

48 

50.5 

5.3 

Laminat  ( 1 )=[ 

0°  ) 8 

(2) 

=Standardlam 

• 

(3) 

=  [0 

•|90’]2s 

Hierbei  ist  s  =  +\/s*die  Standardabweichung,  die  aus  der  Varianz 

10 

S’  =  k 

i  =  l 

bestimmt  wurde.  (Ohne  den  "AusreiBer"  £^=20  bei  AE231  ist  dort  die 
Standardabweichung  ebenfalls  s=5 . 3  wie  bei  AE230). 

Die  Bestimmung  des  Auftretens  erster  Schadigungen  mittels  SEA  ist 
mitzlich  vor  allem,  weil  sich  diese  bei  Laminaten  aus  Endlos  - 
Kohlefasern  nicht  durch  deutliche  Anderungen  des 
Spannungs/Dehnungs  -  Verlaufs  im  Sinne  des  "Puck'schen  Knies" 
(vgl.  ( 14]  )  bemerkbar  machen.  Eine  quantitative  Ruckrechnung  von 
den  auf genommenen  Schallwellen  auf  den  verursachenden  Schadi- 
gungsprozeB  (RiBspitze)  wird  wegen  der  komplizierten 
Werkstoff struktur  noch  nicht  beherrscht,  obwohl  vergleichbare 
Ansatze  fur  isotrope  We rks to ffe  existieren  (z.B.  [15]). 


4.2  Verfolgung  des  Schadensfortschritts  bis  zum  Bruch 

Unter  linear  steigender  Zugbelastung  wachsen  die  Schadigungen  im 
Innern  der  Probe  rasch  an.  Durch  Messung  insbesondere  der  in  den 
Schallemissionen  enthaltenen  Energien  kann  die  insgesamt  freige- 
setzte  elastische  Energie  und  damit  das  AusmaB  des  Schadenswachs- 
tums  beurteilt  werden. 


Ein  Beispiel  zeigt  das  folgende  Bild: 


I  AE  081 


Bruch  t 


Bild  8.  Schallemissionsenergien  einer  [0  ]g-Probe  unter  Zugbe¬ 
lastung  (in  log.  Energy-Counts ) 


Verteilungen  wie  die  in  Bild  8  dargestellte  konnen  noch  um  Gi  b- 
Benordnungen  genauer  aufgelost  und  ausschnittwei se  dargestellt 
werden.  Ihr  Verlauf  ist  auch  bei  gleichartigen  Proben  unter  iden- 
tischen  Versuchsbedingungen  nicht  reproduzierbar  gleichf ormi  q, 
was  auf  den  starken  EinfluB  individueller  Ei genschaf ten  der  e  i  n- 
zelnen  Proben  bzw.  die  I rreversibi li tat  des  Schadigungsprozesres 
hinweist.  Ein  alle  Laminattypen  iibergrei  f  endes  Charakteri  sti  hum 
ist  al  lerdings  die  im  unteren  Lastbereich  (F^F-  F  ,  .  )  hauf  q 
exponentiel le  Grundform  der  Ereignis-  und  Energiesummenkurven . 
Dies  ist  konsistent  mit  der  Annahme,  dafl  das  Fortschrei ten  d«r 
Schaden  und  damit  der  Schallemissionen  hier  vor  allem  vom  AusnnB 
der  schon  vorhandenen  Defekte  abhangt: 


I  I 

.80 


2.80  4.08  6.00  8.00 


Bild  10.  SE-Energiesunune  einer  gekerbten  (0#)R-Probe 


Je  nach  Laminat  setzen  die  zunachst  angenahert  exponentiell  ver- 
laufenden  Energiesummenkurven  bei  unterschiedlich  hohen  Bela- 
stungen  ein  (vgl.4.1).  Mit  zunehmender  Last  wird  die  Rate  der 
f reigesetzten  akustischen  Energien  immer  unregelmaBiger,  wenn 
auch  global  grolier .  Dies  muB  als  typisch  fur  das  Material  CFK  hin- 
genommen  werden.  Die  Verfolgung  des  Schadensf ortschritts  mit 
hoher  Auflosung  ist  bei  Kurzzeit  -  Zugversuchen  nur  im  jeweiligen 
Einzelfall  sinnvoll. 

Wesentlich  einheitlicher  laBt  sich  die  Zerrtittung  der  Proben  bei 
langsamer  Ermtidungsbelastung  charakteri  sieren .  Gegeniiber  der 
linear  steigenden  Zugbelastung  werden  zwei  weitere  Mechanismen 
wichtig : 

Der  Dehnungsverlauf  c(o,t)  von  Faserverbunden  wird  jenseits  des 
elastischen  Bereichs  von  zwei  e -Komponenten  bestimmt,  die  zum 
vi skoelasti schen  Kriechen  bzw.  zum  viskoplasti schen  FlieBen  fiih- 
ren  (vgl.  z.B.  [16]).  Zeitproportional  fortschrei tendes  FlieBen 
spielt  auBer  bei  Reinharzen  keine  Rolle  in  Kurzzeitversuchen. 
Wenn  jedoch  o(t)  ein  relatives  Maximum  durchlauft  oder  festgehal- 
ten  wird,  verursacht  das  vi skoelasti sche  Kriechen  des  Werkstoffs 
in  eine  asymptotisch  angestrebte  Ruhelage  weiteres  Schadens- 
wachstum  und  damit  Schal lemi ssionen .  Mit  zunehmender  Zerrtittung 
tragt  RiBuf erreibung  als  zweiter  Mechanismus  zu  den  Schallemis- 
sionen  unterhalb  der  jeweils  erreichten  Maximallast  bei. 

Das  in  Bild  13  erklarte  Felicity-Verhaltnis  R  (Index  z=Zyklus) 

Rz  _  FAE,  z^max,  z  “  °AE,  z^max,  z  "  * 

ist  ein  phanomenologisches  MaB  ftir  den  Schadigungs-  bzw.  Zerrtit- 
tungszustand  der  Probe.  Awerbuch  fand  ftir  Belastungsverlauf e  mit 
wachsendem  F_,„  _  eine  monotone  Abnahme  von  R=1  bis  R=0.8  ftir  alle 
getesteten  Laminate  und  stellte  eine  Abhangigkeit  der  Schall- 
ereignisrate  pro  Zyklus  von  der  erreichten  Maximallast  F__„ 

iTldX  /  Z 

bzw.  -spannung  o  und  der  Steigerung  gegeniiber  Fmo„  ,  fest 
[  10]  . 

Urn  diesen  Zusammenhang  interpretieren  zu  konnen,  wurden  Versuche 
zur  Verfolgung  der  im  z-ten  Zyklus  f reigesetzten  Energie  Ez  in 
Abhangigkeit  von  Fmax  z  und  z/Fmax  ^  durchgef uhrt .  Die 


Bild  13.  Felicity  -  Effekt 


Energie  ist  aus  in  Kap.2  diskutierten  Griinden  ein  geeigneterer 
Parameter  als  die  Schal lereigni srate . 

Da  diese  Messungen  zunachst  unsystemati sell  durchgefuhrt  wurden, 
urn  geeignete  Versuchsbedingungen  zu  erproben,  steht  bisher  erst 
ein  Satz  von  4  zusammenhangenden  Tests  ( AE224 , 225 , 226 , 228 )  fur 
die  rechnerische  Auswertung  zur  Verfugung.  Diese  stimmen  beziig- 
lich  Probentyp  (  [0°|90°  J^s'  / 5 , 10, 12 ) ,  Belastungsrate  ( In- 
stron  Priif geschwindigkeit  5  mm/min)  und  AE-Versuchsparametern 
uberein.  Das  Verhaltnis  „/F„-v  r,  i  wurde  variiert  (siehe 
folgende  Seite). 

Es  stellte  sich  heraus,  daB  die  Energierate  pro  Belastungszyklus 
ohne  statistisch  verif izierbare  Tendenz  schwankt,  solange 

<Fnax,z-Fmax.z-l»/fmax,z4l  ist-  Der  der  naeh  30  Zyklan 
gemittelten  Differenz  5E  ist  klein.  Die  Annahme  liegt  nahe,  daB 
hier  Re i bgerausche  dominieren,  die  durch  zunehmende  "Glattung" 
der  FiBufer  usw.  schwacher  werden.  Dagegen  ist  AE  bei  10%  Last- 
steigerung  pro  Zyklus  (AE228)  deutlich  positiv  und  zeigt 


Test  Nr. 


AE  224 


AE  225 


AE  226 


AE  228 


Fmax,z  1.01  1.02  1.05  1.10 

F 

max, z-1 


Weitere  Tests  mit  mehr  Zyklen  sind  erforderlich,  um  die  Mechanis- 
men  trennen  und  quantitativ  beurteilen  zu  konnen,  die  bei  derar- 
tigen  Versuchen  zur  akustischen  Energiesumme  beitragen. 


I 


4.3  Faserbiindelbriiche  und  Probenversagen 

Die  Schadigungsmechanismen  in  CFK  sind  mit  verschiedenen  Berei- 
chen  der  SE-Amplitude  verkniipft  [17,18,19,20].  Der  gemessene 
Amp  litudenab  stand  zwischen  Faserbiindelbruchen  und  matrixkon- 
trollierten  Mechanismen  ist  so  groB  (bis  zu  50  dB),  daB  auch  bei 
Nichtberiicksichtigung  der  Dampfung  eine  Unterscheidung  moglich 
ist.  (Da  die  Schallquellen  i.allg.  raumlich  verteilt  in  der  Probe 
liegen  und  daher  vielfaltige  Schal liibertragungswege  durch  das 
stark  anisotrope  Laminat  berucksichtigt  werden  muBten,  ist  eine 
analytische  Beriicksichtigung  der  Dampfung  fiir  die  verschiedenen 
Wellenmoden  beim  technischen  Einsatz  der  SEA  ausgeschlossen ) . 

Bei  unidirektionalen  of f-axi s-Proben  [  0 ] q ,  deren  Versagen  rein 
matrixkontrol liert  und  identisch  mit  'first  ply  failure'  ist, 
trat  kein  Auslaufer  mit  Werten  >  70dB  in  der  Ampli tudenhauf ig- 
kei tsvertei lung  auf  . 

Der  Bruch  von  unidirektionalen  Probestaben  bzw.  multidirektiona- 
len  Laminaten  mit  tragenden  [ 0° } -Lagen  vol lzieht  sich  dagegen  als 
DurchreiBen  aller  Fasern  in  sehr  kurzer  Zeit.  Die  gemessenen 
Scha 1 1 emi ssionen  deuten  auf  DurchriBzeiten  von  unter  50ms  bis 
he;  ah  zu  <  1ms  hin  (zur  Abschatzung  siehe  unten)  . 

Fei  e  l  n>'  i  Probenbreite  von  25mm  ist  die  geschatzte  DurchriBge- 
hw.  ndiakeit  damit  noch  immer  sehr  klein  gegeniiber  der  (longitu- 
•1 1  na  )  f -ha  1 1 geschwindi gkei t  c^  im  Harz  und  erst  recht  in  den 
-  •<  *• !  n  !  •*!  I  :>  hiiB  der  F  robe  als  ganzes  erfolgt  viel  langsamer 
•  !*-•  “Uirelimn  Fasern.  Fur  diese  kann  der  Ansatz  von  Phil- 

:■!!  fin  isrtrrpes  elastisches  Medium  (Rechnung  mit 
•  c  >  Al-rhaf  rung  der  DurchriBzeit  benutzt  werden 
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Bild  14.  Einsetzen  des  Faserbiindelbruches  (  [  0° )  90°4  }  g-Probe ) 

Die  stark  unterschiedliche  DurchreiBgeschwindigkei t  in  den  Fa- 
sern  und  der  Gesamtprobe  bedingt  einen  komplizierten  und  diskon- 
tinuierlichen  Verlauf  des  kollektiven  Faserbruchs  beim  Proben- 
versagen.  Die  Einleitung  dieses  Vorgangs  erfolgt  sehr  schnell 
innerhalb  von  etwa  5-20  ps,  wie  die  Signalanstiegszei t  der  ent- 
sprechenden  Schallereignisse  belegt  (vgl.  auch  Bild  14). 


Diese  GroBe  (RISE  TIME)  kann  als  Ober grenze  der  Anstiegszeit  der 
mechanischen  Schwingungsampli tude  am  Schadensort  akzeptiert  wer- 
den,  da  die  Dispersion  stets  zu  einer  Verbrei terung  des  Signals 
in  der  Zeitdomane  fiihrt.  Analog  gilt  fur  die  DurchriBzeit  tp  der 
Probe 

Tp  <  DURATION 

Direkte  Beobachtung  des  Signalverlauf s  iiber  der  Zeit  fiihrte  zu 
den  o.a.  Werten  im  msec-Bereich,  die  den  Zeitraum  abgrenzen,  in 
denen  das  Signal  noch  nicht  abklingt.  Die  Problematik  des  Uber- 
tragungsverhaltens  und  der  Reflexionen  in  der  Probe  wird  im  Prin- 
zip  beherrscht  [23],  die  rechnerische  Anwendung  bei  Schallemis- 
sionssignalen  in  Faserverbunden  aber  noch  nicht. 

Beim  Standardlaminat  [ O2 I +45 | O2 I -45 | 0 | 90 ] s  sowie  bei  den  [0°]- 
und  den  [  0°  |  90° ] -Proben  wurden  vereinzelt  Schallereignisse  regi- 
striert,  die  von  Faserbriichen  vor  dem  Totalversagen  herriihren 
konnen,  z.B.  EVENT  2179  im  Versuch  AE211  (Daten  siehe  Anhang)  mit 
hoher  Amplitude  im  Vergleich  zum  Energieinhalt  und  deutlich  vom 
Normalfall  abweichender  Frequenzvertei lung . 

Zwei  Eichversuche  an  [  ±45° ] ^ s~Proben  wurden  mit  maximaler  Auflo- 
sung  gefahren  (AE  195-4,-5,  Proben  270/10,  276/14).  Schallereig¬ 
nisse  mit  fur  Faserbruch  typischer  Form  wurden  nicht  gefunden. 


5.  ZUSAMMENFASSUNG  DER  ERGEBNISSE 


Mit  Hilfe  der  Schallemissionsanalyse  lassen  sich  die  drei  Typen 
von  Schadigungsmechanismen,  die  in  CFK  -  Laminaten  auftreten 
konnen,  detektieren: 

•  Kohasionsbruch  der  Fasern 

•  Kohasionsbruch  der  Matrix  (Risse,  Delaminationen) 

•  Adhasionsbruch  Faser-Matrix  (Debonding) 

Der  zeitliche  Ablauf  dieser  Schadigungsmechanismen  kann  verfolgt 

werden.  Hierzu  wurde  eine  Reihe  von  Versuchen  an  CFK-Probestaben 
durchgefiihrt,  die  das  Versuchsprogramm  von  Awerbuch  weiterfiihren 
und  dabei  das  hohe  Auf losungsvermogen  der  PAC  3000/3004  -  Anlage 
ausnutzen  sollten: 

•  Der  Punkt  des  ersten  Auftretens  von  Schallemissionen  iiber 

ansteigender  Zugbelastung  wurde  bestimmt.  Je  nach  Laminat 
liegt  dieser  Punkt  zwischen  15%  (  [  0®  | 90° ) 2S )  und  >90%  ( [ 20° ] g ) 
der  Bruchlast.  Er  zeigt  den  Beginn  matrixkontrollierter  De- 
fekte  an.  Gekerbte  [ O2 I +45 | O2 I -45 | 0 | 90 ) -  Proben  zeigen  Emi s- 
sionen  ab  Belastungsbeginn.  Bei  [0B|90®)_  —  und  [±45“]- 

Proben  weichen  die  Ergebnisse  von  den  von  Awerbuch  erhaltenen 
betrachtlich  ab. 

•  Bei  gekerbten  Proben  lassen  sich  die  ersten  Emissionen  erwar- 

tungsgemaB  an  der  Kerbe  lokalisieren.  Bei  ungekerbten  Proben 
existiert  kein  bevorzugter  Ort  fur  das  Entstehen  erster 
Risse.  Auch  nachdem  die  ersten  mit  SEA  nachgewiesenen  Schaden 
aufgetreten  sind,  ist  hier  keine  Konzentration  der  nachfol- 
genden  Schadensereigni sse  auf  den  Ort  der  ersten 

f estzustel len . 

•  Die  SE-Energiesummenkurven  haben  zunachst  durchweg  einen 
angenahert  exponentiel len  Verlauf.  Mit  steigender  Belastung 
wird  die  Streuung  zwischen  ihnen  auch  bei  gleichartigen  Pro¬ 
ben  unter  identischen  Versuchsbedingungen  sehr  groB,  was  auf 
die  Einmaligkeit  und  Irreversibilitat  jeder  Probengeschichte 
hinweist:  Wahrend  sich  die  individuel len  Schadigungsverlauf e 


anhand  e.xakter  MeBdaten  verfolgen  lassen,  konnen  fur  eine 
Klasse  qleichartiger  Froben  nur  Beschreibungen  gegeben  wer- 
den.  Dies  ist  unbef riedigend  und  erfordert  weitere 
Lntwi cklungsarbei t . 

Faserbruche  wurden  bei  alien  untersuchten  Laminaten  auBer 
den  off-axis-  und  den  [ ±45 0 j -Proben  f estgestel 1 t .  Jedoch  tre- 
ten  nur  vereinzelt  Schal lereigni sse  mit  der  fur  Faserbruch 
typischen  Form  vcr  dem  Totalversagen  auf,  sonst  nur  wahrend 
des  kat astrophalen  Bruches. 

Das  Versagen  der  Proben  mit  tragenden  [0°]-Lagen  vollzieht 
sich  di skontinuierl ich  als  DurchreiBen  von  deren  Fasern  in 
sehr  kurzer  Zeit.  Die  entsprechenden  Schallemissionen  deuten 
auf  DurchriBzei ten  von  unter  50  ms  bis  herab  zu  S  1  ms  bin. 
Der  Ubergang  zwischen  dem  vorhergehenden  Schadigungsver lauf 
und  dem  Einsetzen  des  DurchreiBens  ist  auBerordentlich 
schar f  (5-20  us)  . 

Bei  langsamer  zyklischer  Belastung  der  Proben  konnen  auf- 
grund  des  Feiicity-Verhaltni sses  und  der  in  den  einzelnen 
Belastungszyklen  f  reigesetzten  akustischen  Energie  Riick- 
schliisse  auf  das  MaB  der  Schadigungen  gezogen  werden.  Weitere 
Untersuchungen  zu  diesem  Punkt  sind  erf orderlich. 
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ANHANG 


Beispiel  eines  vollstandigen  Datensatzes 

Im  folgenden  sind  die  Schallemissionsdaten  eines  Zugversuchs  an 
einem  [  0^  |  +45 1 02  | -45 1 0 1  90  ]  g  -  Probestab  mit  GFK-Schlauchmantel 
aufgelistet,  und  zwar  die  ersten  100  sowie  die  letzten  213  der 
insgesamt  2213  registrierten  Schal lereigni sse . 

•  Die  Schallereignisse  (EVENTS)  sind  fortlaufend  numeriert. 

•  TIMING  bezeichnet  die  Lauf zei tdif f erenz  der  Signale  zu  den 
beiden  Transducern  in  ysec.  *HIT*  ist  die  Signatur  des  zuerst 
bzw.  einzig  angesprochenen  Kanals. 

•  CH  ist  die  Kanalnummer  (es  wurde  zweikanalig  gemessen). 

•  TIME  ist  die  Echtzeit.  Auflosung  hier  0.01  sec 

•  DURATION  ist  die  Ereigni sdauer  in  ysec. 

•  COUNTS  sind  die  Uberschreitungen  des  Schwellwerts 
(Threshold)  durch  die  Signalspannung.  Angegeben  ist  ihre 
Zahl  pro  Kanal  und  Event  minus  1. 

•  ENERGY  ist  das  relative  EnergiemaB  der  Schallereignisse. 

(-*  siehe  nachste  Seite) 

•  AMPLITUDE  ist  die  Spi tzenampli tude  der  Schallereignisse.  Der 
Nullpunkt  der  Skala  ist  definiert  als: 
Transducerausgangsspannung  1  jjV  =  0  dB 

•  RISE  TIME  ist  die  Anstiegszeit  der  Signalspannung.  (Die  hohen 

10 

Werte  >1024=2  ,  die  mit  Duration=0  einhergehen,  beruhen  auf 

einem  meBtechni schen  Effekt  bei  hohen  Datenraten.  Er  ist  fur 
die  Schadigungsbeurtei lung  irrelevant). 


Zur  Abschatzung  der  Energie  der  Schallereignisse 


Die  Energy-Counts  sind  ein  relatives  MaB  fur  die  in  die  Sensoren 
dissipierte  Energie.  Deren  absoluter  Betrag  kann  unter  vereinfa- 
chenden  Annahmen  wie  folgt  abgeschatzt  werden: 

Die  maximale  Signalspannung  am  Sensorausgang  erhalt  man  aus  der 
GroBe  AMPLITUDE: 

0  dB  20  dB  40  dB  60  dB  80  dB  100  dB 
1  yV  10  yV  100  yV  1  mV  10  mV  100  mV 

Die  Impedanz  des  Eingangskrei ses  (typ.  400  pF  Kapazitat  des  Sen¬ 
sors,  lOkfi  ||  15pF  Eingangsimpedanz  des  Vorverstarkers)  liegt  im 
150kHz-Bereich  bei  ca.  10k$2.  Naherungswei  se  ergibt  also 
U2  [  Volt2 1/101*  [  fi  ]  die  max.  Momentanlei stung.  Sie  liegt  bei  40dB- 
Events  in  der  GroBenordnung  Picowatt. 

Mit  der  Ereigni sdauer  D  laBt  sich  so  die  Energie  in  Pico-  bzw. 
Nanojoule  abschatzen.  Genauere  Ergebnisse  erfordern  die  Beruck- 
sichtigung  der  Frequenzverteilung  (d.h.  die  Berechnung  der 
"spektralen"  Energiedichte ) . 
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HEADER  INFORMATION 


FILE  B:AE211.DTA  140  KBYTES  2213  EVENTS  IN  #:  0:31.4 
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44.05  47. 
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Event  2199  stellt  das  Probenversaqen  dar.  Anders  als  bei 
unidirektionalen  Proben  fallt  die  Last  nicht  unmittelbar 
auf  Null  ab,  sondern  wird  wahrend  des  AuseinanderreiBens 
iiber  einen  Zeitraum  von  0.08  sec  abgebaut.  Die  Dehnunq 
geht  unmittelbar  zur  Zeit  0:0:36.66  gegen  unendlich  (rea. 
Soannung  10  Volt  =  overflow). 

Der  Schallaufweg  zu  den  Transducern  1  und  2  war  extrem 
ung.leich  lang  (Lauf zeitdif f erenz  808.4  usee).  Dieser 
Unterschied  bedingt  die  starke  DSmofunq  aut  Kanal  2. 

Der  Event  ist  dort  nur  noch  halb  so  lang  wie  auf  Kanal  1 
und  die  Amplitude  ist  um  7.4  dB  (rund  Faktor  2)  vermin- 
dert.  Die  in  die  Transducer  dissipierten  Eneroien  ver- 
halten  sich  folgerichtig  wie  4:1. 

In  der  Mehrzahl  der  Falle  ist  dieses  Ungleichqewicht  weit 
schwacher  ausqepragt  (Versagen  in  Probestabmitte ,  oeringere 
Damofunn  bei  anderen  Laminaten) . 

Die  Damofung  im  Laminat  mul3  stets  beriicksichtict  werden, 
wenn  quantitativ  auf  mechanische  Schadigungsenergien  riick- 
aeschlossen  werden  soli. 


‘v'V  -  ■- J 

►  V  -■  -»  - 


r  » 

.  •*. 


'  V  ■/ 


V  >.■ 


•  •.  V’  •_ 

.  .*  .-  .» 
.*•  -V.-. 

v  v  v 


APPENDIX  B 


On-Line  Measurement  of  Onset  and  Growth  of 
Edge- Delaminations  in  CFRP- Laminates  by  an 
Optical  Grating  Reflection  Method 

by 

R.  Schiitze  and  H.  C.  Goetting 
DFVLR  -  I  ns  ti  tut  filr  Strukturmechanik 


DEUTSCHE  FORSCHUNGS-  UND  VERSUCHSANSTALT 
FUR  LUFT-  UND  RAUMFAHRT  E.V. 
INSTITUT  FUR  STRUKTURMECHANIK 


Interner  Bericht 
IB  131-84/49 


On-Line  Measurement  of  Onset  and  Growth  of 
Edge-Delaminations  in  CFRP-Laminates  by  an 
Optical  Grating  Reflection  Method 


Braunschweig,  Oktober  1984 

Dieser  Bericht  urnfaBt: 

20  Seiten  mit 
9  Bildern 


Die  Bearbeiter: 

(Dipl.-Ing.  R.  Schutze) 

1/  ■  /, 

(Dipl.-Phys.  H.C.Goetting) 


Contents 


1.  Introduction 

2.  Structure  mechanical  aspects  of  edge-delamination  development 

3.  Description  of  the  optical  grating  reflection  method 

4.  Edge-delamination  growth  under  quasi-static  tensile  load 

5.  Evaluation  of  surface  deformations 

6.  Conclusions 

7.  References 

bstract 

.  jie  onset  and  growth  of  edge-delaminations  in  multidirectional 
'FRP-specimens  under  quasi-static  and  fatigue  loading  was  deter¬ 
mined  by  the  application  of  a  new  optical  grating  reflection  me¬ 
thod.  This  procedure  allows  the  on-line  observation  of  surface 
deformations  corresponding  to  the  actual  size  of  propagating  ed¬ 
ge-delaminations  . 

The  evaluation  of  the  grid  patterns  described  the  growth  of  the 
the  delamination  versus  increasing  load  or  number  of  cycles,  re¬ 
spectively,  and  permitted  the  measurement  of  out-of-plane  defor¬ 
mations  with  relatively  high  resolution. 

Optisches  Raster- Reflexionsverfahren  zur  on-line  Messung  des  Randdelamina- 
tionsfortschritts  an  CFK-Laminaten 

Ubersicht 

Im  quasi stati schen  und  Schwingf estigkeitsversuch  wurden  die  Ent- 
stehung  und  der  Fortschritt  von  Randde 1 aminationen  in  multidi- 
rektionalen  CFK-Laminaten  mit  einem  optischen  Raster-Ref lexions- 
verfahren  bestimmt.  Das  Verfahren  gcstattet  wahrend  des  Versuchs 
die  direkte  Beobachtung  der  Oberf lachenverf ormungen,  die  durch 
die  Ausdehnung  der  sich  ausbildenden  Randdelaminationen  bedingt 
sind . 

Die  Auswertung  der  Rasterbilder  ermoglicht  die  Bestimmung  des 
Delaminationsfortschri tts  in  Abhangigkeit  von  der  Last  bzw. 
Lastspielzahl ,  sowie  die  Messung  von  out-of-plane  Verformungen 
mit  relativ  hoher  Auflosung. 
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Composite  materials  like  carbon-fibre-reinforced  plastics  (CFRP) 
are  laminated  from  unidirectional  layers  with  preferred  fibre 
directions  of  0°,  ±45°  and  90°.  Certain  loading  conditions  can 
lead  to  local  separations  between  adjacent  layers  characterized 
as  delaminations.  Under  tensile  loading  of  laminated  test  speci¬ 
mens  the  development  of  such  delaminations  at  the  free  edges  is 
frequently  observed. 

The  observation  and  measurement  of  the  delamination  growth  under 
increasing  load  or  number  of  cycles,  respectively,  is  usually 
performed  by  means  of  ultrasonic  C-scans  or  X-ray  radiography 
[1],  [2],  and  in  some  cases  crack-opening  displacement  (COD)  mea¬ 
surement  methods  are  applied  (3]  .  The  former  two  methods  have  the 
disadvantage  of  permitting  off-line  evaluations  only,  requiring 
unnecessary  interruptions  of  loading  and  testing,  while  COD-mea- 
surements  can  only  refer  to  specific  points  of  the  delamination 
area . 

Therefore  an  optical  grating  reflection  method  was  introduced 
suitable  for  the  on-line  determination  of  the  entire  delaminated 
zones  [4).  This  simple  procedure  allows  the  measurement  of  the 
delamination  size  during  testing  and  the  evaluation  of  curvatures 
and  out-of-plane  deformations. 

Structure  Mechanical  Aspects  of  Edge-Delamination  Development 

Because  of  the  different  stiffnesses  of  unidirectional  layers 
oriented  under  0°,  ±45°  and  90°  relative  to  the  loading 
direction,  different  stresses  develop  in  the  individual  layers  of 
multidirectional  laminates. 

Tensile  loading  of  laminates  with  a  stacking  seguence  of 
[  02°/+45o/02O/-45°/0o/90o  )  g  generates  at  the  specimen  edges  high 
tensile  stresses  0^3  in  the  direction  of  the  laminate  thickness  - 
so  called  peeling  stresses  -  which  cause  the  separation  of  dela¬ 
minated  parts  of  a  specimen,  see  Fig.  1  [5).  These  high  stresses 
in  the  90°  layers  at  the  0°/90°-interf aces  stem  from  the  trans¬ 
verse  contraction  constraints  of  the  90°  layer  on  the  neighboring 
0°  and  45°  layers. 
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Corresponding  to  the  results  of  stress  analyses  performed  by  Roh- 
wer  [5]  edge  delaminations  develop  at  the  interface  between  0° 
and  90°  layers  with  intermittent  "jumps"  to  the  adjacent  interfa¬ 
ce.  This  typical  transition  from  one  interface  to  the  other  is 
aided  by  the  existence  of  a  matrix  crack  pattern,  called  the  "cha¬ 
racteristic  damage  state"  (CDS),  which  is  described  in  detail  by 
Reifsnider  et  al .  [6J. 

Depending  on  the  load  and  the  actual  size  of  the  edge-delamina¬ 
tion,  a  certain  amount  of  bending  of  the  delaminated  sections 
occurs.  The  correspor ding  curvature  changes  can  be  evaluated  by 
the  optical  grating  reflection  method. 

Description  of  the  Optical  Grating1  Reflection  Method 

Fig.  2  shows  the  principle  of  this  measurement  method  as  proposed 
by  Ritter  et  al.  [7].  An  edge-delaminated  CFRP- specimen  serves  as 
"object"  the  cross-section  of  which  is  displayed  by  dashed  lines 
for  the  delaminated  state  and  by  straight  lines  for  the  un-dela- 
minated  state.  A  special  specimen  surface  treatment  of  the  speci¬ 
men  reflecting  like  a  mirror  is  necessary.  This  mirror  surface 
allows  the  observation  of  the  virtual  image  of  an  illuminated 
line-grating. 

The  lines  of  this  grating  are  parallel  and  perpendicular  to  the 
picture  plane,  that  means  the  image  lines  are  parallel  to  the  lon¬ 
gitudinal  direction  of  the  test  specimen.  During  loading  the 
cross-  section  changes  its  shape  into  the  dashed  line  figure  and 
therefore  point  P  moves  into  the  position  P'  .  Because  of  the  chan¬ 
ges  of  the  surface  curvature  the  observer,  looking  from  the 
camera  position  C,  identifies  line  5  instead  of  line  3  .  These 
distortions  of  the  grating  lines  can  be  directly  observed  or  re¬ 
gistered  by  a  camera.  By  means  of  the  geometrical  relations  it  is 
possible  to  evaluate  these  distortion  patterns  quantitatively. 

Fig.  3  shows  the  arrangement  of  a  grating  reflection  measurement 
in  a  servohydraulic  testing  facility.  The  main  part  is  a  light¬ 
box  with  two  gratings  attached  symmetrically  to  the  transverse 
axis  of  the  specimen,  and  adjustable  by  small  micro-meter  screw 
drives.  The  camera  support  is  also  adjustable  in  order  to  facili- 


tate  the  exact  positioning  of  the  light-box  as  well  as  the  camera 
objective  for  quantitative  curvature  evaluations.  Two  gratings 
are  used  to  simultaneously  observe  the  front  and  back  side  of  the 
specimen. 

The  mirror-like  preparation  of  the  test  specimen  must  be  applied 
such  that  it  does  not  influence  the  imposed  deformations  of  the 
specimen . 

An  evaluation  of  different  preparation  procedures  showed  that 
mirror-like  specimen  surfaces  are  best  obtained  by  attaching  a 
thin  epoxy-resin  film  under  pressure  against  a  plate-glass  of 
perfect  smoothness.  This  thin  film  also  smoothes  the  roughness  of 
the  specimens  surface,  and  cannot  significantly  affect  the  speci¬ 
men  deformation  because  of  its  low  Young's  modulus.  The  optical 
contrast  in  the  image  between  the  dark  grating  lines  and  the  re¬ 
flected  light  provided  sufficient  resolution. 

This  simple  technique  permits  the  preparation  of  extensive  series 
of  test  specimens  with  an  effective  preparation-  test  time  ratio. 

Edge-Delamination  Growth  under  Quasi-Static  Tensile  Load 

Pre-test  results  showed  the  onset  of  edge-delaminations  during 
quasi-static  loading  to  occur  only  at  relative  high  load  levels 
near  the  ultimate  load.  In  order  to  avoid  the  difficulty  of  detec¬ 
ting  and  measuring  the  delamination  growth  in  such  a  narrow  range 
the  specimens  were  first  cycled  under  tension-tension  fatigue 
loading  with  R=0 . 1  and  at  a  frequency  of  5  Hz  until  onset  of  dela¬ 
minations  could  be  observed. 

All  of  the  tests  were  conducted  in  a  load-controlled  servo-hy¬ 
draulic  testing  machine  at  constant  ambient  temperature  of  22°C 
and  at  20%  R.H.  Fatigue  loads  as  well  as  static  loads  were  applied 
by  means  of  a  preprogrammed  microprocessor-controlled  function 
generator  aiding  fatigue  interruptions. 

The  specimens  were  first  cycled  at  a  maximum  load  of  500  N/mm2 
followed  by  cycling  at  600  N/mm2,  with  2000  cycles  at  each  level. 
An  increase  of  the  maximum  load  to  700  N/mm2,  equivalent  to  ca. 
70%  of  the  ultimate  load,  induced  the  onset  of  delaminations  af¬ 
ter  120-160  cycles.  After  this  initial  damage  state  step-by-step 
quasi-static  loading  was  applied  using  a  loading  rate  of  5  N/s.  At 


each  level  the  load  was  retained  for  a  period  of  300  s  in  order  to 
measure  time-dependent  phenomena  of  delamination  growth  [8]. 

Fig.  4  shows  a  series  of  observed  grating  patterns  of  a  typical 
test  specimen.  The  first  picture  of  this  series  displays  the  in¬ 
itial  damage  state  after  cycling  at  700  N/mm2  indicating  the  onset 
of  a  delamination  at  the  right-hand  edge  by  the  distortion  of  the 
grating  lines.  The  next  picture  represents  the  rapid  increase  of 
the  delaminated  area  along  the  entire  test  length  due  to  a  small 
load  increment  of  12.5  N/mm2;  the  left-hand  edge  is  now  partially 
delaminated. 

Following  the  progression  of  the  delamination  over  the  complete 
left  side  of  the  specimen,  the  remaining  pictures  show  the  de¬ 
crease  of  the  delamination  growth  rate  which  can  be  observed  by 
regarding  the  undistorted  lines  in  the  center  of  the  specimens 
surface . 

In  Fig.  5  the  growth  of  left-hand  and  right-hand  edge-delamina¬ 
tions  of  three  specimens  versus  load  is  depicted.  The  lower  va¬ 
lues  represent,  at  each  load  level,  the  delamination  size  after 
load  increase,  whereas  the  higher  values  reflect  the  influence  of 
time  under  constant  load.  These  curves  demonstrate  a  rapid  growth 
of  the  delamination  size  immediately  after  onset  due  to  load  in¬ 
crease  as  well  as  to  holding  time.  They  also  depict  the  typical 
decrease  of  the  growth  rate  with  higher  loads. 

Because  of  the  acute  angle  between  specimen  surface  and  camera 
axis,  see  Fig.  1,  the  grating  picture  always  displays  simulta¬ 
neously  both  an  edge  and  a  side.  Therefore  the  edge-opening  at 
higher  loads  and  corresponding  higher  deformations  of  the  delami¬ 
nated  layers,  can  be  observed  as  demonstrated  in  Fig. 6  for  the 
right-hand  edge.  Especially  the  "jumping"  of  delamination  bet¬ 
ween  the  0°/90°-interf aces  can  be  clearly  recognized  in  the 
wave-like  shape  of  the  distorted  grating  lines.  This  behaviour  is 
explained  by  the  different  thicknesses  of  the  delaminated  parts 
leading  to  distinct  bending  stiffnesses  affecting  this  wave-like 
crack- opening . 

In  order  to  quantify  the  observed  delaminated  areas 
independently.  X-ray  radiography  was  applied  after  the  damage 
growth  tests.  The  actual  damage,  i.e.  delaminations  as  well  as 


matrix  cracks,  can  be  made  visible  with  t 
dye-penetrant.  Fig.  7  shows  a  comparison  of  the  lai 
ture  of  a  typical  test  specimen  and  X-ray  radi 
tablishing  good  correlation  between  the  areas  of 
ng  lines  and  the  dark  areas  of  the  X-ray  radiogr 
no  the  delaminations. 


A  more  detailed  description  of 
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>ile  load  of  825  N/mm* .  This  relati’ 


Conclusions 


After  an  overwiew  of  the  structure  mechanical  aspects  of  edge- 
delaminations  in  multidirectional  CFRP-laminates,  a  simple  opti¬ 
cal  grating  reflection  method  is  described  which  offers  the  de¬ 
tection  of  onset  and  growth  of  delaminations.  The  evaluation  of 
the  distortion  of  grating  lines  observed  on  the  mirror-like  spe¬ 
cimen  surface  allows  the  identification  of  delaminated  areas 
during  the  test. 

The  interpretation  of  grating  reflection  pictures  after  step-by- 
step  quasi-static  loading  led  to  typical  curves  for  the  onset  and 
propagation  of  edge-delaminations  versus  load,  including  time- 
dependency.  X-ray  radiographic  records  showed  good  agreement 
with  the  optically  observed  delamination  size. 

The  calculation  of  curvatures  obtained  by  the  evaluation  of  gra¬ 
ting  pictures  after  fatigue  loading  yielded  high-resolution  data 
of  out-of-plane  deformations  in  the  transverse  direction. 

This  method  described  seems  to  be  a  valuable  tool  for  the  observa¬ 
tion  and  evaluation  of  delaminations  and  surface  deformations  of 
CFRP-laminates  because  of  its  simple  preparation  technique,  easy 
handling,  and  relatively  high  resolution. 
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due  to  axial  tensile  load,  calculated  by  Rohwer  [5]  . 


Specimen  width  ,  mm  - - 


Fig.  9a:  Deformation  in  transverse  direction  due  to  left-hand 
edge-delamination . 


Specimen  width  .  m m 


Fig.  9b:  Deformation  in  transverse  direct  •>  m  due  to  right-hand 
edge-delaminat ion . 
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Bezeichnungen 

Es  gilt  die  Summenkonvention,  bei  der  liber  zwei  gleiche,  nicht 
eingeklammerte  Indizes  summiert  wird.  Wenn  es  nicht  explizit 
anders  geregelt  ist,  laufen  lateinische  Indizes  von  i  =  1,2,3 
und  griechische  Indizes  von  a  =  1,2. 


Skalare  GroBen: 


(  ) 

Anf angswert 

(•) 

Inkrement 

(  )T 

Transponi erte 

<  >'i 

partielle  Ableitung  nach  der 

Koordinate  x1 

1  O 

X  ,  X 

mm 

Kartesische  Koordinaten 

r ,  0 

mm,  1 

Polarkoordinaten 

dV 

mm 3 

Volumeninkrement 

dF 

mm 2 

Flacheninkrement 

F1,  F1 

N 

Einzelkraf te 

Fu 

mm2 

Oberflache  mit  Wegrandbedingungen 

F 

0 

mm2 

Oberflache  mit  Kraf trandbedingungen 

F 

c 

mm2 

RiBf lache 

6  F 

c 

mm2 

Flache  des  RiBf ortschri tts 

1 

c 

mm 

RiBlange 

c 

mm 

Lange  eines  fiktiven  Risses 

Ac 

mm 

Teillange  von  c 

G 

N/mm 

Gesamte  Energief rei setzungsrate 

C,(K,2) 

N/mm 

Energief rei setzungsrate  fiir  die 

RiBof fnung 

G2(K22) 

N/mm 

Energief  rei  setzungsrate  fiir  die 

RiBglei tung 

K 

a 

N/mm3/2 

Spannungsintensi tatsf aktoren 

KC 

a 

N/mm3/2 

kri ti sche  Spannungsintensi tatsf aktoren 

R 

o 

mm 

Bezugsradius 

S' 

H 

Nmm 

Mater i a lparameter  fiir  die  Verfestigung 

Y 

1 

FlieBbedingung 

J<  JO. 


p 

e 

i 

Plastische  Vergleichsdehnung 

°E 

N/mm2 

Vergleichsspannung  an  der  Elastizi tats 

°U 

N/nuu2 

grenze 

Grenzf estigkei t  des  Werkstoffs 

P  =  r/Ro 

1 

bezogener  Radius 

* 

N 

Spannungsfunktion 

w<sij> 

N/mm2 

Verzerrungsenergie  je  Volumeneinheit 

W(  o ) 

N/mm2 

konjugierte  Verzerrungsenergie  je 

Tensorielle 

GroBen: 

Volumeneinheit 

1 

J 

u . 

i 

Verschiebungen 

u . 

i 

Randverschiebungen 

Au . 
l 

RiBof f nungsverschiebung 

-i 

Pv 

Volumenbelastung 

-i 

Pp 

Oberf lachenbelastung 

r . 

j. 

Komponenten  des  Einhei tsvektors  normal 

zur  Oberflache 

Verzerrungen 

o ij 

Spannungen 

V  .  . 

Gradient  der  FlieBflache 

IJ 

a  .  . 

Koef f izienten  des  Stof f gesetzes 

ij 

Vektoren,  Matrizen: 

(Nachgiebigkei t ) 

z 

Losungsvektor 

z  . 

i 

Losungsvektor  der  interne  Variablen 

z 

0 

Losungsvektor  der  externe  Variablen 

A.  . 

Untermatri zen  der  Elementmatrix 

hZ2'  Ke 

reduzierte  Elementmatrix 

P2'  Pe 

reduzierter  Lastvektor 

P 

Lastvektor 

L 

obere  Dreiecksmatrix  von  A 

y 

vorreduzierter  Losungsvektor 

1 .  Auf gabenstellung 

Laminate  aus  kohlef aserverstarkten  Kunststoffen  (CFK-Laminate) 
bestehen  aus  dtinnen  Einzelschichten  mit  unidirektionalen  Faser- 
lagen.  Die  Festigkeiten  und  die  Steif igkeiten  der  unidirektiona¬ 
len  Schichten  (UD-Schichten)  sind  in  Faserrichtung  mehr  als  eine 
Zehnerpotenz  groBer  als  quer  zur  Faser.  Optimiert  man  die  Faser- 
winkel  in  den  Einzelschichten,  kann  -  bei  minimalem  Gewicht  - 
die  Festigkeit  der  Laminate  weitgehend  an  vorgegebene  Beanspru- 
chungen  angepaBt  werden.  Dadurch  entstehen  sehr  leichte  und 
feste  Strukturen,  wie  sie  in  der  Luft-  und  Raumfahrt  und  neuer- 
dings  auch  im  Fahrzeugbau  verwendet  werden. 


Aufgrund  der  geringen  Festigkeit  der  Kunststof fmatrix  reiGen  die 
Laminate  bevorzugt  in  Faserrichtung  bzw.  zwischen  den  Schichten 
auf.  Diese,  am  Anfang  meist  kleinen  Schaden,  konnen  sich  mit 
fortschrei tender  Belastung  oder  mit  wachsenden  Lastzyklen  aus- 
weiten  und  den  globalen  Bruch  einleiten.  Versuche  an  Laminaten 
zeigen  eine  Vielfalt  von  unterschiedlichen  Briichen  und  Bruchkom- 
binationen.  Zur  Zeit  kann  ein  lokaler  Schaden  und  der  sich 
daraus  entwickelnde  Bruch  des  Laminates  rechnerisch  noch  nicht 
erfaBt  werden.  Urn  diesem  Ziel  naher  zu  kommen,  werden  Versuche 
an  UD-Schichten  durchgefiihrt  und  parallel  dazu  tiber  Vergleichs- 
rechnungen  Bruchkri terien  ermittelt. 

Derartige  Berechnungen  sind  sehr  aufwendig,  da  sich  mit  fort- 
schreitendem  RiB  die  Geometrie  der  Struktur  andert.  Urn  hier  zu 
vertretbaren  Rechenzeiten  zu  kommen,  wurde  eine  Losungsmethode 
entwickelt,  bei  der  die  linear  elastische,  ungerissene  Struktur 
durch  Zwangsbedingungen  eingeschrankt  wird,  die  den  Bruch  be- 
schreiben.  Bei  jedem  RiBf ortschri tt  werden  nur  die  Gleichungen 
fiir  den  neu  auftretenden  RiB  gelost  und  damit  die  bekannten  Er- 
gebnisse  fiir  die  vorangegangene  RiBlange  modifiziert.  Durch  eine 
neuartige  Tei lstrukturtechnik  wird  der  iiberwiegende  Teil  der 
Struktur  abgelost,  so  daB  sich  die  eigentliche  Berechnung  des 
RiBf ortschri tts  nur  auf  die  unmittelbare  Umgebung  des  Risses 
(oder  der  geschadigten  Zone)  erstreckt. 


2.  Die  Grundqleichungen  fur  linear  elas tische  Strukturen 

Eines  der  allgemeinsten  Prinzipe  der  Kontinuumsmechanik  ist  das 
Prinzip  der  virtuellen  Verruckung.  Es  beschreibt  nur  die  Gleich 
gewichtsbedingungen,  die  auch  bei  nichtlinearen  Problemen  er- 
fiillt  sein  mussen.  In  diesem  Prinzip  sind  die  Spannungen  und  di 
Belastungen  die  wirklich  auftretenden  GroBen,  wahrend  die  virtu 
ellen  Verriickungen  Variationen  der  WeggroBen  sind,  die  die  Weg- 
randbedingungen  erfiilien  mussen.  Erganzt  man  das  Prinzip  der 
virtuellen  Verruckung  urn  die  Nebenbedingungen  fur  die  Wegrandbe 
dingungen,  die  RiBoffnung  und  den  RiBf ortschri tt,  entsteht  das 
erweiterte  Prinzip: 

6JE=  [la  '*  ■  <Su,-  -  p  '  -Su.ldV 

J  rj  v  \ 

V 

-/ Pp  -  Su.-dF 
F<r 

^  1  j  -  6  J  a'*  •  r .  •  (u  j  -  u  j)  d  F 

-s  J  •  Aui  ,dF 

Fc 

-  J  ^  c'1’  r.  •  5A  u  .  •  dF  =  0 

«Fc 


Gleichgewicht 

Wegrandbedingung 

RiBoffnung ,  siehe  [1] 

RiBfortschriit,  siehe  [21,  S.  207 


Darin  beschreiben  o 1 -1  die  Spannunoen,  u^  die  Verschiebungen  und 
Au.  die  Relativverschi ebungen  der  beiden  RiBufer.  p.1.  und  p^ 

1  V  r 

sind  die  Volumen-  bzw.  Oberf lachenlasten,  u_^  die  vorgegebenen 

Wegrandbedingungen  und  r^  die  Komponenten  des  nach  auBen  zeigen- 

den  Einhei tsvektors  normal  zum  Rand.  F  und  F  sind  die  Oberfla- 

o  u 

chen,  auf  denen  die  Spannungs-  bzw.  Wegrandbedingungen 
vorgegeben  sind  und  Fc  ist  die  Flache  des  Risses.  liber  doppelt 
vorkommende,  nicht  eingeklammerte  lateinische  Indices  ist  von 
i  =  1,2,3  zu  summieren. 

Der  erste  Term  der  Gleichung  (1)  wird  mit  der  Legendre schen 
Transformation 


(2) 


rij  -syij  =  5W(yjj)  =  6(crij  -  u;  j- W(oij )) 


umgeformt,  wobei  W(3\j)  die  Verformungsenergie,  W(o^^.)  die  kon- 
jugierte  Verformungsenergie  und  die  Verzerrungen  kennzeich- 
nen.  Damit  geht  die  Gleichung  (1)  liber  in  ein  erweitertes  Prin- 
zip  von  Hellinger-Reissner  [3],  das  zusatzlich  Terme  fur  die 
RiBoffnung  enthalt: 

6  JR  =  5  <  f  ( or ' j  •  u .  . -  W { a ' j )  - p '  -u.)  dV 

I  J  t.J  v  I 

\  V 

(3)  -/ pF  ■  u  j  -  d  F  -  /  a ,J  •  r  -(Uj-UjIdF 

'  '  f  c 

ra  ru 

- /.r'i-rj  Au.dF  j  V«4u.-dF  =0 

F  /  <  F 


Darin  werden  neben  den  Spannungen  und  Verschiebungen  auch  die 
Flache  6Fc  flir  den  RiBf ortschri tt  unabhangig  variiert.  Das  Funk- 
tional  erflillt  alle  Grundgleichungen  flir  das  Kontinuum.  Bei  ei- 
ner  numerischen  Approximation  der  ZustandsgroBen  miissen  die  An- 
satzfunktionen  lediglich  gewisse  Stetigkeitsbedingungen  befrie- 
digen . 


Mit  dem  Prinzip  von  Hellinger-Reissner  wurde  ein  geschichtetes 
WeggroBenelement  zur  Berechnung  von  Laminaten  entwickelt,  das  an 
den  Schichtgrenzen  auch  die  Obergangsbedingungen  flir  die  Span¬ 
nungen  erflillt.  Mit  diesem  Element  wurde  ein  Teil  der  in  [13] 
angegebenen  Vergleichsrechnungen  durchgeflihrt .  Da  die  folgenden 
Untersuchungen  nicht  an  ein  bestimmtes  Elementmodell  gebunden 
sind,  wird  hier  auf  eine  Herleitung  des  in  [4]  beschriebenen 
Elementes  verzichtet. 


3.  Nichtlineare  Effekte 

Die  hier  betrachteten  nichtlinearen  Einfliisse  werden  je  nach  Art 
der  numerischen  Behandlung  unterschieden  in 


V\  *“ 


yyjvy.y 

W 


Werkstof fnicht linear i tat , 
Strukturnichtlinearitat  und 
geometrische  Nichtlinearitat. 


Ein  Beispiel  fur  die  Werkstof fnichtlineari tat  ist  das  Plastifi- 
zieren  eines  Korpers.  Die  Variablen,  die  das  plastische  FlieBen 
beschreiben,  sind  an  den  Werkstoff  und  damit  an  einen  bestimmten 
Punkt  der  Struktur  gebunden.  Sie  lassen  sich  bei  einer  finiten 
Ubersetzung  einzelnen  Elementen  in  der  plastischen  Zone 
zuordnen. 

Ein  typisches  Beispiel  fur  die  Strukturnichtlinearitat  ist  ein 
fortschrei tender  RiB.  Durch  den  RiBf ortschri tt  andert  sich  das 
statische  System,  wodurch  die  Berechnung  nichtlinear  wird.  Die 
Variablen,  die  den  RiB  beschreiben,  sind  Strukturvariable,  d.h. 
sie  konnen  nur  der  Struktur  als  Ganzes  zugeordnet  werden  aber 
nicht.  einem  bestimmten  Punkt  innerhalb  der  Struktur.  So  ist  z.B. 
der  Spannungsintensi tatsf aktor  immer  an  die  RiBspitze  gebunden 
und  wandert  mit  dieser  mit. 

Bei  der  geometri schen  Nichtlinearitat  wird  der  Berechnung  die 
verformte  Struktur  zugrunde  gelegt.  Dadurch  miissen  bei  einer  fi¬ 
niten  Ubersetzung  die  Elemente  fortlaufend  an  die  sich  andernde 
Geometrie  angepaBt  werden.  Die  Verformungen  beeinflussen  die 
gesamte  Struktur  und  nicht  nur,  wie  in  den  vorangegangenen  Bei- 
spielen,  einige  Tei lbereiche .  Entsprechend  aufwendig  gestaltet 
sich  auch  die  Berechnung. 

In  den  folgenden  Abschnitten  wird  die  Methode  der  Randgleichun- 
gen  erlautert,  mit  der  die  hier  angesprcchenen  nichtlinearen 
Probleme  mit  geringerem  numerischen  Aufwand  gelost  werden 
konnen . 


3 . 1  Die  Werkstof fnichtlineari tat 

Von  den  moglichen  Werkstof fnichtlinearitaten  soli  hier  nur  die 
Plastizitat  naher  betrachtet  werden.  Andere  Einfliisse  wie  z.B. 
Mikrobriiche,  Coulombsche  Reibung  u. a.  lassen  sich  analog  erfas- 
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sen.  Eine  anschauliche  Deutung  der  plastischen  Verzerrungen  ist 
in  [5]  angegeben. 


Grundlage  der  Plastizi tatstheorie  ist  die  Existenz  einer  skala- 
ren  FlieBbedingung,  die  den  Spannungszustand  begrenzt.  Sie  ist 
im  einfachsten  Fall  nur  von  den  Spannungen  und  bei  verfestigen- 
den  Werkstoffen  noch  von  den  einachsigen  plastischen  Vergleichs- 
dehnungen  abhangig.  Eine  einfache  FlieBbedingung  hat  die  Form 

(4)  V  =  Y  (<rij.  £P  ;  crE.  ou)<0  . 

Darin  sind  die  Spannung  an  der  Elastizitatsgrenze  und  die 
Grenzfestigkei t  des  Werkstoffes,  beide  gemessen  im  einachsigen 

p 

Zugversuch.  Die  plastische  Dehnung  z  beschreibt  die  Verfesti- 
gung  des  Materials  und  damit  die  Aufweitung  der  FlieBflache  un- 
ter  wachsender  Beanspruchung,  Bi Id  1 . 

Das  totale  Differential  der  FlieBbedingung  lautet 

(5)  Y  =  y..  •  b"  -H  •  €P  =  0 
mi  t 


(6) 

H  =--^-1 

9eP  Jep--£p 

Darin  sind  (  ^  )  die  bekannten  ZustandsgroBen  der  vorangegangenen 
Laststufe  und  (  )  die  einer  Laststeigerung  zugeordneten  Inkre- 
mente  (Differentiale) .  Fur  einen  Werkstoff  ohne  Verfestigung  ist 
H  =  0,  so  daB  die  Spannungsinkremente  in  der  Tangentialebene  der 

A 

FlieBflache  liegen.  Der  Verf estigungsparameter  H  beschreibt  die 
Aufweitung  der  FlieBflache.  Dadurch  dreht  sich  der  Spannungsten- 
sor  o1J  nach  auBen  und  beschreibt  einen  Punkt  auf  der  neuen 


10 


FlieBflache . 


Wahrend  die  FI ieBbedingung  (4)  erne  nichtlineare  Funktion  der 
Spannungen  und  der  plastischen  Verzerrungen  ist,  beschreibt  das 
Differential  (5)  einen  linearen  Zusammenhang  zwischen  den  Inkre- 
menten.  Daher  wird  auch  das  Funktional  (3)  fur  die  ungerissene 
elastische  Struktur  fur  die  Inkremente  anaeschr ieben  und  nach 
der  Methode  von  Lagrange  mit  der  FlieBbedingung  (5)  einge- 


schrankt . 


;p  ;r  r  ' p  ( *  '  <i  i  -  •  p\ 

J  =J  -Jl  -(y,,  -a 


dF  =  stat. 


In  Gleichung  (7)  wurde  bereits  beriicksicht igt ,  daB  nach  [5]  der 
Lagrange-Faktor  mit  dem  plastischen  Dehnungsi nkrement  e^  iden- 
tisch  ist.  Die  Euler-Gleichungen  des  Funktionals  liefern  u.a. 
die  konsti tutiven  Gleichungen  nach  Prandtl-ReuB .  Dies  beweist, 
daB  die  FlieBbedingung  (5)  der  Pragerschen  Normalenregel 


•  p  ~  ;p 

£;j  =  y.j  •  £ 


Equivalent  ist. 

Ausaehend  von  einem  bekannten  Anf angszustand  z  wird  das  Funktio¬ 
nal  (7)  finit  ubersetzt.  Es  entsteht  ein  algebraisches  Glei- 
chungssystem,  das  fiir  eine  kleine  endliche  Laststei gerung  den 
Zuwachs  z  der  ZustandsgroBen  beschreibt.  Mit  f ortschrei tendem 
FlieBen  weitet  sich  die  plastische  Zone  auf.  Die  dort  auftreten- 
den  plastischen  Dehnungsinkremente  werden  wie  die  ZustandsgroBen 
elementweise  approximi ert .  Um  das  teilweise  Plasti f izieren  eines 
Elementes  naherungswei se  zu  erfassen,  werden  bei  der  numerischen 
Integration  der  Arbei tsausdriicke  nur  die  GauB-Punkte  innerhalb 
der  plastischen  Zone  beriicksichtigt .  Die  mit  f ortschreitender 
Flastif izierung  sukzessive  auftretenden  Freiwerte  fiir  die  pla¬ 
stischen  Dehnungsinkremente  werden  fortlaufend  am  Ende  des 
Losungsvektors  angefiigt.  Dadurch  wird  das  vorhandene  Gleichungs- 
system  mit  der  FlieBbedingung  (5)  gerandert,  Bi Id  2 . 


Die  Ausgangsgleichungen  fiir  z  bleiben  dabei  unverandert  und  wer- 
ri n n  nur  ein  einziqes  trial  aufqeldst.  Die  Randgl ei chungen  werden 
fur  j  ede  Laststufe  iterativ  solange  verbessert,  bis  der  Span- 
nuii'ti'punkt  fiir  den  benachbarton  Zustand 


auf  der  FlieBflache  liegt. 

Da  bei  der  Iteration  nur  wenige  Randgleichungen  verandert  und 
neu  gelost  werden  miissen,  ist  die  Berechnung  sehr  schnell.  Sie 
hat  allerdings  den  Nachteil,  daB  die  Programmorganisation  fur 
eine  noch  unbekannte  Zahl  von  Randgleichungen  komplizierter  ist 
als  fur  ein  normales  Gleichungssystem  mit  festen  Grenzen. 

Als  Beispiel  fur  die  Berechnung  zeigt  Bi Id  3  die  Last-Verfor- 
mungs-Kurven  fur  den  oberen  Riegel  eines  Stockwerksrahmens .  Der 
EinfluB  der  Normalkrafte  auf  die  plastischen  Verformungen  der 
FlieBgelenke  ist  gering.  Sie  werden  bei  iiberwiegend  auf  Biegung 
beanspruchten  Tragwerken  meist  vernachlassigt . 

3.2  Die  Strukturnichtlinearitat 

Bei  einem  f ortschrei tenden  RiB  andert  sich  -  zunachst  lokal  - 

die  Geometrie  der  Struktur.  Mit  wachsender  RiBlange  treten  dann 

in  zunehmendem  MaBe  Spannungsumlagerungen  auf,  die  das  Tragver- 

halten  der  Struktur  verschlechtern  und  schlieBlich  zum  Bruch 

fiihren.  Sind  dabei  die  Verformungen  noch  hinreichend  klein,  kann 

man  sie  vernachlassigen .  Erst  der  RiBf ortschri tt,  der  das  stati- 

sche  System  verandert,  bringt  einen  nichtlinearen  Effekt  in  die 

Berechnung.  MaBgebend  fur  den  RiBf ortschri tt  ist  die  Spannungs- 

konzentration  an  der  RiBspitze,  die  durch  die  Spannungsintensi- 

tatsfaktoren  K  beschrieben  wird.  Uber schrei ten  sie  einen 
a 

kritischen  Wert  KC  oder  bei  kombinierter  Beanspruchung  die  In- 

a 

teraktionsbedingung 
(10)  f  lKa;  KCJ  £  0  , 

weitet  sich  der  RiB  aus.  Die  Schwierigkei t  bei  der  Berechnung 
der  Spannungsintensitatsf aktoren  besteht  darin,  daB  die  Span- 
nungssingulari tat  an  der  RiBspitze  durch  finite  Elemente  nur 
unvollkommen  approximiert  wird.  Im  folgenden  werden  verschiedene 


A  »A 


Methoden  zur  Berechnung  der  Spannungsintensi tatsf aktoren  erlau- 
tert . 


3.2.1  Extrapolationsmethode 

Fiir  die  unendlich  ausgedehnte,  gerissene  Scheibe  existieren 
exakte  Losungen  iiber  den  Verlauf  der  Spannungen  und  Verschiebun- 
gen.  Die  Spannungen  haben  an  der  RiBspitze  eine  Singulari tat, 
die  mit  1//r  abklingt,  Bild  4.  Die  singulare  Losung  erfiillt  bei 
einer  endlichen  Struktur  nicht  die  Randbedingungen,  so  daB  zur 
Korrektur  der  Randwerte  stetige  Funktionen  iiberlagert  werden. 

Die  singulare  Losung  wird  durch  die  Gleichungen 


Uct=  Ua  >?•  aij)  -v/^T-Kx 
o“P=  s“p(y,  Qjj )  /  V2m  -K  x 


X  =  1.2 


approximiert .  Darin  sind  und  bekannte  Funktionen,  die 

nur  vom  Richtungswinkel  0  und  den  Koef f i zienten  des  Stoffge- 

setzes  abhangig  sind. 

Fiir  0=0  wird  die  Losung  (11)  fiir  die  Querzug-  und  Schubspan- 
nungen  entlang  der  x^Achse  besonders  einfach: 


K ,  =  v  2  it  r  •  or ' 


K  2  =\//2~tt7  ■  a  12 


fur  ip  =  0 


Fiir  beliebige  Radien  r  liefern  d;.e  Gleichungen  (12)  fiir  die  sin¬ 
gulare  Losung  konstante  Spannungs intensitatsf aktoren.  Setzt  man 

CL  2 

daaegen  naherungswe i se  die  Spannungen  c  aus  einer  finiten  Ele- 
mentberechnung  ein,  sind  die  ermittelten  Spannungsintensi tats- 
faktoren  vom  Radius  abhangig  (untere  Linie  in  Bild  5 )  . 


Die  Ergebnisse  der  finiten  E lementberechnung  enthalten  nicht  nur 
den  singularen  Term,  sondern  auch  die  homogene  Losung  zur  Kor¬ 
rektur  der  Randwerte.  Da  die  homogene  Losung  im  gesamten  Gebiet 

stetig  und  schwach  veranderlich  ist,  muB  die  Funktion  K  (r)  na- 

a 

hezu  lineare  Bereiche  aufweisen.  In  der  Nahe  der  RiBspitze  tre- 
ten  Abweichungen  von  der  Geraden  auf,  weil  die  finite  Elementlo- 


i 


sung  zu  ungenau  ist.  Ahnlich  verhalt  es  sich  bei  groBen  Radien, 

da  die  singulare  Losung  schnell  abklingt  und  der  RandeinfluB  an 

Bedeutung  gewinnt.  Extrapoliert  man  den  geraden  Mittelbereich 

der  Funktion  K  (r)  bis  zur  RiBspitze,  erhalt  man  eine  Naherung 
o 

fur  die  Spannungsintensitatsfaktoren. 

Leider  ist  die  finite  Elementberechnung  haufig  so  ungenau,  daB 
die  Funktion  K  (r)  keine  erkennbare  Gerade  enthalt  (obere  Linie 
in  Bild  5).  Tests  an  unidirektionalen  Zugproben  aus  CFK  haben 
gezeigt,  daB  oft  auch  bei  extrem  feiner  Elementteilung  keine 
Extrapolationsgerade  entsteht.  AuBerdem  hat  die  Extrapolations- 
methode  den  Nachteil,  daB  der  Spannungsverlauf  gezeichnet  und 
per  Hand  ausgewertet  werden  muB .  Tafel  1  zeigt  ein  Beispiel  fur 
die  Berechnung  der  Spannungsintensitatsfaktoren  in  einer  gekerb- 
ten  und  angerissenen  Zugprobe  aus  unidi rektionalem  CFK-Laminat. 


3.2.2  Die  Methode  der  Linienintegrale 

Von  Rice  [6]  u.a.  wurden  Linienintegrale  angegeben,  die  unabhan- 
gig  vom  Integrationsweg  sind.  Das  sogenannte  J-Integral  von  Rice 
ist  ein  MaB  fur  die  Beanspruchung  an  der  RiBspitze.  Es  basiert 
auf  dem  Prinzip  der  virtuellen  Verriickungen  (1).  Im  folgenden 
sind  keine  Volumenbelastungen  zugelassen.  Die  an  den  RiBufern 
austretenden  Spannungen  sollen  verschwinden  und  die  Verschiebun- 
gen  sollen  die  Wegrandbedingungen  erftillen.  Ersetzt  man  auBerdem 
die  variierten  GroBen  durch  die  Differentials,  entsteht  die 
Gleichung 

dj  =  J,oIJ  *duijdV-J’  pp’dU|  dF 

(13)  v  F* 

=  \  /®U  ‘  rj  ’  d  AUj  d F  . 
dFc 

Darin  beschreibt  die  erste  Zeile  die  gesamte  potentielle  Energie 
der  Struktur  und  die  zweite  Zeile  die  bei  einem  RiBf ortschritt 
dF  dissipierte  Energie.  Die  dissipierte  Energie  laBt  sich  durch 
die  Energief rei setzungsrate  G  (Energy  release  rate)  multipli- 


< i 

I  =  -  d  Fc  •  G  I 


"i  v«  F.rf  ■  l  m  n.i 


ziert  mit  dem  Zuwachs  dFc  der  Rififlache  ersetzen.  Etwas  umge- 
formt,  lautet  Gleichung  (13) 

d  J  r  dW  ( -y  jj 


(14) 


a  j  f  aw  i  - 

-  =  G  =  -  /  - 1 

i  dF, 


dV  *  /  Pi  •  Hr-  dF«  ■ 

C  p  U  rc 

a 


Im  folgenden  wird  nur  der  ebene  Spannungszustand  betrachtet 
CL  3 

(o  =0).  An  die  Rifispitze  wird  ein  Koordinatensystem  angehef- 

1 

tet,  das  mit  der  Rifispitze  mitwandert  und  dessen  Koordinate  x 
in  die  Verlangerung  des  Risses  zeigt,  Bild  6. 

In  die  Gleichung  (14)  wird  die  Beziehung 

(15)  d Fc  =  -  d x  1 


fc**-  I**  '.*•  .“•*)» ‘ 

k>V-  *“*  V  •  V 


m 


eingefiihrt  und  das  Volumenintegral  durch  das  entsprechende  Rand- 
integral  ersetzt. 

(16)  J.,  --  j  (W  (y  .)  •  r,  -  pe  -u  ,i  dS 

s 

Der  Integrationsweg  in  Gleichung  (16)  mufi  nicht  dem  Rand  der 

wirklichen  Struktur  folgen,  sondern  er  kann  auch  eine  beliebige 

Teilstruktur  umfahren,  Bild  6.  Aufierdem  kann  der  Rifif ortschri tt 

2 

und  damit  die  Ableitung  in  Richtung  der  Koordinate  x  erfolgen. 
Damit  geht  die  Gleichung  (16)  in  die  allgemeinere  Form 

(17)  J.«  =  /  •  rtt  -  p £-ue  JdF 

s 

liber . 


In  das  J-Integral  (17)  wird  die  singulare  Losung  nach  Gleichung 
(11)  eingesetzt.  Die  Auswertung  liefert  eine  Beziehung  zwischen 
den  J-Integralen  und  den  Spannungsintensi tatsf aktoren . 


Die  Konstanten  Rp^  sind  nur  von  den  Koef f izienten  des  Stoffae- 

a  = 

setzes  abhangig.  Sie  nehmen  bei  isotropera  Material  die  folgen¬ 
den,  von  Null  verschiedenen  Werte  an: 


15 


Mit  der  singularen  Ldsung  fiir  den  geraden  RiB  wurde  auch  das 
Integral  J ,  ausgewertet,  bei  dem  der  RiBf ortschritt  rechtwih 
lig  zur  urspriinglichen  RiBrlchtung  verlauft.  Die  Gleichung  fii 
J,2  in  (18)  beschreibt  keine  mechanisch  sinnvolle  Losung.  Auf 
Grund  dieses  Widerspruches  konnen  die  beiden  Gleichungen  (18) 
und  (19)  auch  nicht  benutzt  werden,  um  die  Spannungsintensita 
faktoren  iiber  die  J-Integrale  zu  bestimmen.  Einen  Ausweg  aus 
diesem  Dilemma  bietet  die  Trennunq  der  Enerqief rei setzunqsrat* 


Bei  der  Ableitung  der  J-Integrale  (17)  wurde  vorausgesetzt ,  daB 
das  virtuelle  Verschiebungsf eld  6u^  affin  ist  zu  den  wirklichen 
Verschiebungen  (Su^wurde  durch  du^  ersetzt).  Es  gibt  daher 
keine  Moglichkei t ,  daB  J, ^-Integral  explizit  in  die  beiden  Terme 
und  aufzuspalten .  Im  folgenden  soil  gezeigt  werden,  daB 
diese  Trennung  naherungswei se  mit  Hilfe  der  finiten  Elementme- 
thode  moglich  ist. 

Bei  einem  linear  elastischen  Korper  ist  die  innere  Energie 

W(J  .  )  halb  so  groB  wie  das  Potential  der  auBeren  Krafte.  Unter 
P  " 

dieser  Voraussetzung  laBt  sich  die  Gleichung  (14)  umformen  in 


(24)  G  .  1/ p'  ^  dS  . 

2l  He 


Bei  dem  hier  betrachteten  ebenen  Spannungszustand  und  der  Schei- 

bendicke  d  =  1  wird  dF  durch  die  Lange  dl  des  RiBf ortschri tts 

c  c 

ersetzt.  Mit  den  in  Bild  7  angegebenen  Bezeichnungen  wird  die 
Gleichung  (24)  als  Dif f erenzenausdruck  fur  die  Einzelkrafte  an- 
geschrieben: 


(25) 


£P'- Aui 
2Ac 


mi  t 

AUj  =  Uj(lc  +Cn) -Uj  dc  ♦  Cn-1  )  , 

(26) 

Ac  -  Cp  -  Cn_i 

Darin  sind  P1  die  einer  RiBlange  1  zugeordneten  Belastungen. 
Diese  werden  bei  einer  virtuellen  Verlangerung  des  Risses  um  die 
Lange  c  konstant  gehalten.  Der  virtuelle  RiBf ortschri tt  c^  wird 
in  n  Stufen  aufgebracht,  wobei  Ac  iiblicherweise  einer  Element- 
lange  entspricht. 

Die  Gleichung  (24)  beschreibt  die  gesamte,  an  der  RiBspitze 
freigesetzte  Energie.  Wenn  fur  den  virtuellen  RiBf ortschri tt  nur 


eine  RiBoffnung  quer  zur  RiBachse  zugelassen  wird  und  die  Tan- 
gentialverschiebungen  der  angrenzenden  Elemente  weiterhin  gekop- 
pelt  bleiben,  liefert  die  Gleichung  (25)  nur  die  Energiefrei- 
setzungsrate  fur  .  Analog  wird  fur  nur  die  Schubverbindung 
zwischen  den  angrenzenden  Elementen  gelost,  wahrend  die  Element- 
knoten  quer  zum  RiB  gekoppelt  bleiben. 

Der  Vorteil  dieses  Verfahrens  besteht  darin,  daB  die  Auswertung 
vollstandig  im  Rechner  erfolgen  kann.  AuBerdem  ist  das  Ergebnis 
genauer  als  bei  der  Extrapolationsmethode  da  die  mittleren  Ener- 
giewerte  bei  der  finiten  Elementmethode  erheblich  genauer  appro- 
ximiert  werden  als  einzelne  Spannungskomponenten  in  der  Nahe  der 
RiBspitze . 

Bild  8  zeigt  ein  Beispiel  fur  die  Energief reisetzungsrate  eines 
gekerbten  [  90 0  ]  -.-Laminates  aus  CFK .  Die  Energief  rei setzung  nimmt 

O 

mit  wachsender  RiBlange  zu,  so  daB  in  diesem  Beispiel  der  RiB 
bis  zum  Bruch  der  Struktur  weiterlauft  (plotzlicher  Bruch).  Nur 
wenn  die  Energief reisetzung  mit  wachsender  RiBlange  abnimmt, 
stoppt  der  RiB,  sobald  ein  kritischer  Wert  GC  unterschri tten 
wird.  Die  J-Integrale  sind  quadratische  Funktionen  der  Span- 
nungsintensitatsf aktoren.  Die  Funktionen  Ka(r)  verlaufen  in  der 
Nahe  der  RiBspitze  linear,  Bild  5,  so  daB  die  Energief rei set- 
zungsrate  G  naherungswei se  durch  eine  quadratische  Funktion  der 
RiBlange  beschrieben  wird.  In  der  Regel  tragt  man  daher  die 
Funktion  /G  auf,  die  linear  extrapoliert  werden  kann. 

Eine  andere  Methode  zur  Bestimmung  der  Energief rei setzungsraten 
G^(K^)  und  G^fK^)  beruht  auf  den  in  (7)  angegebenen  Formeln 


(27) 


G, 


1 

=  lim  -=-r-  /  cr 22  ( r  =  x  ,  ip  =  0)- Au2(r=  6c-x,  »p  =  in)  dx 

6c  —  0  ‘oc  n 


g2  = 


lim 
6c— 0 


1 

26c 


6c 


y"  0 12  ( r  =  x,f  =  0)'  iu,(r  =  6c-x,<p  =  +  n)  dx 
o 


Die  zugeordneten  Bezeichnungen  sind  in  Bi Id  9  beschrieben.  Die 
Gleichungen  (27)  liefern  nur  fur  die  formelmaBig  angegebene, 
exakte  Losung  die  zugeordneten  Energi ef rei setzungsraten .  Die  so 


ermi ttelten  G  -Werte  sind  meist  sehr  ungenau,  wenn  fur  einen 
a 

endlichen  RiBfortschri tt  die  numerisch  ermittelten  Spannungen 
und  Verschiebungen  aus  einer  finiten  Elementberechnung  einge- 
setzt  werden.  Dies  liegt  an  der  i.a.  sehr  ungenauen  Approximate 
on  der  Spannungen  in  der  Nahe  der  RiBspitze.  In  [8]  wird  eine 
modifizierte  RiBschluBintegralmethode  (modified  crack  closure 
integral  method)  erlautert,  bei  der  die  in  Gleichung  (27)  ent- 
haltenen  Funktionen  fur  die  Spannungen  und  Verschiebungen  durch 
die  entsprechenden  Knotenkrafte  und  -verschiebungen  der  finiten 
Elementberechnung  ersetzt  werden.  Mit  den  Bezeichnungen  nach 
Bild  10  entstehen  fur  die  Energief rei setzungsraten  die  Differen 
zenausdrucke 

ue, 

G2=  —-(Fr.AuJ.F^.Au,3), 

die  nach  [8]  beim  WeggroBenverf ahren  nur  um  wenige  Prozent  von 
der  genauen  Losung  abweichen. 


3.2.3  Sinoulare  finite  Elemente 


Die  Genauigkeit  einer  finiten  Elementberechnung  laBt  sich  erheb 
lich  steigern,  wenn  die  Singularitat  an  der  RiBspitze  durch  be- 
sondere  Elemente,  sogenannte  singulare  Elemente,  erfaBt  wird. 
Fur  diese  Elemente  wird  als  Ansatz  die  genaue  Losung  der  Schei- 
bengleichungen  herangezogen .  Sie  konnen  daher  beliebig  groB 
sein.  Als  Beispiel  wird  hier  nur  das  Element  von  Hu  (10)  erlau¬ 
tert,  das  auf  der  Spannungsf unktion  von  Williams  [9]  aufbaut. 
Williams  Spannungsfunktion 

♦'I/5'’  {'V[s'n(?-')  '  V  "  nfl  '  s'n  (-J-  •  1  )  '  ?  ] 

*  bn  '  [C0S(  2  '  1  )  '  T  -  (  J  •  ')  •  f]  } 


(29) 


ist  eine  Reihenentwicklung,  deren  erstes  Glied  die  Singularitat 
an  der  RiBspitze  beschreibt.  Die  hoheren  Reihenglieder  sind  ste- 
tige  Funktionen,  mit  denen  die  Randbedingungen  erfiillt  werden 
konnen.  Aus  Gleichung  (29)  laBt  sich  das  Verschiebungsf eld 

(30)  ua  =  u£  (r,ip)  Kx  Fn(r,<p)  una 

n=  2 

ableiten,  das  als  Freiwerte  fur  die  finite  Ubersetzung  die  Span- 
nungsintensi tatsf aktoren  und  die  Verschiebungen  in  den 
Randknoten  enthalt.  Bild  11  zeigt  das  singulare  Element,  einge- 
bettet  in  ein  normales  finites  Elementnetz. 

Mit  den  singularen  Elementen  kann  ein  vorhandener  RiB  nachge- 
rechnet  werden.  Die  Untersuchung  eines  RiBf ortschri ttes  ist  nur 
mit  groBem  numerischen  Aufwand  moglich,  da  das  singulare  Element 
und  damit  das  Elementnetz  sukzessive  an  den  RiBf ortschri tt  ange- 
paBt  werden  miissen. 

3.  RiBbeschreibung  durch  Strukturvariable 


Damit  bei  f ortschrei tendem  RiB  das  Elementnetz  unverandert  blei- 
ben  kann,  werden  Strukturvariable  eingefiihrt.  Sie  sind  an  die 
Struktur  -  in  unserem  Fall  an  den  RiB  -  gebunden  und  nicht  ein- 
zelnen  Elementen  zugeordnet.  Fur  die  Strukturvariablen  werden 
folgende  Ansatzfunktionen  vorgeschl agen : 


Ua  =  Lj  (S.nl-Ua  * 

*  •  ( g 1 ' 2  -  g 2 )  •  U  i  ( <P ,  a  i  j )  •  K  x  ♦ 

+  Dm  ( 5c  ,  'I  c  )  ’  A  U U 


stetiger  Ansatz 

Singularitat  an 
der  Riflspitze 

RiHoffnung 


0  <  g  <  1 


0  <  <  1  • 


Darin  beschreiben  der  erste  Term  die  iiblichen  Ansatze  innerhalb 
ernes  finiten  Elementes  und  die  folgenden  Terme  die  Strukturva- 


r  m 

/{.w. 

.  *>  V  V 

<  v  s.  * 


h'.-.V.V 


>■ 

»  .  - 


v- 

->v< 


V  V  V  v 

It/ 


20 


riablen  fur  die  Singularitat  an  der  RiBspitze  und  fur  die  RiB- 
of fnung . 

Der  Ansatz  fur  die  Singularitat  enthalt  die  exakte  Losung  nach 
Gleichung  (11),  die  mit  stetigen  Funktionen  iiberlagert  wird. 
Dadurch  verschwindet  die  Ansatzfunktion  auf  einem  Kreis  mit  dem 

A 

Radius  R  .  Die  Freiwerte  K  sind  die  Spannungsintensitatsf akto- 
o  a 

ren.  Bild  13  zeigt  die  Ansatzfunktionen  fur  die  RiBspitze.  Bei 
der  finiten  Ubersetzung  werden  die  Arbei tsausdrucke  fur  die 
Strukturvariablen  der  RiBspitze  nur  in  den  Bereichen  der  Elemen- 
te  ausgewertet,  die  innerhalb  eines  Kreises  mit  dem  Radius  Rq  urn 
die  RiBspitze  liegen,  Bild  12.  Naherungswei se  werden  bei  einer 
numerischen  Integration  nur  die  GauBpunkte  eines  Elementes  aus¬ 
gewertet,  die  innerhalb  des  Kreises  Rq  liegen. 

In  Gleichung  (31)  wird  die  RiBoffnung  durch  die  Sprungfunktionen 

D  beschrieben,  die  in  Richtung  des  Risses  quadratisch  approxi- 
m 

miert  werden,  Bild  12.  An  der  RiBspitze  miinden  die  Sprungfunk¬ 
tionen  tangential  in  die  RiBachse  ein,  so  daB  dort  keine  Singu¬ 
laritat  auftritt. 

Die  Strukturvariablen  entstehen  sukzessive  mit  wachsender  RiB- 
lange.  Sie  werden,  wie  die  plastischen  Dehnungsinkremente  in 
Bild  2,  am  Ende  des  Losungsvektors  angeordnet.  Dadurch  randern 
die  Gleichungen  zur  Beschreibung  des  Risses  das  algebraische 
System.  Die  Randgleichungen  miissen  nur  im  Einf luBbereich  der 
singularen  Funktionen  iterativ  verbessert  werden,  da  nur  hier 
Veranderungen  durch  den  RiBf ortschri tt  eintreten. 

Die  Genauigkeit  der  beschri ebenen  Verfahren  ist  unterschiedlich . 
Bild  14  zeigt  den  Vorteil,  den  die  singularen  Elemente  bieten, 
besonders  deutlich.  Sie  liefern  bei  wesentlich  groBeren  Element- 
abmessungen  genauere  Ergebnisse  als  die  ersten  beiden  hier  be- 
schriebenen  Verfahren.  Die  Methode  der  Strukturvariablen  diirfte 
den  singularen  Elementen  gleichwertig  sein. 
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3.3  Die  qeometri sche  Nichtlineari tat 

Die  bisher  beschriebenen  physikalisch  nichtlinearen  Probleme 
lassen  sich  fur  die  Inkremente  lineari sieren .  Dadurch  konnen  sie 
mathematisch  exakt  uber  Randgleichungen  erfaflt  werden.  Anders 
liegen  die  Verhaltnisse  bei  stark  verformten  Strukturen.  Hier 
mussen  die  Verformungen  im  Gleichgewicht  berucksichtigt  werden. 
Dadurch  andert  sich  in  der  finiten  Ubersetzung  der  Struktur  die 
Koef f izientenmatrix .  Bei  physikalisch  und  geometrisch  nichtli¬ 
nearen  Froblemen  ist  es  daher  i.a.  nicht  sinnvoll,  die  Werk- 
stof fnichtlinearitat  uber  Randgleichungen  zu  erfassen,  wenn 
auBerdem  aufgrund  der  groBen  Verformungen  die  Koef f izienten¬ 
matrix  geandert  werden  muB.  Urn  auch  hier  die  Methode  der 
Randgleichungen  anwenden  zu  konnen,  wird  im  folgenden  ein  Nahe- 
rungsverf ahren  vorgeschlagen  und  am  Beispiel  eines  einfachen 
Balkens  erlautert. 

Fur  den  in  Bild  15  skizzierten  Balken  soil  der  Anf angszustand 
rnit  dem  Ldsungsvektor  uq  und  die  zugeordnete  finite  Ubersetzung 
bekannt  sein.  AuBerdem  soil  der  Anf angszustand  die  Gleichge- 
wichtsbedingungen  erfullen.  Das  Gleichungssystem 

(33)  SuT  (  A  I  u„)  •  u  -  p)=  0 

liefert  fur  eine  kleine  endliche  Laststeigerung  p  den  zugeordne- 
ten  Losungsvektor,  der  als  1 .  Glied  einer  Reihenentwicklung  mit 
u^  bezeichnet  wird.  Aufgrund  der  uberpropcrtional  anwachsenden 
Verformungen  erfullen  die  inneren  Krafte  im  System  fur  den  Nach- 
barzustand 

(34)  ut  =  u0 ♦ u, 

nicht  die  Gleichgewichtsbedingungen .  Sie  halten  nur  dem  Lastan- 
tei  1 

(35)  b  =  A ( u , )  •  u , 


die  Waage .  In  den  folgenden  Schritten  werden  jeweils  die  Last- 
differenzen  p  -  b  aufgebracht  bis  die  Losung  konvergiert.  Fur 


die  folgende  Laststufe  ist  der  berechnete  Nachbarzustand  der 
neue  Anfangswert. 

Ublicherweise  wird  die  Koef f izientenmatrix  zu  Beginn  einer  neuen 
Laststufe  angepaBt  und  die  Iteration  innerhalb  einer  Laststufe 
iiber  die  Lastspalten  durchgefiihrt .  Im  Gegensatz  dazu  soil  hier 
die  Losung  durch  eine  Reihenentwicklung  angenahert  werden, 

Bild  16. 

•  *  )|r  • 

(36)  u  =  X1-u1  +  u2 
mi  t 

(37)  u,  =  u ,  /  II u ,  II  - 

Darin  sind  die  neue  Variable,  u^  die  Losung  des  ersten  Itera- 
tionsschri ttes  und  ||*||  die  Euklidische  Norm.  Setzt  man  diese 
Losung  in  das  Funktional  (33)  ein,  entsteht  das  in  Bild  16  skiz- 
zierte  geranderte  Gleichungssystem .  Dessen  Losung  liefert  in 
einem  Schritt  den  zum  Gleichgewicht  erf order lichen  Anteil 
u1  =  X  *  u,|  .  Der  Vektor  u 2  ist  zu  u1  orthognonal  und  korrigiert 
die  Losung  u^.  In  den  folgenden  I terationsschri tten  wird  die  Rei¬ 
henentwicklung  (36)  jeweils  urn  ein  Glied  verlangert,  bis  die 
Losung  konvergiert  (u^  -*  0).  Zur  Verbesserung  der  Konvergenz 
wird  nach  einem  Vorschlag  von  Crisfield  [11]  die  Belastung  so 
reduziert,  daB  die  Nachbarpunkte  auf  einem  Kreisbogen  liegen. 

Das  vorgeschlagene  Naherungsverf ahren  beruht  auf  der  Annahme, 
daB  die  Beulformen  durch  die  Uberlagerung  weniger  orthogonaler 
Biegelinien  beschrieben  werden  konnen.  Numerische  Untersuchungen 
mussen  zeigen,  ob  diese  Annahme  gerechtf ertigt  ist. 


4.  Tei lstrukturtechnik 


4.1  Anf orderungen  an  den  Gleichunqsloser 

An  einen  Gleichungsloser  werden  sowohl  vora  mathematisch-mechani- 
schen  Modell  als  auch  von  der  Programmorganisation  her  bestimmte 
Anforderungen  gestellt.  Fur  die  hier  anstehenden  Probleme  sind 
die  wichtigsten  im  folgenden  aufgelistet  und  erlautert: 

•  Gleichungssystem  symmetrisch: 

Bei  der  Losung  symmetri scher  Gleichungssysteme  wird  nur  auf 
die  untere  Halfte  der  Koef f izientenmatrix  zugegriffen.  Da- 
durch  wird  der  Rechen-  und  Speicheraufwand  gegeniiber  einem 
unsymmetri schen  System  etwa  auf  die  Halfte  reduziert.  Diese 
Forderung  laflt  sich  bei  geeigneter  Wahl  des  Funktionals  fast 
immer  erfiillen. 

•  Gleichungssystem  semidefinit: 

Bei  den  gemi schen  Methoden  sind  die  entstehenden  Gleichungs¬ 
systeme  semidefinit,  d.h.  bei  der  Dreieckszerlegung  der  Ko¬ 
ef  f  izientenmatrix  konnen  Pivot-Elemente  mit  PIV  <  0 
auftreten.  Eine  Pivot-Suche  und  die  damit  verbundene 
Spalten-  und  Zei lenvertauschung  muB  moglichst  vermieden  wer¬ 
den,  urn  aufwend^ge  Indexrechnungen  und  Speicherumordnungen 
zu  vermeiden. 

•  Randwerte  variabel: 

Bei  den  anstehenden  Problemen  wachsen  die  Randwerte  nicht 
nur  proportional  an,  sondern  sie  konnen  sich  auch  unabhangig 
voneinander  in  den  einzelnen  Laststufen  andern.  Daher  miissen 
die  in  der  Koef f izientenmatrix  zu  streichenden  Spalten  auf- 
gehoben  werden,  damit  sie  bei  einem  spateren  Restart  mit 
veranderten  Randwerten  wieder  zur  Verfiigung  stehen. 

•  Speicherung  zeilenweise: 

Die  untere  Halfte  der  Koef f izientenmatrix  wird  zeilenweise 
in  einem  eindimensionalen  Feld  abgelegt.  Die  Speicherung 
einer  Zei le  beginnt  beim  ersten  besetzten  Element  und  endet 
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mit  dem  Hauptdiagonalenelement  ( MD-Element ) . Um  unnotige  In- 
dexrechnungen  zu  vermeiden  und  um  moglichst  lange  Vektoren 
fur  die  bei  der  Auflosung  entstehenden  Skalarprodukte  zu 
erhalten,  muB  das  Gleichungssystem  zeiienweise  gelost 
werden. 

•  Blocken  des  Gleichungssystems : 

Bei  groBen  Gleichungssystemen  paBt  die  Koef f izientenmatrix 
nicht  mehr  vollstandig  in  den  Kernspeicher  hinein.  Daher 
werden  aufeinanderfolgende  Zeilen  zu  Zei lenblocken  zusammen- 
gefaBt  und  auf  dem  Hintergrundspei cher  abgelegt.  Um  die  An- 
zahl  der  zei taufwendigen  Zugriffe  auf  den  Hintergrundspei- 
cher  klein  zu  halten,  mtissen  die  Zeilenblocke  moglichst  groB 
gewahlt  werden. 

•  Verandern  des  Gleichungssystems : 

Die  anstehenden  nichtlinearen  Probleme  lassen  sich  nur 
schrittweise  iterativ  losen.  Dabei  miissen  teilweise  modifi- 
zierte  Gleichungssysteme ,  die  bei  der  vorqeschlagenen  Rande- 
rungsmethode  eine  variable  Anzahl  von  Unbekannten  aufweisen, 
mehrfach  gelost  werden. 

In  den  folgenden  Abschnitten  wird  ein  Glei chungsloser  vorge- 
stellt,  der  diesen  Forderungen  weitgehend  entspricht.  AuBerdem 
lassen  sich  damit  direkt  die  reduzierten  Elementmatrizen  und 
Lastspalten  fiir  die:  Tei  1  strukturen  aufbauen. 

4.2  Der  Aufbau  von  Tei 1 strukturen 

4.2.1  _Er lauterungen  zur  Tei 1 s trukturtechnik 

GroBe  Strukturen  werden  haufig  in  Tei  1  systems  zerlegt  und  iiber 
die  gemeinsamen  Koppelknoten  miteinanaer  verbunden.  Dadurch  ent- 
steht  bei  einer  finiten  Ubersetzung  nicht  ein  einziges  groBes 
Gleichungssystem,  sondern  es  bilden  sich,  entsprechend  der  An¬ 
zahl  der  Tei 1 strukturen ,  mehrere  kleine.  Die  Teilsysteme  konnen 
unabhangig  voneinander  berechnet  und  auf  dem  Hintergrunaspeicher 


abgelegt  werden.  Von  dort  werden  sie  bei  Bedarf  wie  ein  normales 
Element  abgerufen  und  zu  der  gewunschten  Struktur  zusammenge- 
setzt . 


Besonders  vorteilhaft  ist  die  Tei lstrukturtechnik  bei  der  Losung 
von  physikalisch  nichtlinearen  Problemen.  Hier  interessieren 
haufig  nur  die  Spannungen  und  Verschiebungen  in  der  unmittelba- 
ren  Umgebung  der  gestorten  Zone.  Als  "gestort"  wird  hier  der 
Bereich  einer  Struktur  bezeichnet,  in  dem  der  linear  elastische 
Spannungszustand  durch  Zwangsbedingungen  (FlieB-,  Bruch-,  Reib-, 
Kontaktbedingungen  u.a.)  eingeschrankt  wird.  Nur  dieser  Bereich 
wird  entsprechend  dem  Storungsf ortschri tt  jedesmal  neu  approxi- 
miert,  wahrend  die  restlichen  Tei lstrukturen  nur  einmal  berech- 
net  und  dann  nur  noch  zugeladen  werden. 

Bild  17  zeigt  ein  Versuchsstiick  aus  CFK  mit  einer  kleinen  Dela¬ 
mination.  Nur  die  unmittelbare  Umgebung  der  geschadigten  Zone 
wird  durch  aufwendige  Volumenelemente  approximi ert .  Die  restli- 
che  Struktur  wird  durch  einfache  Schalenelemente  angenahert.  Sie 
ubertragt  die  durch  die  Struktursteif igkeit,  Randbedingungen  und 
Belastungen  hervorgeruf enen  Beanspruchungen  auf  den  zu  untersu- 
chenden  Teilbereich. 

Der  in  Bild  17  schraffierte  Teilbereich  wird  durch  eine  Sub- 
struktur  ersetzt,  in  der  als  Freiwerte  nur  noch  die  Unbekannten 
in  den  Koppelknoten  und  evtl.  in  einigen  Randknoten,  in  denen 
veranderliche  Belastungen  eingeleitet  werden,  auftreten.  Die 
Unbekannten  in  den  hier  durch  schwarze  Punkte  markierten  Innen- 
knoten  werden  eliminiert.  Die  finite  Ubersetzung  der  Teilstruk- 
tur  liefert  das  unten  dargestellte  Gleichungssystem,  in  dem  nach 
der  Elimination  der  inneren  Variablen  z^  das  reduzierte  Element- 
system 

(38)  Ke?p  =  Pe 

entsteht.  Bevor  die  Berechnung  der  reduzierten  El ementmatri zen 
mathematiscl  abgeleitet  wird,  soil  die  vorgeschl agene  Substruk- 
turtechnik  anschaulich  erlautert  werden. 
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Angenommen,  das  in  Bild  17  dargestellte  Gleichungssystem  ist 
regular.  Dann  ist  die  normale  Dreieckszerlegung  der  Koeffizien- 
tenmatrix  und  die  Vorreduktion  des  Lastvektors  moglich.  Schnei- 
det  man  nun  vom  hinteren  Teil  des  zerlegten  Gleichungssystems 
ein  beliebiges  Teilsystem  ab,  enthalt  dieses  die  bereits  in 
Dreiecksmatrizen  zerlegte  reduzierte  Elementmatrix  Kg  und  die 
mit  der  zerlegten  Elementmatrix  vorreduzierten  Lastspalten  pe- 

Der  Gleichungsloser  ist  of f ensichtlich  einen  Schritt  zu  weit 
gegangen : 

Er  hat  nicht  nur  die  reduzierte  Elementmatrix  aufgebaut, 
sondern  sie  auch  in  Dreiecksmatrizen  zerlegt. 

Es  ist  nun  sehr  einfach,  den  Gleichungsloser  zu  modif izieren,  so 
daB  direkt  die  reduzierten  Elementmatrizen  und  die  reduzierten 
Lastspalten  entstehen.  Voraussetzung  ist  nur,  daB  im  Losungsvek- 
tor  zuerst  die  zu  eliminierenden  inneren  Variablen  angeordnet 
weraen  und  dann  erst  die  verbleibenden  externen  Variablen. 

Zum  Aufbau  der  Gesamtstruktur  werden  die  reduzierten  Element¬ 
matrizen  und  die  Lastspalten  aus  den  einzelnen  Teilsystemen 
"herausgef ischt"  und  in  das  Gleichungssystem  fur  die  Gesamt¬ 
struktur  eingesetzt.  Die  Losung  dieses  Gleichungssystems  liefert 

die  externen  Variablen  z  ,  die  bei  der  Ruckreduktion  in  der  Sub- 

e 

struktur  an  die  Stelle  der  Elemente  p  treten.  Die  normale 

e 

Ruckreduktion,  beginnend  mit  dem  letzten  Element  der  inneren 
Variablen  z^,  liefert  den  kompletten  Losungsvektor .  Hierbei  wer¬ 
den  die  externen  Variablen  zg  wie  inhomogene  Randbedingungen 
behandelt . 

Der  groBe  Vorteil  des  hier  vorgeschlagenen  Verfahrens  besteht 
darin,  daB  die  Substrukturen  ohne  besonderen  Mehraufwand  direkt 
beim  Loser,  der  Gleichungen  aufgebaut  werden.  Als  Beispiel  zeigt 
Tafe'  .  die  Losung  eines  semidef initen  Gleichungssystems,  in  dem 
die  ersten  m  Unbekannten  das  Elementsystem  fur  eine  Teilstruktur 
reprasentieren.  Das  Gleichungssystem  wurde  so  gewahlt,  daB  bei 
der  Auflosung  nur  aanze  Zahlen  auftreten,  die  leicht  kontrol- 
liert  werden  konnen.  Die  Gleichungen  wurden  nach  dem  Verfahren 
von  Cholesky  [12]  gelost,  das  als  bekannt  vorausgesetzt  wird. 


Bei  semidef ini ten  Systemen  treten  bei  der  Auflosung  auch  negati¬ 
ve  Pivot-Elemente  auf.  Dadurch  entstehen  in  der  Dreiecksmatrix 
auch  imaginare  Spaltenvektoren.  Da  bei  der  Zerlegung  nur  Produk- 
te  mit  rein  reellen  oder  rein  imaginaren  Zahlen  auftreten,  muB 
in  den  Skalarsummen  nur  das  richtige  Vorzeichen  beriicksichtigt 
werden. 

Bei  der  Auflosung  des  in  Tafel  2  angegebenen  Gleichungssystems 
entsteht  in  der  Zeile  i  =  3  eine  Pivotnull.  Um  eine  Pivotsuche 
zu  vermeiden,  wird  das  Hauptdi agonalelement  A^  willkurlich  um 
die  Zahl  1  vergroBert  und  die  Anderung  der  Ausgangsgleichungen 
durch  die  angefiigte  Randgleichung  n  wieder  riickgangig  gemacht. 
Beginnt  man  die  Dreieckszer legung  mit  der  Randgleichung,  laBt 
sich  leicht  zeigen,  daB  die  hinzugefugte  Diade  die  Ausgangsglei¬ 
chungen  nicht  beeinfluBt. 

Betrachtet  man  in  Tafel  2  die  letzten  4  Unbekannten  des  Glei¬ 
chungssystems  als  externe  Variablen,  laBt  sich  mit  den  zerlegten 
Teilmatrizen  das  in  Tafel  3  angegebene  reduzierte  Elementsystem 
bestimmen.  Zur  Kontrolle  wird  in  Tafel  4  das  reduzierte  Element- 
system  durch  blockweises  Auflosen  des  Gleichungssystems  berech- 
net . 


4.2.2  Losungsverf ahren  zum  Aufbau  von  Tei 1 strukturen 


Um  den  Beweis  der  propagierten  Substrukturtechnik  zu 
erleichtern,  wird  vorab  auf  einen  besonderen  Aspekt  der  Rand- 
gleichungsmethode  eingegangen.  Das  in  Bild  2  erlauterte  Randern 
des  Gleichungssystems  laBt  sich  nicht  ohne  weiteres  in  bestehen 
de  finite  Elementprogramme  implementieren,  da  diese  meist  von 
festen  GleichungsgroBen  ausgehen.  Damit  die  Methode  dennoch  ge- 
nutzt  werden  kann,  mtissen  die  Randgleichungen  liber  zusatzliche 
Lastspalten  erfaBt  werden. 

Ausgangspunkt  dieser  Untersuchung  ist  das  in  Bild  18  angegebene 
Gleichungssystem  A  *  z  =  r  ,  dessen  Variable  durch  n  Nebenbedin' 
gungen  eingeschrankt  werden.  Der  Lagrange- Faktor  fur  die  Neben- 
bedingungen  ist  z.  Die  erste  Matrizengleichung  wird  nach  z  auf- 


gelost,  wobei  die  Koef f izientenmatrix  A  durch  das  Produkt  der 

T 

beiden  Dreiecksmatrizen  L  *  L  ersetzt  wird.  Die  Losung  von  z 

wird  in  die  zweite  Matrizengleichung  eingesetzt.  Dabei  entsteht 

ein  Gleichungssystem  fur  den  Lagrange-Faktor  t,  das  nur  Terme 
T  - 1  T  - 1 

der  Form  g  *  L  bzw.  r  *  L  enthalt.  Diese  Terme,  in  Bild  18 

mit  v  und  y  bezeichnet.  entstehen  auch  bei  einer  normalen  Vor- 
r  J  g 

reduktion  des  Gleichungssystems  A  *  z  =  r,  wenn  die  Lasrspalten 
fur  r  urn  zusatzliche  Spalten  fur  die  Matrix  g  erweitert  werden. 


Mit  den  Losungen  y  und  y  wird  das  Gleichungssystem  fiir  e  auf- 

r  g 

gebaut,  wobei  jetzt  nur  noch  eine  rechte  Seite  auftritt.  Dieses 
Gleichungssystem  ist  i.a.  recht  klein,  so  daB  es  direkt  im  Kern- 
speicher  gelost  werden  kann.  Mit  der  Losung  e  wird  die  vorredu- 
zierte  Lastspalte  y^  modifiziert,  so  daB  die  Ruckreduktion  des 
urspriinglichen  Gleichungssystems  A  *  z  =  r  den  Losungsvektor  z 
lief ert . 


In  dieser  Form  kann  die  Methode  der  Randgleichungen  ohno  r>n- 
nenswerten  Mehraufwand  auch  in  normalen  finiten  Element prc-r?  air¬ 
men  realisiert  werden.  Mit  Hilfe  der  in  Bild  18  gezeigten  ;; 
chungen  wird  die  in  Bild  19  vorgeschlagene  Substrukturtecnn ;  r. 
erlautert.  Damit  beide  Bilder  direkt  vergleichbar  sind,  warden 
dieselben  Bezeichnungen  gewahlt.  Der  Vektor  z  enthalt  die  r,u 
el iminierenden  inneren  Variablen  und  e  die  verbleibenden  exter- 
nen  Variablen. 


in  Bild  19  liefert  die  normale  Zerlegung  des  Gleichungssystems 

A  *  z  =  r  die  Dreiecksmatrix  L  und  die  vorreduzierte  Lastspalte 

y  .  Die  Vorreduktion  der  Lastspalte  r  lauft  nach  den  gleichen 

Rechenvorschrif ten  ab,  wie  die  Bildung  einer  Spalte  der  oberen 

Dreiecksmatrix  L.  Daher  kann  sich  die  normale  Zerlegung  auch  auf 

die  Matrix  g  erstrecken,  so  daB  in  der  oberen  Nebendi agonalma- 

trix  y  entsteht.  Da  bei  einer  Auflosung  nach  Choleskv  die 
9 

Symmetric  der  Gleichungen  erhalten  bleibt,  taucht  die  transpo- 
T 

nierte  Matrix  y  auch  in  der  unteren  Nebendiagonalmatrix  auf. 


Die  zweite  Gleichung  in  Bild  18  enthalt  nur  e  als  Unbekannte. 
Damit  steht  dort  bereits  das  reduzierte  Elementsystem  fur  die 
externen  Variablen  e.  Analog  wird  in  Bild  19  die  reduzierte  Ele- 


mentmatrix  aufgebaut,  indem  von  der  bereits  vorhandenen  Matrix  h 

T 

das  Skalarprodukt  y  *  y  abgezogen  wird. 

Die  reduzierte  Elementmatrix  entsteht  auch  bei  einem  normalen 
Gleichungsloser ,  wenn  dort  die  Skalarprodukte  nicht  iiber  die 
gesamte  Zeilenlange  sondern  nur  tiber  den  Bereich  der  inneren 
Variablen  z  ausgefuhrt  werden  und  wenn  alle  Operationen  mit 
Hauptdiagonalelementen  im  externen  Bereich  unterbleiben . 
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Das  Losungschema  fiir  die  Lastspalten  ist  identisch  mit  dem  fur 
die  obere  Dreiecksmatrix ,  so  daB  hier  die  gleichen  Einschrankun- 
gen  gelten. 

Die  letzte  Gleichung  in  Bild  18  zeigt  die  modi f izierten  Last¬ 
spalten  fur  die  Rtickreduktion.  Sie  entstehen  mit  dem  normalen 
Gleichungsloser  autoraatisch,  wenn  der  reduzierte  Lastvektor  mit 
den  jetzt  bekannten  externen  Variablen  e  iiberspeichert  wird  und 
bei  der  Ruckreduktion  alle  Operationen  mit  Elementen  der  redu- 
zierten  Elementmatrix  unterdriickt  werden. 

Der  normale  Gleichungsloser  wird  geringfugig  modifiziert,  so  daB 
damit  auch  Substrukturen  aufgebaut  werden  kdnnen.  Selbstver- 
standlich  konnen  mit  der  modif izierten  Routine  die  Gleichungssy- 
steme  auch  vollstandig  gelost  werden,  indem  die  Anzahl  der 
externen  Variablen  Null  gesetzt  wird. 


Modifizierte  Cholesky  Routinen 


Im  Rechner  wird  die  untere  Halfte  der  Koef f izientenmatrix  zei- 
lenweise  gespeichert.  Urn  unnotige  Indexrechnungen  zu  vermeiden, 
werden  die  Gleichungen  auch  zeilenweise  aufgelost.  Bei  der  Re- 
duktion  der  einzelnen  Elemente  wird  nur  auf  zuriickl iegende ,  be¬ 
reits  zerlegte  Zeilen  zurtickgegrif  fen.  Vor  diesem  Element  ist 
die  unberuhrte  Matrix  vorhanden.  Dadurch  konnen  nach,  aber  auch 
wahrend  der  Dreieckszerlegung  beliebige  Randgleichungen  hinzuge- 
fugt  werden. 


Bei  dem  hier  vorgeschlagenen  Verfahren  miissen  fur  die  Randbedin- 
gungen  nicht  mehr  die  Zeilen  und  Spalten  in  der  Koef f izienten- 
matrix  gestrichen  und  zur  Lastspalte  addiert  werden.  Fur  die 
homogenen  Randbedingungen  werden  die  entsprechenden  Zeilen  und 
Spalten  beim  Auflosen  der  Gleichungen  iibersprungen.  Bei  den  in- 
homogenen  Randbedingungen  werden  die  zugeordneten  Zeilen  und 
Spalten  in  der  Koef f izientenmatrix  teilweise  mit  zerlegt,  so  daB 
sie  beim  Auflosen  der  Lastspalten  liber  die  Skalarprodukte  auto- 
matisch  zum  Losungsvektor  addiert  werden.  Dadurch  konnen  die 
Randwerte  auch  nach  der  Zerlegung  der  Koef f izientenmatrix  belie- 
big  geandert  werden,  da  die  zu  streichenden  Spalten  in  der 
Koef f izientenmatrix  erhalten  bleiben.  Dies  ist  bei  vielen  fini- 
ten  Elementprogrammen  nicht  der  Fall. 

Die  bei  einer  finiten  Elementberechnung  entstehenden  Gleichungs- 
systeme  sind  i.a.  zwar  symmetrisch,  aber  nicht  immer  positiv  de- 
finit.  Um  bei  einer  Pivotnull  eine  Zeilen-  und  Spaltenver- 
tauschung  zu  vermeiden,  wird  die  in  Abschnitt  4.2.1  (Tafel  2) 
erlauterte  Randerungsmethode  angewendet.  Da  sich  auch  die  Nicht- 
linearitat  liber  Randgleichungen  erfassen  last,  kommt  es  auf  eine 
zusatzliche  Randgleichung  fur  die  Pivotnull  nicht  an.  Die  prak- 
tische  Erfahrung  zeigt  auCerdem,  daB  Pivotnullen  nur  auBerst 
selten  auftreten. 

Um  die  einzelnen  Falle  beim  Auflosen  des  Gleichungssystems  un- 
terscheiden  zu  konnen,  wird  in  einem  Vektor  S  die  Pivotcharakte- 
ristik  abgespeichert ,  Bild  20.  Vor  der  Zerlegung  der  Koeffizien- 
tenmatrix  enthalten  die  Elemente  des  Vektors  S  die  Werte  von 
1-4.  Damit  werden  die  normalen  Gleichungen,  die  homogenen  und 
die  inhomogenen  Randbedingungen  gekennzeichnet .  S^  =  2  tritt 
nur  bei  wiederholter  Losung  des  Gleichungssystems  auf  und  kenn- 
zeichnet  die  Randgleichungen  zur  Korrektur  der  Pivotnullen.  Nach 
der  Zerlegung  der  Koef f izientenmatrix  enthalt  S  die  Zahlen  von 
5-8,  wobei  das  Vorzeichen  von  mit  dem  Vorzeichen  des  Pivot- 
elementes  identisch  ist. 


Nach  der  Auf losung  des  Gleichungssystems  kann  eine  beliebige 
Anzahl  von  Pivotkennzeichen  auf  die  urspriinglichen  Werte  zurlick- 
gesetzt  werden.  Das  soil  hier  z.B.  fur  die  Zeilen  N  und  1  ae- 
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schehen  sein.  Bei  einer  erneuten  Zerlegung  der  Koef f izienten- 
matrix  werden  dann  nur  die  Zeilen  mit  den  zuriickgesetzten 
Kennzeichen  und  der  zugeordnete  Hauptdiagonalblock,  beginnend 
mit  der  kleinsten  MD-Nummer,  neu  aufgelost.  Dadurch  kann  z.B. 
eine  iterative  Verbesserung  der  Randgleichungen  mit  minimalen 
Rechenaufwand  im  Gleichungsloser  erfaBt  werden.  In  den 
Bildern  21  -  23  sind  die  FluBdiagramme  fur  die  modif izierten 
Cholesky  Routinen  angegeben.  Aus  Griiden  der  Klarheit  wurde  in 
den  FluBdi agrammen  auf  eine  Blockung  der  Gleichungssysteme  fur 
die  Substrukturtechnik  verzichtet.  Tafel  5  zeigt  ein  numerisches 
Beispiel  fur  die  Substrukturtechnik. 


5.  Proaramm-  und  Datenorqani sation 


5 . 1  Glob a  1  e  _  Speicherorganisation 

Die  Methcde  der  Randgleichungen  hat  bei  den  hier  angesprochenen 
nichtl inearen  Problemen  gegeniiber  anderen  Verfahren  erhebliche 
numerische  Vorteile.  Da  die  Anzahl  der  Randgleichungen  vorab 
nicht  bekannt  ist,  setzt  sie  eine  flexible  Speichertechnik  vor- 


Im  Kernspeicher  werden  alle  anfallenden  Daten  in  einem  eindimen- 
sionalen  Feld  abgelegt,  Bild  24.  Aus  diesem  Feld  werden  bei  Be- 
darf  beliebige  Teilfelder  herausgeschni tten,  die  durch  ihre  An- 
fangsadresse  und  durch  ihre  GroBe  festgelegt  sind.  Beide  Werte 
und  bei  mehrdimensionalen  Feldern  auch  die  Anzahl  der  Matrizen- 
zeilen  und  -spalten  werden  im  Adressenfeld  KON  abgelegt, 

Bx Id  25 .  Alle,  die  Struktur  beschre i benden  Daten  und  Steuerfel- 
der  bleiben  permanent  im  Kernspeicher  und  sind  dort  in  dichter 
Folge  gespeichert.  Da  diese  Felder  zum  Teil  erst  mit  fortschrei- 
tender  Storung  (Schadigung)  der  Struktur  entstehen  oder  sich 
vergroBern,  schlieBt  sich  an  die  Steuerfelder  ein  leerer  Spei- 
cherpuffer  an.  Der  verbleibende  Platz  im  Speicherfeld  wird  in 
zwei  gleichgroBe  Arbeitsspeicher  aufgeteilt.  Die  GroBe  der  Ar- 
bei tsspeicher  bestimmt  die  maximal  mdgliche  BlockgroBe  in  den 
G1 ei chungs sy s t emen . 
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Die  Gleichungssysteme  werden  blockweise  aufgebaut,  gelost  und 
auf  den  Hintergrundspeicher  abgelegt.  Nach  AbschluB  einer  Last- 
stufe  werden  auch  die  zugeordneten  Steuerfelder  auf  den  Hinter¬ 
grundspeicher  geschrieben.  Parallel  dazu  wird  die  Adresse  fur 
das  erste  Element  der  Steuerfelder  und  der  erforderliche  Spei- 
cherbedarf  in  das  Adressenfeld  LOST  eingetragen,  das  der  exter- 
nen  Speicherung  vorangestellt  ist. 

In  den  folgenden  Laststufen  werden  nur  die  sich  andernden  Glei- 
chungsblocke  modifiziert  und  am  freien  Ende  des  Hintergrundspei- 
chers  abgelegt.  Die  neue  Speicheradresse  dieser  Blocke  wird  in 
die  im  Kernspeicher  vorhandenen  Steuerfelder  emgetragen,  die  am 
Ende  dieser  Laststufe  wieaer  in  den  Hintergrundspeicher  ge¬ 
schrieben  werden. 

Bei  der  vorgeschlagenen  Speicherorganisation  ist  ein  Restart 
nach  jeder  beliebigen  Laststufe  moglich,  indem  die  zugeoraneten 
Steuerfelder  eingelesen  werden.  Darin  sind  die  aktuellen  Adres- 
sen  fur  die  Gleichungsblocke  emgetraoen,  die  dieser  Laststufe 
zugrunde  liegen.  Die  Daten  fur  die  hoheren  Laststufen  werden  im 
Verlauf  der  Neuberechnung  uberschneben . 

5 .2 _ Dateneingabe  zur  Elementkopplung 

Die  vorgeschlagene  Substrukturtechnik  setzt  voraus,  daB  die  ex- 
ternen  Variablen  einer  Teilstruktur  immer  am  Ende  des  Losungs- 
vektors  stehen.  Daher  miissen  bei  einer  tiefen  Tei  1  struktur- 
schachtelung  die  externen  Unbekannten  mehrfach  umaeordnet 
werden.  Im  folgenden  wird  die  Eingabe,  die  die  Umordnung  der 
Variablen  steuert,  in  groben  Zugen  skizziert. 

In  den  KOPL-Karten  werden  die  Elemente  einer  Struktur  einander 
zugeordnet,  ohne  daB  bereits  interne  Bindungen  beriicksi  chti  gt 
werden  miissen.  In  Bf  Id  2  6  ist  em  Beispiel  angegeben,  in  dera  aus 
einem  Einzelelement,  beschneben  durch  die  SLAY-Karte,  eine  Ge- 
samtstruktur  aufgebaut  wird.  Jede  KOFL-Karte  mit  NREF  *  0  dcp- 
pelt  eine  in  den  vorangegangenen  KOFL-  oder  El ementkarten 
beschriebene  Struktur  iNREFi  mal.  Dabei  werden  in  Strukturen 


hoherer  Ordnung  automatisch  alle  Substrukturen  bis  hin  zum  Ein- 
zelelement  mi tgedoppelt .  Dadurch  konnen  auch  in  den  gedoppelten 
Strukturen  die  internen  Variablen  vollstandig  bestimmt  werden. 
Dies  unterbleibt  fur  Elemente  mit  NREP  <0,  in  denen  dann  nur  die 
externen  Variablen  zur  Verfiigung  stehen.  Mit  diesen  Elementen 
lassen  sich  z.B.  Symmetry-  und  Antimetriebeziehungen  ohne  nen- 
nenswerten  Mehraufwand  simulieren. 

Enthalten  die  KOPL-Karten  keinen  Wiederholungsf aktor  (NREP  =  0), 
werden  mehrere  Elemente  gleicher  Tief enstaf felung  zu  einem  Ele¬ 
ment  der  nachsthoheren  Ordnung  zusammengef aBt .  Hierbei  kann  nur 
auf  jungf rauliche ,  noch  nirgends  eingebaute  Elemente  zuriickge- 
griffen  werden. 

Die  eigentliche  Verknupfung  der  Knotenvari ablen  erfolgt  in  den 
KNUM-Karten .  Jeder  KOPL-Karte  ist  eine  KNUM-Karte  zugeordnet, 
die  den  Einbau  des  in  der  KOPL-Karte  beschriebenen  Elementes  in 
die  nachsthohere  Struktur  regelt.  Dabei  werden  mehrfach  gedop- 
pelte  KOPL-Karten  mitaezahlt. 

In  Bild  27  ist  die  Verknupfung  einer  Teilstruktur  aus  Schichten- 
elementen  beschrieben,  die  den  KOPL-Karten  5-8  zugeordnet 
sind.  Da  alle  Schi chte lemente  die  gleiche  lokale  Knotennumerie- 
rung  haben,  ist  hier  nur  ein  Element  dargestellt.  Sowohl  im 
Schichtelement  als  auch  in  der  Teilstruktur  werden  die  umrande- 
ten  Einzelkncten  zu  einem  Hauptknoten  zusammengef  a!3t .  In  den 
Schichtelementen  wird  der  mittlere  Hauptknoten  eliminert.  Die 
verbleibenden  externen  Variablen  werden  in  die  Teilstruktur  ein- 
gebaut  und  entsprechend  der  neuen  Hauptknotennumeri erung  umge- 
ordnet.  Dabei  konnen  auch  einzelne  Variablen  eines  lokalen 
Hauptknotens  unterschiedl ichen  Hauptknoten  der  globalen  Struktur 
zugeordnet  werden.  Urn  unnotige  Eingabedaten  zu  vermeiden,  ist 
die  vollstandige  Verknupfung  aller  Unbekannten  eines  lokalen 
Hauptknotens  mit  denen  eines  globalen  Knotens  durch  Voreinstel- 
lung  geregelt.  Im  DO-LOOF-Bereich  der  KNUM-Karte  werden  nur  die 
Abweichunaen  von  dieser  Regel  angegeben. 


Mit  den  Angaben  der  KOFL-  und  KNUM-Karten  lassen  sich  die  Steu- 
erfelder  zur  Beschreibung  der  Gleichunqsstruktur  auffauen. 


Bild  28.  Die  Blockmatrix  enthalt  Angaben  iiber  die  einzelnen 
Gleichungsblocke  und  deren  Anf angsadressen  auf  dem  Hintergrund- 
speicher.  In  der  DOF-Matrix  wird  die  Gleichungsstruktur  inner- 
halb  eines  Blockes  beschrieben  und  die  COPY-Matrix  regelt  den 
Einbau  des  vorliegenden  Elementsystems  in  die  nachsthdhere  Teil- 
struktur.  Damit  auch  die  in  einer  Laststufe  zu  erwartenden  Rand- 
gieichungen  eingetragen  werden  kbnnen,  miissen  die  Steuerfelder 
hinreichend  groB  gewahlt  werden. 

Um  eine  sparsame  Belegung  des  Hintergrundspei chers  zu  erreichen, 
wird  der  Programmablauf  nach  der  Transformation  der  Elemente  und 
dem  Einbau  der  Randbedingungen  optimiert.  Gleiche  Elemente  wer¬ 
den  dabei  erkannt  und  nur  einmal  auf  dem  Hintergrundspeicher 
abgelegt . 


G.  Ausblick 

i'as  skizzierte  Losungsschema  fur  die  Berechnung  von  physikalisch 
nichtl inearen  Froblemen  wurde  bereits  programmiert .  Bisher  ist 
al lerdings  erst  die  Incore-Vet  si  on  fur  die  Berechnung  linear 
olastischer  Systeme  f ert 1 ggost e 1 1 t .  Ein  f or tschrei tender  RiB 
wird  zur  Zeit  noch  durch  das  Losers  von  Kcppl ungen  simuliert. 
Einige  der  mit  diesem  Progrnmm  bcrechneten  Beispiele  sind  in 
[13  1  entha 1  ten . 

Hum  SchluB  mochte  ich  mich  mit  einer  Bitte  um  Hilfe  an  Sie  wen- 
den.  Es  gibt  moglicherweise  noch  eine  Methode,  um  die  Spannungs- 
singularitat  an  der  RiBspitze  zu  bestimmen.  Die  Grundlage  dieser 
Idee  soil  am  Beispiel  einer  Ei nf luB f 1  ache  fur  die  Spannungen 
skizziert  werden. 

In  einer  Scheibe  liefert  z.B.  das  Verschi ebunasf eld  fiir  v^  die 

33 

F:  l  nf  luB  f  lache  fiir  die  Spannuna  o  ,  wenn  im  Aufpunkt  ein  infi- 
nrtesimaler  Schlitz  in  die  Schoibe  gesc'nnitten  wird  und  die 
Schnittufer  um  die  Relati vver schi ebunq  Av  = 

B i Id  29. 


■v  %  ,-w  k>>' 


.  *  •  A  *  -  i  Jl *  i 


1  gespreizt  werden, 


Ersetzt  man  bei  einer  finiten  Ubersetzung  den  inf initesimalen 
Schlitz  durch  einen  endlicher  Lange  und  bringt  die  erforderliche 
Spreizung  an,  erhalt  man  nur  eine  grobe  Naherung  fur  die  Ein- 
fluBflache.  Der  auftretende  Fehler  liegt  bei  30  -  50%.  Legt  man 
dagegen  den  Aufpunkt  in  das  Innere  eines  finiten  Elementes  und 
ersetzt  die  Spreizung  durch  eine  aquivalente  eingepragte  Verzer- 
rung  ^^3'  liefert  die  finite  Ubersetzung  eine  iiberraschend  ge- 
naue  Naherung  fur  die  Einf luBf lache .  Selbst  in  den  Randknoten 
des  Elementes,  in  dem  der  Aufpunkt  liegt,  betragt  der  Fehler  nur 
wenige  Prozent. 
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In  Bild  29  wurden  nach  dieser  Methode  die  Einf luBf lache  [a  J 
fur  den  Mittelpunkt  einer  Quadratscheibe  errechnet  und  fur  ver- 
schiedene  Randbelastungen  ausgewertet.  Die  Einf luBf lache  und  die 
normale  finite  Elementberechnung  liefern  praktisch  identische 
Ergebni sse . 

Es  ist  iiberraschend,  wie  genau  hierbei  die  Singularitat  der  Ein- 
f luBf lache  im  Aufpunkt  approximert  wird.  Bei  der  Einf luBf lache 
geht  die  Singularitat  mit  1 /r  ~  und  nicht  wie  an  der  RiBspitze 
mit  1/Vr  -» 

Meine  Frage  lautet  nun: 

Gibt  es  eine  eingepragte  Verzerrung,  die  als  Freiwerte  vor- 
aussichtlich  noch  die  Spannungsintensitatsf aktoren  enthalt, 
mit  der  die  Singularitat  an  der  RiBspitze  approximiert  wer- 
den  kann? 


Wenn  das  zutrifft,  lieBe  sich  die  Berechnung  gerissener  Struktu- 
ren  erheblich  vereinf achen,  da  keine  Sonderelemente  fiir  die  RiB¬ 
spitze  mehr  gebraucht  wurden. 
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Loaded  cracked  structure 


4  Typische  Spannungsverteilung  in 
einer  gerissenen  Struktur. 


Extrapolationsmethode  zur  Bestimmung 
der  Spannungsintensitatsfaktoren. 


Bild  6  Bezeichnungen  zum  J- Integral. 


Bild  8  Energief rei setzungsrate  fur  einen 
virtuellen  RiBfortschritt . 


Energy  release  rate 


Ut  |U«  J-i  J 


Standard 

elements 


.Singular  elements 


Crack  tip  element  size  /  crack  size 


Bild  14  Nach  verschiedenen  Methoden  zu  erwartende  Fehler 

bei  der  Bestimmung  der  Spannungsintensi tatsf aktoren 
(Quelle:  Programm  PAM-FISS,  Fa.  ESI,  Paris). 
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Functional  Balance  forces 

6  uT  (A  (uo)-u  -  p  )  =  0  6j  =  A  (G-,  Hi 


Approximation  of  the  solution  vector 


1.  Step  u  =  u  i 

2. Step  ii  =  AfU!  ♦  u2 

3. Step  ii  =  X1*u  -j  +  A2’U2+u3 


bj  =  A  (u i  ) - u i 

Uj  =  Uj  /llUjll 


Solution  scheme 


Bild  16  Zwangsbedingungen  bei  groBen  Verformungen 


Constrained  algebraic  equations 
Az*g-E  =  r 
gT*z  +Ite=  s 


m  equations 
n  constraints 


Triangulansation  of  A 
A  =  LT  •  L 

Blockwise  solution 

z  =  I'M  L1'1-  r  -  LT* 1  -  g - e  ) 

( h  -  gT '  l_1  •  LT_1  g  )  •  e  =  s  -  gT  L _1  -  L1'1  r 

' - - - ^ - A - v - 

Solution  scheme  ^9  ^ 

1  Forward  reduction 

LT •  f  yr .  yg  ]  =  (  r,  g  ] 

(m:m)  ( m  :  n ♦  1 )  ( m  :  n+1  ) 

2. Equations  for  the  constraints 

(h-yj*yg)-e  *  s  -  yg-  yr 

(  n  :  n  )  I  n :  1 )  ( n :  1 ) 

3  Backwards  reduction 

L  •  z  =  yr  -  yg-E 

(m:m)  (m:1)  (m:1) 


Bild  18  Einschr ankung  von  algebraischen  Gleichungen 
iiber  zusatzliche  Lastvektoren. 
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x  Resolved  part 

S  i 


Pivot  characteristic 

Before  decomposition  After  decomposition 
5;  =  1  :  in  aeneral  S;  =  ±8  :  Pivot  £  0 


in  general 

Si  =  ±8  : 

Pivot 

§  0 

Pivot  =  0 

7  : 

Pivot 

■  0 

BC  =*=0 

6  : 

BC 

+  0 

BC  =0 

5  : 

BC 

-  0 

BC  =  Boundary  conditions  for  I<  R  only 


Bild  20  Bezeichnungen  fur  die  Cholesky-Routine. 


...a 


I  =  M/-1/1 


gw 
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Influence  surface  ( ov3 3 ]  in  the  field  point  F  for  vertical  load 


Bild  29  E; nf tuBf lache  [o^  ]  fur  den  Mittelpunkt  F  einer 
Qwadratschei.be  unter  vertikaler,  symmetr  i  sober 
Be  I  as  hung. 


Solution  scheme 
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Fig  2,  the  stress-strain  diagram  is  non-linear.  In  order  to  sup¬ 
port  the  actual  test  conditions  the  elasticity  constants,  listed 
in  Tab.  2,  were  proportionally  reduced  until  the  calculated  dis¬ 
placements  correspond  to  the  measured  displacements  at  the  elas¬ 
tic  limit  load.  The  such  modified  constants  were  also  used  for 
higher  load  levels. 

Fig.  3  shows  a  comparison  of  measured  and  calculated  force-dis¬ 
placement  curves  for  the  test  specimen  elongation  wl  and  for  the 
average  of  the  notch  opening  displacement  COD.  Discrepancies 
between  measured  and  calculated  values  occur  only  in  the 
wl -curves  when  the  crack  tips  approach  the  terminals  of  the  dis¬ 
placement  transducers.  The  terminals  are  not  attached  to  discrete 
points  at  the  specimen  center  line  but  to  clamped-on  yokes 
extending  over  the  entire  specimen  width.  Since  in  the  presence 
of  long  cracks  the  displacements  over  the  width  of  the  specimen 
are  not  constant  (see  Fig.  4)  and,  because  of  the  yoke  defor¬ 
mations,  the  clamping  stresses  at  the  edges  of  the  specimen  are 
higher  than  at  the  center,  a  more  pronounced  curvature  in  the  mea¬ 
sured  wl -curve  should  be  expected  than  in  the  calculated  curve 
related  to  the  specimen  center  line.  The  measurements  are  not 
indicative  of  a  possible  plastic  deformation  at  the  crack  tip  as 
is  evident  from  the  unloading  curve  in  Fig.  5,  which  terminates  at 
the  point  of  origin. 

The  majority  of  the  experiments  led  to  a  crack  pattern  which  is 
almost  point-symmetric  relative  to  the  center  of  the  specimens. 
The  consistency  of  this  phenomenon  could  be  confirmed  by  a  finite 
element  calculation  on  a  notched  (°° Je"1  aminate  with  a  single 
10  mm  long  crack  at  one  of  the  notch  tips.  The  calculated  stresses 
along  the  lines  of  potential  crack  propagation  are  given  in 
Fig.  6.  The  fracture  process  is  influenced  predominantly  by  tli^ 
normal  stresses  transverse  to  the  crack  and  hardly  by  the  shear- 
stresses,  as  verified  in  report  [2].  Therefore  only  the  normal 
stresses  must  be  compared,  which  are  quite  similar  around  differ¬ 
ent  crack  tips.  The  aberrations  from  the  point-symmetric  fracture 
mode,  observed  in  ca.  30%  of  the  test  specimens,  may  be  due  to 
scatter  of  the  material  strength  properties.  For  the  analysis  of 
the  measurement  a  point-symmetric  crack  pattern  is  assumed. 


In  UD-laminates  the  cracks  start  always  at  the  notch  tip.  Two 
cracks  occur  at  each  notch  tip  in  on-axis  laminates  and  one  crack, 
declined  to  the  specimen  axis,  in  off-axis  laminates.  Therefore 
four  or  two  force-crack  length  curves,  respectively,  must  be  ana¬ 
lyzed  simultaneously.  According  to  the  location  of  a  crack  in  the 
specimen  the  curve  is  plotted  in  the  correspondent  quadrant  of  a 
coordinate  system  as  illustrated  in  Fig.  7.  The  force-crack 
length  curves,  obtained  for  the  individual  tests  in  a  test 
series,  are  plotted  on  transparent  paper,  demonstrated  in  Fig.  8 
to  13  for  an  [ 0° ] Q-laminate .  Then  the  transparencies  are  put  on  a 
light  table  and  turned  over  and  over  again  at  one  or  both  of  the 
symmetric  axis  of  the  specimen,  until  the  scatter  for  all  of  the 
curves  is  as  small  as  possible.  For  the  | 0°  ] g- laminate  the  super¬ 
position  of  all  of  the  force-crack  length  curves  is  shown  in 
Fig.  14. 

Because  of  the  assumed  point- symmetric  crack  pattern,  the  crack 
lengths  in  diagonal  opposite  quadrants  of  the  specimens  are  aver¬ 
aged  in  each  load  step.  All  of  the  corresponding  load-crack 
length  curves  are  plotted  into  one  diagramm  and  approximated  by  a 
fitting  curve.  Analogically  the  notch  opening  displacements  COD1 
and  COD2  are  averaged  and  plotted  versus  the  load,  whereas  into 
the  load-elongation  diagramm  the  test  values  for  wl  were  intro¬ 
duced  directly. 

For  off-axis  laminates  the  stable  crack  lengths  become  shorter 
and  shorter,  when  the  fiber  angle  increases.  Occasionally  only 
one  or  two  measured  values  are  available.  For  off-axis  laminates 
the  test  value  are  listed  point-by-point  in  the  diagrams  without 
a  connecting  line. 

In  view  of  the  calculation,  described  in  the  third  part  (2J,  the 
fitting  curves  are  picked  up  pointwise  as  functions  of  the  crack 
lengths  and  assembled  in  tables  for  the  different  fiber  angles 
and  specimen  geometries. 


3.  ANALYSIS  OF  THE  TESTS 


3.1  [0°]n-Laminates 

For  simple  precalculations  a  symmetric  crack  pattern  is  assumed, 
such  that  only  a  quarter  of  the  specimen  must  be  analysed  numer¬ 
ically.  In  this  case  the  four  different  crack  length  within  a  spe¬ 
cimen  are  averaged  in  each  load  step.  For  the  mean  crack  lengths 
1  the  average  stresses  at  the  grip-ends  are  plotted  in  Fig.  15. 
The  chosen  fitting  curve  is  marked  by  a  fat  solid  line. 

As  pointed  out  in  the  first  part  [1],  misaligned  fibers  or  layers 
act  as  a  crack  stopper.  This  might  be  the  reason,  that  in  some 
specimen  the  point-symmetric  crack  pattern  is  blocked  and  a  more 
symmetric  one  appears.  In  the  figures  of  this  chapter  a  symmetric 
crack  pattern  is  marked  in  the  plots  by  a  dashed  line.  They  indi¬ 
cate,  that  for  a  symmetric  crack  pattern  the  critical  stresses, 
where  the  cracks  start  to  propagate,  are  often  much  higher  than 
those  for  the  point-symmetric  mode. 

In  Figs.  16  and  17  the  crack  lengths  lj  and  I2  are  plotted  for  a 
point-symmetric  crack  pattern  and  different  mean  stresses  at  the 
grip-end.  The  mean  values  of  both  curves  lead  to  the  fitting  curve 
in  Fig.  15.  Fig.  18  shows  the  elongation  wl  of  the  specimen  and 
Fig.  19  the  notch  opening  displacements  COD  as  function  of  the 
mean  stresses  respectively. 

In  order  to  justify  the  assumption  of  the  point-symmetric  crack 
pattern  in  Figs.  20  and  21  the  notch  opening  displacements  are 
plotted  for  the  individual  notches  and  Fig.  22  shows  the  superpo¬ 
sition  of  both  figures.  In  all  three  plots  the  idealization, 
determined  in  Fig.  19  for  a  point- symmetric  crack  pattern  is 
introduced. 

In  Figs.  15  to  19  the  chosen  fitting  curves  for  a  [0° ]g-laminate 
are  drawn  and  the  corresponding  numerical  values  are  assembled  in 
Tables  3  and  4 . 


In  Figs.  23  to  27  the  crack  lengths,  specimen  elongations  and  the 
notch  opening  displacements  for  a  [  0° ] ^-laminate  are  plotted  as 
functions  of  the  mean  stresses. 

3.2  [5°]g-Laminates 

In  off-axis  laminates  cracks  start  always  at  the  notch  tips  and 
arise  only  in  diagonal  opposite  quadrants  of  the  specimen.  The 
width  of  the  specimen  tested  are  50  mm  and  20  mm  respectively. 

For  the  50  mm  wide  specimen  the  crack  lengths,  the  specimen  elon¬ 
gations  and  the  notch  opening  displacements  are  plotted  versus 
the  load  in  Figs.  28  to  31.  The  values,  picked  up  for  the  fitting 
curves,  are  collected  in  Table  5.  Similarly,  for  the  20mm  wide 
specimen  the  test  results  and  the  fitting  curves  are  given  in 
Figs.  32  to  35  and  the  corresponding  numerical  values  are  listed 
in  Table  6 . 

In  the  50  mm  wide  specimen  unbroken  fibers  run  from  one  clamped 
end  to  the  other.  When  the  cracks  propagate  the  net  cross  section 
is  reduced  and  predominantly  the  unbroken  fibers  take  over  the 
stress  increase.  Because  of  the  growing  crack  length  the  strains 
at  the  crack  tips  are  released  and  the  propagation  stops.  There¬ 
fore  the  load-crack  lengths  curve  (Fig.  28)  shows  a  significant 
increase  of  the  load  for  a  growing  crack. 

In  the  20  mm  wide  specimen  the  relation  is  quite  different.  All  of 
the  fibers  at  least  at  one  side  end  at  a  free  surface.  Therefore, 
in  front  of  the  crack,  the  load  must  be  carried  completely  by  the 
resin.  Theoretically  the  cracks  will  propagate  unstably  right 
from  the  beginning,  because  the  stresses  at  the  crack  tip 
increase  with  the  crack  length.  The  small  load  increase  between 
the  crack  initiation  and  the  ultimate  state,  shown  in  Fig.  32, 
might  be  due  to  prior  damage  introduced  into  the  specimen  by  the 
notch  fabrication. 
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3.3  [ 10° ]g- Laminates 

For  specimen  with  50  mm  width  the  crack  lengths,  the  specimen 
elongations  and  the  notch  opening  displacements  are  given  as 
functions  of  the  load  in  Figs.  36  to  39.  The  corresponding  numer¬ 
ical  values  for  the  fitting  curves  are  listed  in  Table  7  for  dif¬ 
ferent  crack  lengths.  For  the  [ 10° ] g- laminates  the  crack 
initiation  could  not  be  observed.  After  a  short  stable  crack  pro¬ 
pagation  the  crack  moved  off  unstably  and  the  specimen  separates. 


3.4  [ 15° ]g- Laminates 

For  specimen  with  50  mm  width  the  crack  lengths,  the  specimen 
elongations  and  the  notch  opening  displacements  are  plotted  ver¬ 
sus  the  load  in  Fig.  40  to  43.  The  corresponding  numerical  values 
of  the  fitting  curves  are  listed  in  Table  8.  The  largest  fiber 
angle,  where  a  stable  crack  propagation  could  be  observed  in  the 
wide  specimen,  is  =  15°,  whereas  the  small  specimen  failed 
suddenly. 


3.5  [30°]g-Laminates 

In  tests  on  UD-laminates  with  fiber  angles  of  aF  £  30°  the  failure 
occured  suddenly  without  any  indication  of  a  crack  initiation  and 
the  meassured  load-elongation  curves  were  straight  lines.  In 
Table  9  the  averaged  test  results  were  given  for  side  notched 
( 30° ] g-laminates  with  specimen  widths  of  20  mm  and  50  mm,  respec¬ 
tively.  In  addition  the  limit  values  for  unnotched  specimen  are 
I  listed  for  comparison. 

3.6  [45°]n-Laminates 

i 

In  Table  10  the  test  results  are  given  for  notched  and  unnotched 
!  [ 45 ° ] q-  and  [  45° 1 16-laminates  . 

k 

! 
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3.7  [90° ]n- Laminates 

For  [90° ]q-  and  [90° ) 16-laminates  with  different  notch  forms  the 
test  values,  measured  under  ultimate  static  load,  are  summarized 
in  Table  11.  Compared  with  the  sharp  notch  tip  according  to 
Fig.  1  a  more  blunt  saw  cut  will  give  a  significantly  higher  ulti¬ 
mate  load.  Contrary  to  metals  the  laminate  is  quite  insensitive 
to  surface  defects  as  the  test  with  a  one-sided  blunt  saw  tap 
indicates . 


4.  RESULTS 


According  to  Table  1  only  the  | aP ] Q- laminates  with  sharp  notches 
and  50  mm  specimen  width  are  tested  thoroughly  for  all  fiber 
angles.  Smaller  specimens,  stacking  sequences  with  16  layers  or 
blunt  notch  forms  are  tested  only  occasionally.  These  tests  are 
used  to  verify,  that  in  the  [ ap ] g- laminates  with  sharp  notches  the 
first  cracks  initiate  at  mean  stresses  in  the  net  cross  section, 
which  are  lower  than  those  in  specimen  of  other  configurations. 
This  phenomenon  is  illustrated  in  Figs.  44  to  46  for  on-axis  lami¬ 
nates  of  different  thicknesses.  In  this  case  the  misalignment  of 
fibers  and/or  layers  may  act  as  a  crack  stopper,  which  happens  in 
thick  laminates  more  often  than  in  thin  one. 

In  Fig.  47  the  load-crack  lengths  curves  are  plotted  for  sharp 
notched  { ap ] g-laminates  with  different  fiber  angles.  The  length 
of  the  stable  crack  becomes  shorter  and  shorter,  when  the  fiber 
angle  referred  to  the  load  decreases.  For  ap  =  15°  the  crack  ini¬ 
tiation  can  just  be  observed  before  the  specimen  fails  cat¬ 
astrophically. 

The  measured  mean  stresses  for  the  crack  initiation  are  drawn  in 
Fig.  48  versus  the  fiber  angles  ap.  In  the  plot  the  ap-axis  is 
extended  in  a  logarithmic  scale  (ap  [°]  =  100xlog(ap°)  [mm)). 
For  ap  =  5°  to  45°  the  measured  values  are  about  in  alignment.  Sim¬ 
ilar  to  a  Wohler's  curve  the  straight  line  is  extended  at  both  end 
by  lines  of  constant  level  to  fit  the  test  values  for  ap  =  0°  and 


90°.  In  order  to  support  the  chosen  approximation  additional 
tests  are  desirable  for  <ip  =  2°,  3°  and  60°.  According  to  the 
approximation  curve  the  mean  stresses  for  the  crack  initiation 
will  stay  constant  for  fiber  angles  of  ap  <  2°  and  otp  >  52° 
respectively . 

For  unidirectional  laminates  the  results  received  from  the  tests 
directly  can  be  summarized  by  the  following  statements: 

•  The  cracks  develop  always  parallel  to  the  fibers. 

•  The  crack  propagation  is  significantly  delayed  by  inevitable 
misalignment  of  fibers,  fiber  bundles  or  layers. 

•  In  thin  laminates  the  crack  initiation  starts  earlier  but 
propagate  with  less  speed  than  in  thick  ones. 

The  second  point  together  with  Fig.  48  indicate,  that  in  a  packet 
of  equal  orientated  UD-layers  an  alternating  misalignment  of 
about  Actp  =  ±2°  in  the  single  layers  will  reduce  the  crack  propa¬ 
gation  significantly  without  any  loss  in  the  laminate  strength. 
This  statement  is  quite  important  for  the  practical  design  of  new 
laminates  and  should  be  verified  by  further  tests. 
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( 0° ) q- laminate  along  the  expected  crack  directions 
calculated  by  finite  elements. 
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Figure  7  Quadra 
with  t 


Crack  length  in  mm 


Figure  8 


Measured  crack  length  under  static  load 
Test  specimen  no.  781/1 
Specimen  width  50  mm 

Material  T300/914C 


Stacking 


sequence  [0°]g 


Crack  length  in  mm 


Figure  9  Measured  crack  length  under  static  load 
Test  specimen  no.  781/2 
Specimen  width  50  mm 
Material  T300/914C 

Stacking  sequence  [0°Jo 


Tensile  force  in 
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Figure  12  Measured  crack  length  under  static  load. 
Test  specimen  no.  781/6 
Specimen  width  50  mm 
Material  T300/914C 

Stacking  sequence  [0°]g 


Average  stress  Fc/A  at  the  grip-end  of  the  specimen 


20  Measured  notch  opening  displacements  C0D1 
Test  specimen  no.  781/1,2,4-7 
Specimen  width  50  mm 

Material  T300-914C 

Stacking  sequence  [0°]g 


(TMT|2  I  Antisymmetric  v  Y/Z 


crack  pattern 

_ J _ 


■Symmetric 
crack  pattern 


Idealization 


Crack  initiation 


Elastic  limit 


v=  const  / 


Notch 


Hi— 

3 ♦ C0D2 


0  0.05  0.10  0.15  0.20  mm 

Notch  opening  displacements  C0D2  - 


Figure  21  Measured  notch  opening  displacements  COD2 
Test  specimen  no.  781/1,2,4-7 
Specimen  width  50  mm 
Material  T300-914C 

Stacking  sequence  [G°]g 
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Figure  23  Measured  average  crack  length  1 
for  symmteric  crack  pattern. 
Test  specimen  no.  778/2-7 
Specimen  width  50  mm 

Material  T300-914C 

Stacking  sequence  (0°J1£ 


Average  stress  Fc/A  at  the  grip-end  of  the  specimen 


Idealization 
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Figure  24  Measured  average  crack  length  lj 
for  antisymmetric  crack  pattern. 
Test  specimen  no.  778/2-7 
Specimen  width  50  mm 
Material  T300-914C 

Stacking  sequence  [ 0 ° ^ ^ 


•*.  1 
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Figure  25  Measured  average  crack  length  I2 
for  antisymmetric  crack  pattern. 
Test  specimen  no.  778/2-7 
Specimen  width  50  mm 
Material  T300-914C 

Stacking  sequence 


Idealization 


Crack  initiation 

V 


Elastic  limit 


v  =const 


0  0.05  0.10  0.15  0.20  mm 
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Figure  27  Measured  average  notch  opening  displacements 
Test  specimen  no.  778/2-7 
Specimen  width  50  mm 

Material  T300-914C 

Stacking  sequence  [0°]1C 


m  '  « 


Average  stress  Fc/A  in  N/mm2 


[Mia 
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_ J  .02  .04  .06  ,08  .11 

Notch  opening  displacements  C0D1 ,  C0D2  in  mm 


Figure  31  Measured  notch  opening  displacements  for 
both  notches. 

Test  specimen  no.  891/1-4 
Specimen  width  50  mm 
Material  T300-914C 

Stacking  sequence  [ 5  0  J Q 


Average  stress  Fc/A  in  N/mm2 


in  N/m 


Average  notch  opening  displacements  COD  in  mm 


gure  34 


Measured  average 
Test  specimen  no. 
Specimen  width 
Material 

Stacking  sequence 


notch  opening  displacement 
891/1-4 
20  mm 
T300-914C 


Average  stress  FcrA  in  N/mm2 


Average  notch  opening  displacements  COD  in  mm 


Figure  38  Measured  average  notch  opening  displacements 
Test  specimen  no.  890/1-4 
Specimen  width  50  mm 

Material  T300-914C 


Stacking  sequence  f  1C0],-. 


Notch  opening  displacements  C0D1>  C0D2  in  mm 

Figure  39  Measured  notch  opening  displacements  for 
both  notches. 

Test  specimen  no.  890/1-4 
Specimen  width  50  mm 

Material  T300-914C 

Stacking  sequence  |10°]n 


Average  stress  Fc/A  in  N/mm2 
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Figure  44  Comparision  of  crack  lenths  for 
(  0°  ]8-  and  1  0°  ]  laminates . 
Material  T300-914C 
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7.  TABLES 


Stacking 

sequence 

Specimen  width  in  mm 

50  I  20 

No.  of  tests 

[0°], 

[0°],6 

6 

6 

I5°]8 

4 

1 

3 

1 

H0°]B 

4 

1 

1 

[15°]8 

4 

1 

4 

| 

130o18 

4 

1 

4 

MM 

■EH 

5 

5 

2 

1 

(90°18 
190°), 6 

2 

3 

3 

1 

3 

Notch  form: 
A 

°7  _l 

B 

1 

C 

J  J  I 

u 

Unnotched 

x: _ 

. 

10  40 

to 

o 

1  50  r 

n  OJL 

in 

cn 

1  c 

d 

o 

Table  1  Test  program  for  side  notched  specimen 
under  quasi-static  tensile  load. 
Material  T300-914C 


«*, 

n 


Fiber  direction 


Cartesian  coordinates 


1 

El 

Sym. 

-  y?1 

IE L 

1 

E  2 

0 

0 

1 

G 

*22  ♦  *22  'T 

_  12  n 


Constitutive  equations 


T 

23 

Temperature  in  °C 

M 

0 

50 

Moisture  in  % 

1784 

1889 

Ultimate  stresses 

IT11 

°D 

-  1086 

-  1  090  ■ 

in  N/mmz 

'  <=* 

^  r-ol 

r-'Oi 

64.8 

55.6 

(  )z  :  Tension 

nll 

°D 

-  215 

-  212 

(  )0  :  Compression 

«U 

al2 

92.4 

94.4 

Ei 

143  242 

150  782 

Elasticity  constants 

V21 

0.317 

0.325 

in  N/mm2,  measured 

E; 

10  165 

9  987 

by  30%  of  the 

G 

6  239 

5  859 

ultimate  stresses 

To 

190° 

Glas  transition  point  in 

«11 

-0.8- 

10’6 

"C.Temperatur  expan  = 

« *22 

°q 

GO 

Ovl 

10'6 

sion  coefficient  in  eC_1 

Table  2  Material  constants  for  UD-laminates  of  T300/914C 


v  =  const 


3  ♦COD 


1.25 


2,50 


3.75 


5.0 


7.5 


10,0 


12.5 


15.0 


17.5 


20.0 


Fc/A 

wl 

COD 

N/mm2 

mm 

mm 

110,0 

0,120 

0,057 

122,0 

0,134 

132,0 

0,145 

0,076 

138,0 

0,153 

0,084 

U2,5 

0,161 

0,092 

148,0 

0,171 

0,104 

152,5 

0,180 

0,116 

156,0 

0,188 

0,126 

159,0 

0,195 

0,136 

161,5 

0,201 

0,145 

164,0 

0,207 

mm 

30,0 

171,0 

0,226 

0,190 

40,0 

176,0 

0,239 

0,218 

n  test  values  for 
tic  load  with  sym 
t  specimen  no.  78 
cimen  width  50 
erial  T3 


it 


>> 


^\\\\\\\ 


3  +  C0D2 
125  +  wl  - 


Specimen 

no . 

Width 

mm 

*c 

mm 

Fc/A 

N/mm2 

wl 

mm 

COD 

mm 

891/1-4 

50 

0. 

87.6 

0.103 

0.041 

1.25 

100.4 

0.118 

0.048 

2.5 

103.9 

0.  123 

0.050 

3.75 

106.2 

0.127 

0.053 

5. 

107.5 

0.129 

0.055 

7.5 

109.3 

0.131 

0.058 

10. 

110.6 

0.136 

0.061 

15. 

112.4 

0.  141 

0.065 

20. 

113.9 

0.144 

0.069 

25. 

114.7 

0.  147 

0.071 

30. 

115.5 

0.149 

0.073 

40. 

116.1 

0.153 

0.075 

50. 

116.7 

0.155 

0.077 

60. 

117.3 

0.157 

0.079 

70. 

117.9 

0. 160 

0.081 

Notch:  Ultrasonic  knife  cut, 

5mm  deep,  0.6  mm  wide,  0 . 1  mm  tip  radius 
Young's  Modulus  in  load  direction:  E  =  111000  N/mm2 


Table  5  Mean  test  values  for  [ 5 ° J g- 1 aminates  under 

static  load  with  antisymmetric  crack  pattern. 
Test  specimen  no.  891/1-4 
Specimen  width  50  mm 

Material  T300/914C 


j. 


J  Specimen 

Width 

1 

c 

Fc/A 

wl 

COD 

j  no . 

mm 

mm 

N/mm* 

mm 

mm 

891/1-4 

20 

0. 

O 

CO 

>j0 

0.089 

0.0372 

1 . 25 

70.1 

0.094 

0.0389 

2 . 5 

71.6 

0.096 

0.0402 

3 . 75 

72.3 

0.098 

0.0407 

5. 

72.7 

0.099 

0.0411 

7.5 

73.4 

0.  10  0 

0.0417 

10. 

73.9 

0.101 

0.0421 

15. 

74.3 

0. 102 

0.0426 

>=20. 

74.6 

0.103 

0.0430 

Notch:  Ultrasonic  knife  cut, 

5mm  deep,  0.6  mm  wide,  0.1  mm  tip  radius 
Young's  Modulus  in  load  direction:  E  =  108000  N/mm1 


Table  6  Mean  test  values  for  [ 5 ° ] Q- 1 amina tes  under 

static  load  with  antisymmetric  crack  pattern. 
Test  specimen  no.  891/1-4 
Specimen  width  20  mm 
Material  T300/914C 


xWWWW* 


3  +  C0D1 

v  =  const. 

-1  t~“ 

-  ■  .  1  [  JO 

=  15' 


°  V/  s 


Fc/A 

wl 

N/mm2 

mm 

0. 

j  51.8 

0. 125 

j  53.6 

0.25 

i  54.2 

0.37  5 

54.6 

0.5 

55.0 

0.75 

55 . 6 

1.0 

56.1 

1.25 

56.5 

<D 

in 

e 

<tf 

e 

0) 

\ 

u 

u 

2: 

e 

o 

■H 

o 

U 

o 

nj 

un 

0.0309 

0.0317 

0.0321 

0.0326 

0.0328 

0.0332 

0.0338 

0.0342 

0.0346 

0.0353 

0.0357 

0.0360 

0.0362 

0.0364 

0.0364 


Noth:  Ultrasonic  knife:  cut, 

5mm  deep,  0.6  mm  wide,  0.1  mm  tip  radius 
COD  =  ( C0D1  +  COD2)/2 


Table  8  Moan  test  values  for  [  15° Jg-laminates  under 

static  load  with  antisymmetric  crack  pattern. 
Test  specimen  no.  888/1-4 
Specimen  width  50  mm 

Material  T300/914C 


; 


const . 


3  +  C0D1 


( COD1  +  COD2)/2 


3  +  C0D2 
125  +  wl 


“SHE” 


Specimen 

no . 

Width 

mm 

Fc/A 

N/mm2 

wl 

mm 

889/1-4 

20. 

22.9 

0.131 

50. 

28.3 

0.167 

20. 

25.0 

0.140 

125.0 

0.998 

167.6 

— 

776/6-7 

50. 

40.0 

0.224 

140.0 

0.998 

150.6 

-  - 

Notch 

form 


0.0238 

0.0283 


Notch  A:  Ultrasonic  knife  cut, 

5  mm  deep,  0.6  mm  wide,  0.1  mm  tip  radius 
U:  Unnotched 


Table  9  Mean  test  values  for  [ 30° ] g- laminates 
under  ultimate  static  load. 

Material  T300/914C 


ife  ,  . 


jpecimen 

no . 

Stacking 

sequence 

Fc/A 

N/mm1 

wl 

nun 

COD 

nun 

776/1-5 

[45"], 

17.4 

0.188 

0.0250 

777/1-5 

[ 45" ] 1, 

19.3 

0.208 

0.0289 

776/6-7 

[45"], 

92.6 

1.200 

777/6-7 

[45*  ]x« 

91.8 

1.214 

Notch  A:  Ultrasonic  knife  cut, 

5  nun  deep,  0.6  nun  wide,  0.1  nun  tip  radius 
U:  Unnotched 


Table  10  Mean  test  values  for  [ 45° }  -laminates 

1  '  n 

under  ultimate  static  load. 

Material  T300/914C 


71 


COD  =  (COD1  +  COD2)/2 


Specimen 

no . 

Stacking 

sequence 

Fc/A 

N/mm2 

wl 

mm 

COD 

mm 

— 

785/1-5 

[ 90“ J  * 

12.1 

0.198 

0.0245 

784/1-5 

[90a]lt 

13.9 

0.223 

0.0237 

785/4,5,7 

190"]. 

17.9 

0.262 

0.0350 

785/6 

(90"], 

56.0 

0.840 

784/4-7 

[  90“  J  i* 

56.3  | 

0.844 

Notch 

form 


Notch  A:  Ultrasonic  knife  cut, 

5  mm  deep,  0.6  mm  wide,  0.1  mm  tip  radius 
B:  Blunt  saw  cut,  5  mm  deep,  0.35  mm  wide 
C:  One-sided  blunt  saw  cut,  0.1  mm  deep,  0.35  mm  wide 
U:  Unnotched 


Table  11  Mean  test  values  for  (90°)  -laminates 

'  n 

under  ultimate  static  load. 

Material  T300/914C 
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^  -  %  * • * , 
'  •*  ■  *.  -  .  '  J  -  . 

v  V 

VV//1' 

vV'.*V 

V  \  *w' 


APPENDIX  E 


EinfluB  von  Temperatur  und  Feuchte  auf 
die  Schwingfestigkeit  ungekerbter  CFK  - 
Laminate 

von 

H.  C.  Goetting  und  W.  Ganz 


DFVLR  -  Institut  fur  Strukturmechanik 


DEUTSCHE  FORSCHUNGS-  UND  VERSUCHSANSTALT 
FUR  LUFT-  UND  RAUMFAHRT  E.V. 
INSTITUT  FUR  STRUKTURMECHANIK 


Interner  Bericht 
IB  131-84/52 


EinfluB  von  Temperatur  und  Feuchte  auf 
die  Schwingfestigkeit  ungekerbter  CFK  - 
Laminate 


Braunschweig,  Dezember  1984 


(Dr.-Ing.  H. Wl Bergmann) 


Dieser  Bericht  umfaBt: 

60  Seiten  mit 
48  Bildern 


Die  Bearbeiter: 

v.  r. 

(Dipl.-Phys.  H.C.Goetting) 


( Ing .  W.  Ganz ) 


1.  EINLEITUNG 


Im  Rahmen  des  ESA-Auftrags  "Development  of  Fracture  Mechanisms 
Maps  for  Composite  Materials"  wurde  zur  Untersuchung  des  Einflus- 
ses  von  Temperatur  und  Feuchtigkeit  auf  das  Schwingfestigkeits- 
verhalten  sowie  der  Art  und  Ausbreitung  verschiedener  Schadi- 
gungsmechanismen  unter  Schwingbelastung  ein  sehr  umfangreiches 
Testprogramm  an  ungekerbten  Flachprobestaben  mit  ausgewahlten 
Lagenaufbauten  durchgefiihrt . 

Folgende  Laminataufbauten  wurden  nach  MaBgabe  von  Table  A  des 
ESA-Proposals  ausgewahlt: 


Unidirektional , 

on-axis,  0° 

8 

Lagen 

Unidirektional , 

off-axis,  5° 

8 

Lagen 

Unidirektional, 

off-axis,  15° 

8 

Lagen 

Unidirektional , 

off-axis,  45° 

8 

Lagen 

Multidirektional , 

[02°/+45o/02o/-45o/0°/90o1s 

16 

Lagen 

Multidirektional , 

[ 0°/+45o/-45o/90° ] 

s 

7 

Lagen 

Angle-ply,  ±15 

o 

8 

Lagen 

Angle-ply,  ±30 

o 

8 

Lagen 

Angle-ply,  ±45 

0 

8 

Lagen 

Die  Proben  bestanden,  wie  alle  im  Rahmen  des  ESA-Auftrags  unter- 
suchten,  aus  T300/914  C  Prepregs,  die  in  einem  rechnergesteuerten 
Autoklaven  unter  stets  gleich  gehaltenen  Bedingungen  verarbeitet 
wurden. 
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2.  PROBENVORBEREITUNG  UND  VERSUCHSDURCHFuHRUNG 

Mit  zwei  verschiedenen  Temperaturen  und  zwei  verschiedenen 
Feuchtegraden  wurden  vier  unterschiedliche  Klimazustande  in  fol- 
gender  {Combination  simuliert: 


• 

Normaltemperatur 

-  trocken 

(ad) 

• 

Hohe  Temperatur  (80°) 

-  trocken 

(hd) 

• 

Normaltemperatur 

-  feucht  (95%  RH) 

(aw) 

• 

Hohe  Temperatur 

-  feucht 

(hw) 

j  Diese  Klimazustande  wurden  mit  selbstentwickelten  kleinen  Klima- 

kammern  erzeugt,  wobei  sich  nur  die  ProbenmeBlange  von  140  mm  in 
der  Kammer  befand,  urn  einen  EinfluB  auf  die  Tragf ahigkei t  der 
Probenarmierung  zu  vermeiden  und  die  zu  klimati sierende  Masse 
|  moglichst  klein  zu  halten. 

Fur  die  Untersuchungen  unter  95%  RH  wurden  die  Proben  in  einer 
speziellen  Klimakammer  mit  ihren  MeBlangen  bei  80°  C  und  95%  RH 
solange  ausgelagert,  bis  nach  vorliegenden  Referenzkurven  mit 
|  einem  Sicherhei tsf aktor  von  1.5  die  Sattigungsaufnahme  erreicht 

war.  Die  trockenen  Proben  wurden  vor  dem  Test  6  Wochen  in  der 
Prvifhalle  bei  konstanter  relativer  Luftfeuchte  von  40%  und  kon- 
stanter  Temperatur  von  22°  C  gelagert. 

■  Da  sich  bei  Normaltemperatur  und  95%  RH  durch  den  Wasserdampf  die 

;  Temperatur  in  der  Kammer  auf  28°C  erhohte,  ergaben  sich  folgende 

j  Eckwerte  fur  die  vier  Klimazustande,  die  wahrend  des  Versuchs 

>  konstant  gehalten  wurden: 

i 


ad: 


22  °C 


40%  RH 


hd: 

O 

0 

O 

00 

- 

5% 

RH 

aw: 

28°C 

95% 

RH 

hw: 

U 

0 

O 

03 

- 

95% 

RH 

Die  Schwingfestigkeitsversuche  wurden  durchweg  im  Zugschwellbe- 
reich  bei  einem  Verhaltnis  von  Unterlast  zu  Oberlast  von  R=0 . 1 
durchgef  iihrt .  Die  Belastungsvorgabe  erfolgte  uberwiegend  durch 
pro jrammierbare  mikroprozessorgesteuerte  Funktionsgeber  und 
durch  ProzeBrechner  (PDP  11),  wobei  eine  einheitliche  Bela- 
stungsf requenz  von  5  Hz  eingehalten  wurde.  Bei  ca.  70%  aller 
Schwingf estigkeitsversuohe  wurde  standig  die  Probendehnung  mit- 
gemessen,  urn  die  Steif igkeitsabnahme  bei  konstant  gehaltener 
Last  bis  zum  Bruch  verfolgen  zu  konnen. 

Die  Ermittlung  der  Lebensdauerlinien  (Wohlerkurven)  erfolgte  auf 
die  Weise,  daB  auf  drei  oder  vier  Spannungshorizonten  die  Bruch- 
lastspielzahlen  mit  jeweils  4  bis  6  Proben  gemessen  und  stati- 
stisch  ausgewertet  wurden,  wobei  die  Linien  fur  10%,  50%  und  90% 
Dberlebenswahrscheinlichkei t  in  Abhangigkeit  von  den  Logarithmen 
der  Lastspielzahlen  berechnet  wurden. 

Weitere  statistische  Auswertungen  wurden  mit  der  Ermittlung  der 
Summenhaufigkeitsverteilungen  im  Weibull-Netz  einschlieBlich 
der  Vertrauensbereiche  durchgefiihrt .  Statistische  Signifikanz- 
tests  wurden  immer  dann  angewendet,  wenn  durch  teilweise  sehr 
groBe  Streuungen  die  Erkennung  signif ikanter  Unterschiede  der 
Mittelwerte  eines  Spannungsniveaus  nicht  moglich  war. 

Zur  Charakterisierung  typischer  Bruchbilder  in  Abhangigkeit  von 
Laminataufbau  und  Umgebungsbedingungen  wurden  lichtmikroskopi- 
sche  Aufnahmen  und  rasterelektronische  Detailaufnahmen  von  aus- 
gewahlten  Proben  an  bestimmten  Stellen  angefertigt. 


Im  gesamten  Testprogramm  zur  Untersuchung  des  Einflusses  von  Tem- 
peratur  und  Feuchtigkeit  auf  schwingbeanspruchte  CFK-Laminate 
wurden  ca.  600  Proben  getestet,  wovon  420  im  Zugschwellversuch 
bis  zum  Versagen  beansprucht  wurden,  ca.  150  mit  NDT-Methoden  auf 
den  Schadigungszustand  untersucht  wurden,  und  ca.  30  in  stati- 
schen  Vorversuchen  gemessen  wurden.  In  Tabelle  I  ist  die  Zahl  der 
untersuchten  Proben  fur  jeden  Laminataufbau  aufgelistet. 
Insgesamt  sind  es  ca.  130  Proben  mehr  als  im  ESA-Proposal  vorge- 
schlagen  wurden,  da  sich  die  Notwendigkei t  ergab,  auf  Grund  der 
zwi schenzei tl ich  gewonnenen  Erkenntnisse  und  Ergebnisse  vor  al- 
lem  unter  normalen  Umgebungsbedingungen  groBere  Probenumf ange  zu 
untersuchen,  da  sonst  eine  sinnvolle  statistische  Auswertung  in 
Frage  gesteilt  war. 


3.  VERSUCHSERGEBNISSE 

3.1  Versuche  bei  normalen  Umgebungsbeding’ung’en  (ad) 

Obwohl  in  der  Literatur  schon  umfangreiche  Ergebnisse  von 
Schwingfestigkeitsuntersuchungen  an  ungekerbten  CFK-Laminaten 
mit  verschiedenen  Schichtfolgen  vorliegen  jl],  [2],  |3],  wurden 
diese  Versuche  in  dem  oben  beschriebenen  Umfang  durchgefiihrt ,  urn 
einen  Vergleich  zwischen  Proben,  die  unter  absolut  gleichen  Be- 
dingungen  hergestellt  wurden,  zu  ermoglichen,  da  erf ahrungsgemaB 
der  EinfluB  des  Herstellungsprozesses  nicht  unerheblich  ist. 
Beziiglich  des  Einflusses  der  Orientierung  der  Einzellagen  und  des 
Laminataufbaus  auf  die  Schwingf estigkeit  bzw.  die  Lebensdauer 
wurde  hierbei  gute  Ubereinstimmung  mit  an  anderer  Stelle  erziel- 
ten  Ergebnissen  festgestellt . 

Wie  empfindlich  die  Schwingf estigkeit  im  Zugbereich  schon  durch 
wenige  Grad  "of f-axis-Orientierung"  der  Lagen  beeinfluBt  wird, 
ist  bei  dem  Vergleich  zwischen  den  on-axis  0°-Proben  und  den  off- 
axis-Proben  mit  5°,  15°  und  45°  Winkel  zur  Belastungsrichtung  zu 
erkennen,  Bild  1.  In  diesem  Bild  sind  nur  die  Lebensdauer linien 
fur  50%  Uberlebenswahrscheinlichkeit  aufgetragen,  urn  die  Dar- 
stellung  nicht  zu  uniabersichtlich  zu  machen.  Die  Streubreiten  und 


die  Linien  fur  10%  und  90%  Oberlebenswahrscheinlichkei t  sind  je¬ 
wel  Is  aus  den  einzelnen  Diagrammen  zu  entnehmen,  in  diesem  Fall 
den  Bildern  6,  9,  10  und  14. 

In  ahnlicher  Weise  fallt  der  Vergleich  der  Schwingfestigkeiten 
fur  die  drei  untersuchten  angle-ply-Laminate  aus,  der  in  Bild  2 
dargestellt  ist.  Auch  hier  fallt  die  Schwingf estigkei t  vom  mehr 
f aserkontrollierten  ±15°-Laminat  zum  mehr  matrixkontrollierten  ± 
45°-Laminat  erheblich  ab,  wobei  jedoch  auffallt,  daB  die  Neigung 
der  Wohlerlinie  des  ±15°-Laminats  etwas  groBer  ist  als  bei  den 
anderen.  Auch  hier  sind  aus  Griinden  der  Ubersichtlichkei t  die 
Streubreiten  und  10%-  bzw.  90%-Linien  weggelassen,  diese  sind 
wieder  in  den  Einzeldarstel lungen  in  den  Bildern  15,  19  und  23 
enthalten. 

3.2  EinfluG  von  Temperatur  und  Feuchte  auf  unidirektionale  Laminate 

Die  Auswirkung  von  Temperatur  und  Feuchte  auf  unidirektionale  0°- 
Laminate  ist  in  Bild  3  aufgetragen.  Der  Ubersichtlichkei t  halber 
sind  auch  hier  nur  die  50%-Lebensdauerlinien  ei.  gezeichnet .  Man 
erkennt,  daB  die  Temperatur  offenbar  auf  die  Neigung  der  Wohler- 
linie  einen  EinfluB  hat,  wohingegen  die  Feuchtigkeit  sogar  zu 
einer  giinstigen  Verschiebung  der  Mittelwerte  fuhrt.  Dieses  ist 
sicher  auf  den  Ausgleich  der  durch  den  Curing-ProzeB  eingebrach- 
ten  Vorspannungen  zuriickzufiihren.  Allerdings  ist  auch  die 
Streuung  der  Bruchlastspielzahlen  bei  Feuchtigkeit  erheblich 
hoher,  sodaB  vermutlich  noch  mehr  Versuche  durchgefiihrt  warden 
miissen,  urn  zu  eindeutigen  und  statistisch  hinreichend  gesicher- 
ten  Aussagen  zu  gelangen. 

Besser  als  bei  der  Darstellung  der  Lebensdauerlinien  wird  deut- 
lich,  wie  unterschiedlich  Mi ttelwertverschiebung  und  Streuung 
auftreten,  wenn  z.B.  beim  Vergleich  zwischen  normaler  und  hoherer 
Temperatur  die  Auswertung  der  einzelnen  Spannungshorizonte  im 
Wahrscheinlichkei tsnetz  der  zwei-parametrigen  Weibull-Vertei- 
lung  aufgetragen  ist.  Diese  ist  fur  die  Spannungen  1200  N/mm*, 
1100  N/mm*  und  1000  N/mm*  in  den  Bildern  35,  36  und  37 
dargestellt.  Noch  deutlicher  wird  dieses  bei  Hinzunahme  der 
Feuchtigkeit,  wie  fur  1200  N/mml  in  Bild  38  zu  sehen  ist.  Hieraus 


ergibt  sich,  daB  fiir  eine  verniinftige  statistische  Aussage  bei 
vorliegenden  betrachtlichen  Streuungen  mehr  Versuche  pro  Span- 
nungsniveau  durchgefiihrt  werden  miissen .  Dieses  fiihrt  dann  zu 
einer  Darstellung  im  Weibull-Netz  mit  Angabe  des  Vertrauensbe- 
reichs  fiir  die  jeweils  erforderliche  Vertrauens- 
wahrscheinlichkeit ,  wie  an  dem  Beispiel  in  Bild  39  gesehen  werden 
kann. 

Eine  Auswertung  der  Lebensdauerlinien  fiir  0°-Laminate  bei  Tempe- 
ratur  und  Feuchte  war  nicht  moglich,  da  die  Bruchlastspielzahlen 
fast  alle  iiber  1200000  lagen.  Bei  Erreichen  dieser  Lastspielzah- 
len  wurden  bei  alien  Versuchen  die  Proben  ausgespannt,  um  die 
Gesamtversuchszei t  nicht  unnotig  zu  vergroBern.  Die  Ergebnisse 
zeigen  aber,  dafi  Temperatur  und  Feuchte  gemeinsam,  ebenfalls  wie 
nur  Feuchte  allein,  einen  giinstigen  EinfluB  auf  die  Schwingfe- 
stigkeit  ausiiben. 

Die  Auswirkung  von  Temperatur  und  Feuchte  auf  die  Schwingf estig- 
keit  unidi rektionaler  of f-axis-Laminate  wurde  nur  an  dem  15°- 
Laminat  untersucht.  Sowohl  Temperatur  als  auch  Feuchte  fiihren  zu 
einer  Abminderung  der  Schwingfestig-kei  t,  wie  in  Bild  10  und  Bild 
11  zu  erkennen  ist.  Fiir  das  Spannungsniveau  140  N/mm1  ist  dieses 
auch  sehr  gut  aus  der  Auftragung  im  Weibull-Netz  zu  sehen,  Bild 
40.  Beide  gemeinsam  reduzieren  betrachtlich  die  Schwingf estig- 
keit,  bewirken  aber  auch  ein  Verringerung  der  Neigung  der 
Wohlerlinie,  siehe  Bild  12. 

3.3  EinfluB  von  Temperatur  und  Feuchte  auf  angle -ply- Laminate 

Bild  4  zeigt  den  Vergleich  der  Wohlerlinien  in  Abhangigkeit  von 
den  Umgebungsbedingungen  fiir  die  drei  angle-ply-Laminate .  Deut- 
lich  ist  hier  fiir  die  ±15°-  und  ±30°-  Proben  eine  signifikante 
Abminderung  der  Schwingf estigkeit  zu  erkennen,  die  besonders  bei 
gemeinsamer  Einwirkung  von  Temperatur  und  Feuchtigkeit  auftritt. 
Auch  bei  diesen  Laminaten  treten  auf  einzelnen  Spannungshorizon- 
ten  sehr  groBe  Unterschiede  in  den  Streuungen  auf,  wie  man  in  den 
Bildern  41,  42  und  43  fur  alle  drei  Laminattypen  am  Vergleich  der 
Weibullvertei lungen  sehen  kann. 

Ein  moglicherwei se  vorhandener  EinfluB  der  Probenherstellung  auf 
die  Streuung  ist  nicht  auszuschlieBen,  da  haufig  festgestellt 


V 
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wurde,  daB  bei  Proben  gleichen  Laminataufbaus,  die  aus  verschie- 
denen  Platten  stammten,  typische  Unterschiede  in  den  Bruchlast- 
spielzahlen  auftraten.  Um  hier  mehr  Klarheit  zu  gewinnen,  miiBten 
gezielte  statistische  Untersuchungen  mit  groBeren  Stichprobenum- 
fangen  durchgefiihrt  werden.  Bild  44  und  45  zeigen,  daB  bei  einem 
Stichprobenumf ang  von  nur  3  oder  4  Proben  der  Vertrauensbereich 
fur  die  Verteilung  so  groB  wird,  daB  eine  sinnvolle  Interpreta¬ 
tion  fragwurdig  wird. 


3.4  EinfluG  von  Tempera tur  und  Feuchte  auf  multidirektionale  Laminate 

In  Bild  5  ist  die  Einwirkung  der  Umgebungsbedingungen  auf  die 

Schwingf estigkei t  der  beiden  multidirektionalen  Laminate  mit  den 

Laminataufbauten  [ 02O/+45o/02°/-45o/0o/90o } g  (in  den  Rechner- 

plots  mit  Type  1  bezeichnet)  und  ( 0°/+45o/-45o/90°  ]  (in  den 

s 

Plots  mit  Type  3  bezeichnet)  durch  Vergleich  der  50%-Linien  auf- 
getragen.  Fur  900  N/mm2  und  800  N/mm*  bei  Typ  1  und  fur  445  N/mm2 
bei  Typ  3  ist  zusatzlich  der  Vergleich  der  Verteilungen  im  Wei- 
bull-Netz  in  den  Bildern  46,  47  und  48  zu  sehen. 

Bei  beiden  Laminattypen  ist  eine  leichte  Verbesserung  der 
Schwingfestigkeit  durch  die  Einwirkung  der  feuchten  Umgebung  bei 
Norma  1 t«mperatur  f estzuste 1 len ,  sowie  eine  leichte  Abminderung 
bei  hoherer  Temperatur,  die  bei  dem  [ 0°/+45o/-45o/90° ]  -Laminat 
mit  nur  zwei  0°-Lagen  ausgepragter  ist.  Bemerkenswert  ist,  daB 
bei  Tempt -atur  und  Feuchtigkeit  zusammen  bei  beiden  Laminaten 
zvar  die  gemessenen  Mittelwerte  etwas  niedriger  liegen  und  die 
S^reuungen  uberwiegend  groBer  sind,  die  flache  Neigung  der  Woh- 
lerlimen  jedoch  erhalten  bleibt.  Man  kann  daher  nicht  eine 
signifikante  Degradation  der  Schwingfestigkeit  feststellen. 
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4.  ZUSAMMENFASSUNG 


Die  zusammenf assende  Betrachtung  der  bisher  gewonnenen  Ergebnis- 
se  fiihrt  zu  dem  SchluB,  daB  ein  EinfluB  durch  die  Umgebungstempe- 
ratur  und  - f euchtigkeit  bei  alien  Laminattypen  vorhanden  ist, 
wobei  wegen  der  z.T.  groBen  Streuungen  eindeutige  Abhangigkei ten 
und  Zusammenhange  nur  schwer  feststellbar  sind.  Genauere  Aussa- 
gen  konnen  hier  statistische  Signif ikanztests  liefern,  die 
jedoch  im  Mittel  groBere  Stichprobenumf ange  erfordern  und  daher 
nur  zum  Teil  durchgefuhrt  werden  konnten.  Weiteren  AufschluB 
konnten  rastere lektronenmikroskopi sche  Auswertungen  bringen, 
deren  bisheriger  Stand  jedoch  noch  nicht  zu  eindeutig  verwendba- 
ren  Unterscheidungsmerkmalen  gefiihrt  hat. 

Anhand  der  Lebensdauerlinien  und  der  Betrachtung  ihrer  Neigung 
laBt  sich  jedoch  feststellen,  daB  bei  mehr  matrixkontrol lierten 
Laminaten  wie  den  angle-ply-Laminaten  Temperatur  und  Feuchte 
einen  abmindernden  EinfluB  auf  die  Schwingf estigkei t  haben,  wah- 
rend  bei  mehr  f aserkontrollierten  wie  den  multidirektionalen  und 
den  unidirektionalen  on-axis-Laminaten  die  Feuchtigkeit  sogar  zu 
einer  Verbesserung  fiihren  kann.  Zur  gleichzeitigen  Einwirkung 
von  Temperatur  und  Feuchte  ist  zu  bemerken,  daB  zwar  teilweise 
eine  beachtliche  Erhohung  der  Streuung  auftritt,  die  flache  Nei¬ 
gung  der  Wohlerlinien  bei  mehr  f aserkontrollierten  Laminaten 
jedoch  nur  unwesentlich  geandert  wird. 

Das  bedeutet,  daB  auch  unter  "worst  case"  Umgebungsbedingungen 
die  Extrapolation  der  Lebensdauerlinie  fur  90%  Oberlebenswahr- 
scheinl ichkeit  in  den  Bereich  der  Dauerf estigkeit  im  Vergleich  zu 
herkommlichen  metallischen  Werkstoffen  hervorragende  Werte  der 
Zugschwingfestigkei t  liefert. 


Unidirektional : 


On-axis  0°  8 
Off-axis  5°  8 
Off-axis  15°  8 
Off-axis  45°  8 


Multidirektional : 

[02o/+45o/02o/-45o/0°/90o]s  16 

[ 0°/+45  °/-45 °/90° ]  7 

s 


Angle-ply: 


[±15°]8 
[  ±30°  1 8 
[  145°  ]g 


v 


fc* 


Tabelle  I:  Anzahl  der  untersuchten  Proben  pro  Laminataufbau 


Probability  of  Survival 


unidirektionaler  Laminate 


Probability  of  Survival 


50%  -  Probab 


Probability  of  Survival 


90%,  50%,  1 0%  probob 


10%  probability  of  survival 


90%,  50%,  10%  probab 
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0%  proocbility  of  survival 


90%,  50%,  10%  probability  of  survival 


Probability  of  failure, 
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131/13 


Two-parameter  Weibull  distribution 


Bi Id  35 


Probability  of  failure, 


Two-parameter  Weibull  distribution 


Bi Id  36 


Probability  of  failure, 
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Two-parameter  Weibull  distribution 


Bild  37 


ility  of  failure, 


Unidirectional  laminate,  on-axis 
Tension  — tension  fatigue,  R  =  0.1,  o  =  1100  N/mm2 
Temp.  =  22°  C,  rel.  hum.  =  40% 


Two-parameter  Weibull  distribution 

131/13 
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Confidence  interval  =  90% 

Bild  39 
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Probability  of  failure 


Multidirectional  laminate,  type  1 
Tension-tension  fatigue,  R  =  0. 1,  a  --  900  N/mm 
a  :  Temp.  =  22°  C,  rel.  hum.  =  40% 

o  :  Temp.  =  80°  C,  rel.  hum.  =  05% 

a  :  Temp.  =  28°  C,  rel.  hum.  =  95% 

o  :  Temp.  =  80°  C,  rel.  hum.  =  95% 
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Bild  46 
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Probability  of  failure, 
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